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Chapter 1 


Introduction 


1.1 Origin of word: Algorithm 

The word Algorithm comes from the name of the muslim author Abu Ja Jar Mohammad ibn Musa 
al-Khowarizmi. He was born in the eighth century at Khwarizm (Kheva), a town south of river Oxus in 
present Uzbekistan. Uzbekistan, a Muslim country for over a thousand years, was taken over by the 
Russians in 1873. 

His year of birth is not known exactly. Al-Khwarizmi parents migrated to a place south of Baghdad when 
he was a child. It has been established from his contributions that he flourished under Khalifah 
Al-Mamun at Baghdad during 813 to 833 C.E. Al-Khwarizmi died around 840 C.E. 

Much of al- Khwarizmi’s work was written in a book titled a l Kitab al-mukhatasar fi hisab al-jabr 
wa’l-muqabalah (The Compendious Book on Calculation by Completion and Balancing). It is from the 
titles of these writings and his name that the words algebra and algorithm are derived. As a result of his 
work, al-Khwarizmi is regarded as the most outstanding mathematician of his time 


1.2 Algorithm: Informal Definition 

An algorithm is any well-defined computational procedure that takes some values, or set of values, as 
input and produces some value, or set of values, as output. An algorithm is thus a sequence of 
computational steps that transform the input into output. 


1.3 Algorithms, Programming 

A good understanding of algorithms is essential for a good understanding of the most basic element of 
computer science: programming. Unlike a program, an algorithm is a mathematical entity, which is 
independent of a specific programming language, machine, or compiler. Thus, in some sense, algorithm 
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design is all about the mathematical theory behind the design of good programs. 

Why study algorithm design? There are many facets to good program design. Good algorithm design is 
one of them (and an important one). To be really complete algorithm designer, it is important to be aware 
of programming and machine issues as well. In any important programming project there are two major 
types of issues, macro issues and micro issues. 

Macro issues involve elements such as how does one coordinate the efforts of many programmers 
working on a single piece of software, and how does one establish that a complex programming system 
satisfies its various requirements. These macro issues are the primary subject of courses on software 
engineering. 

A great deal of the programming effort on most complex software systems consists of elements whose 
programming is fairly mundane (input and output, data conversion, error checking, report generation). 
However, there is often a small critical portion of the software, which may involve only tens to hundreds 
of lines of code, but where the great majority of computational time is spent. (Or as the old adage goes: 
80% of the execution time takes place in 20% of the code.) The micro issues in programming involve 
how best to deal with these small critical sections. 

It may be very important for the success of the overall project that these sections of code be written in the 
most efficient manner possible. An unfortunately common approach to this problem is to first design an 
inefficient algorithm and data structure to solve the problem, and then take this poor design and attempt 
to fine-tune its performance by applying clever coding tricks or by implementing it on the most expensive 
and fastest machines around to boost performance as much as possible. The problem is that if the 
underlying design is bad, then often no amount of fine-tuning is going to make a substantial difference. 

Before you implement, first be sure you have a good design. This course is all about how to design good 
algorithms. Because the lesson cannot be taught in just one course, there are a number of companion 
courses that are important as well. CS301 deals with how to design good data structures. This is not 
really an independent issue, because most of the fastest algorithms are fast because they use fast data 
structures, and vice versa. In fact, many of the courses in the computer science program deal with 
efficient algorithms and data structures, but just as they apply to various applications: compilers, 
operating systems, databases, artificial intelligence, computer graphics and vision, etc. Thus, a good 
understanding of algorithm design is a central element to a good understanding of computer science and 
good programming. 


1.4 Implementation Issues 

One of the elements that we will focus on in this course is to try to study algorithms as pure mathematical 
objects, and so ignore issues such as programming language, machine, and operating system. This has 
the advantage of clearing away the messy details that affect implementation. But these details may be 
very important. 

For example, an important fact of current processor technology is that of locality of reference. Frequently 
accessed data can be stored in registers or cache memory. Our mathematical analysis will usually ignore 
these issues. But a good algorithm designer can work within the realm of mathematics, but still keep an 
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open eye to implementation issues down the line that will be important for final implementation. For 
example, we will study three fast sorting algorithms this semester, heap-sort, merge-sort, and quick-sort. 
From our mathematical analysis, all have equal running times. However, among the three (barring any 
extra considerations) quick sort is the fastest on virtually all modem machines. Why? It is the best from 
the perspective of locality of reference. However, the difference is typically small (perhaps 10-20% 
difference in running time). 

Thus this course is not the last word in good program design, and in fact it is perhaps more accurately just 
the first word in good program design. The overall strategy that I would suggest to any programming 
would be to first come up with a few good designs from a mathematical and algorithmic perspective. 

Next prune this selection by consideration of practical matters (like locality of reference). Finally 
prototype (that is, do test implementations) a few of the best designs and run them on data sets that will 
arise in your application for the final fine-tuning. Also, be sure to use whatever development tools that 
you have, such as profilers (programs which pin-point the sections of the code that are responsible for 
most of the running time). 


1.5 Course in Review 

This course will consist of four major sections. The first is on the mathematical tools necessary for the 
analysis of algorithms. This will focus on asymptotics, summations, recurrences. The second element 
will deal with one particularly important algorithmic problem: sorting a list of numbers. We will show a 
number of different strategies for sorting, and use this problem as a case-study in different techniques for 
designing and analyzing algorithms. 

The final third of the course will deal with a collection of various algorithmic problems and solution 
techniques. Finally we will close this last third with a very brief introduction to the theory of 
NP-completeness. NP-complete problem are those for which no efficient algorithms are known, but no 
one knows for sure whether efficient solutions might exist. 


1.6 Analyzing Algorithms 

In order to design good algorithms, we must first agree the criteria for measuring algorithms. The 
emphasis in this course will be on the design of efficient algorithm, and hence we will measure 
algorithms in terms of the amount of computational resources that the algorithm requires. These 
resources include mostly running time and memory. Depending on the application, there may be other 
elements that are taken into account, such as the number disk accesses in a database program or the 
communication bandwidth in a networking application. 

In practice there are many issues that need to be considered in the design algorithms. These include 
issues such as the ease of debugging and maintaining the final software through its life-cycle. Also, one 
of the luxuries we will have in this course is to be able to assume that we are given a clean, fully-specified 
mathematical description of the computational problem. In practice, this is often not the case, and the 
algorithm must be designed subject to only partial knowledge of the final specifications. Thus, in practice 
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it is often necessary to design algorithms that are simple, and easily modified if problem parameters and 
specifications are slightly modified. Fortunately, most of the algorithms that we will discuss in this class 
are quite simple, and are easy to modify subject to small problem variations. 


1.7 Model of Computation 

Another goal that we will have in this course is that our analysis be as independent as possible of the 
variations in machine, operating system, compiler, or programming language. Unlike programs, 
algorithms to be understood primarily by people (i.e. programmers) and not machines. Thus gives us 
quite a bit of flexibility in how we present our algorithms, and many low-level details may be omitted 
(since it will be the job of the programmer who implements the algorithm to fill them in). 

But, in order to say anything meaningful about our algorithms, it will be important for us to settle on a 
mathematical model of computation. Ideally this model should be a reasonable abstraction of a standard 
generic single -processor machine. We call this model a random access machine or RAM. 

A RAM is an idealized machine with an infinitely large random-access memory. Instructions are 
executed one-by-one (there is no parallelism). Each instruction involves performing some basic operation 
on two values in the machines memory (which might be characters or integers; let’s avoid floating point 
for now). Basic operations include things like assigning a value to a variable, computing any basic 
arithmetic operation (+, - , x , integer division) on integer values of any size, performing any comparison 
(e.g. x < 5) or boolean operations, accessing an element of an array (e.g. A[10]). We assume that each 
basic operation takes the same constant time to execute. 

This model seems to go a good job of describing the computational power of most modern (nonparallel) 
machines. It does not model some elements, such as efficiency due to locality of reference, as described 
in the previous lecture. There are some “loop-holes” (or hidden ways of subverting the rules) to beware 
of. For example, the model would allow you to add two numbers that contain a billion digits in constant 
time. Thus, it is theoretically possible to derive nonsensical results in the form of efficient RAM 
programs that cannot be implemented efficiently on any machine. Nonetheless, the RAM model seems to 
be fairly sound, and has done a good job of modelling typical machine technology since the early 60’s. 


1.8 Example: 2-dimension maxima 

Let us do an example that illustrates how we analyze algorithms. Suppose you want to buy a car. You 
want the pick the fastest car. But fast cars are expensive; you want the cheapest. You cannot decide which 
is more important: speed or price. Definitely do not want a car if there is another that is both faster and 
cheaper. We say that the fast, cheap car dominates the slow, expensive car relative to your selection 
criteria. So, given a collection of cars, we want to list those cars that are not dominated by any other. 

Here is how we might model this as a formal problem. 

• Let a point p in 2-dimensional space be given by its integer coordinates, p = (p.x,p.y). 
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• A point p is said to be dominated by point q if p.x < q.x and p.y < 9-y- 

• Given a set of n points, P = {pi , Vi, ■ ■ ■ , Pn} in 2-space a point is said to be maximal if it is not 
dominated by any other point in P. 

The car selection problem can be modelled this way: For each car we associate (x, y ) pair where x is the 
speed of the car and y is the negation of the price. High y value means a cheap car and low y means 
expensive car. Think of y as the money left in your pocket after you have paid for the car. Maximal 
points correspond to the fastest and cheapest cars. 

The 2-dimensional Maxima is thus defined as 

• Given a set of points P = {pi , V 2 , ■ ■ ■ , Pn} in 2-space, output the set of maximal points of P, i.e., 
those points p L such that pi is not dominated by any other point of P. 

Here is set of maximal points for a given set of points in 2-d. 



n i i i i i i i 

2 4 6 8 10 12 14 16 18 


Figure 1.1: Maximal points in 2-d 


Our description of the problem is at a fairly mathematical level. We have intentionally not discussed how 
the points are represented. We have not discussed any input or output formats. We have avoided 
programming and other software issues. 


1.9 Brute-Force Algorithm 

To get the ball rolling, let’s just consider a simple brute-force algorithm, with no thought to efficiency. 

Let P = {pi , P2, • • ■ > Pn} be the initial set of points. For each point p t , test it against all other points p j. If 
Pi is not dominated by any other point, then output it. 
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This English description is clear enough that any (competent) programmer should be able to implement 
it. However, if you want to be a bit more formal, it could be written in pseudocode as follows: 

MAXiMA(int n, Point P[l...n]) 

1 for i <— 1 to n 

2 do maximal <— true 

3 for j <— 1 to n 

4 do 

5 if (i j) and (P[i].x < P[j].x) and (P[i].y < P[j].y) 

6 then maximal <— false; break 

7 if (maximal = true) 

8 then output P[i] 


There are no formal rules to the syntax of this pseudo code. In particular, do not assume that more detail 
is better. For example, I omitted type specifications for the procedure Maxima and the variable maximal, 
and I never defined what a Point data type is, since I felt that these are pretty clear from context or just 
unimportant details. Of course, the appropriate level of detail is a judgement call. Remember, algorithms 
are to be read by people, and so the level of detail depends on your intended audience. When writing 
pseudo code, you should omit details that detract from the main ideas of the algorithm, and just go with 
the essentials. 

You might also notice that I did not insert any checking for consistency. For example, I assumed that the 
points in P are all distinct. If there is a duplicate point then the algorithm may fail to output even a single 
point. (Can you see why?) Again, these are important considerations for implementation, but we will 
often omit error checking because we want to see the algorithm in its simplest form. 

Here are a series of figures that illustrate point domination. 



2 4 6 8 10 12 14 16 18 



2 4 6 8 10 12 14 16 18 


Figure F2: Points that dominate (4,11) 


Figure 1.3: Points that dominate (9, 10) 
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Figure 1.4: Points that dominate (7, 7) Figure 1.5: The maximal points 


1.10 Running Time Analysis 

The main purpose of our mathematical analysis will be measuring the execution time. We will also be 
concerned about the space (memory) required by the algorithm. 

The running time of an implementation of the algorithm would depend upon the speed of the computer, 
programming language, optimization by the compiler etc. Although important, we will ignore these 
technological issues in our analysis. 

To measure the running time of the brute-force 2-d maxima algorithm, we could count the number of 
steps of the pseudo code that are executed or, count the number of times an element of P is accessed or, 
the number of comparisons that are performed. 

The running time depends upon the input size, e.g. n Different inputs of the same size may result in 
different running time. For example, breaking out of the inner loop in the brute-force algorithm depends 
not only on the input size of P but also the structure of the input. 

Two criteria for measuring running time are worst-case time and average-case time. 


Worst-case time is the maximum running time over all (legal) inputs of size rt. Let I denote an input 
instance, let |I[ denote its length, and let T(I) denote the running time of the algorithm on input I. 
Then 


Tworst(n) = max T(I) 
|I|=n 


Average-case time is the average running time over all inputs of size rt. Let p(I) denote the probability 
of seeing this input. The average-case time is the weighted sum of running times with weights 
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being the probabilities: 

Tov fl (n) = ^p(I)T(I) 

|I|=n 


We will almost always work with worst-case time. Average-case time is more difficult to compute; it is 
difficult to specify probability distribution on inputs. Worst-case time will specify an upper limit on the 
running time. 


1.10.1 Analysis of the brute-force maxima algorithm. 


Assume that the input size is n, and for the running time we will count the number of time that any 
element of P is accessed. Clearly we go through the outer loop n times, and for each time through this 
loop, we go through the inner loop n times as well. The condition in the if-statement makes four accesses 
to P. The output statement makes two accesses for each point that is output. In the worst case every point 
is maximal (can you see how to generate such an example?) so these two access are made for each time 
through the outer loop. 


MAXlMA(int n, Point P[1 . . . n]) 

1 for i <— 1 to n [ >• A,:- C] 

2 do maximal <— true 

3 for j <— 1 to n [ ;■ A/,-- - j 

4 do 

5 if (i ^ j)&(P[i].x < P[j].x)&(P[i].y < P[j].y) 

6 then maximal <— false break 

7 if maximal 

8 then output P[i].x, P[i].y ' v, •.> < > I 



Thus we might express the worst-case running time as a pair of nested summations, one for the i-loop 
and the other for the j-loop: 


T(n)=Jj2 + ^4 


i=l 

n 


j=1 


}j4n + 2) 


i=1 


= (4n + 2)n = 4n 2 + 2u 


For small values of n, any algorithm is fast enough. What happens when n gets large? Running time 
does become an issue. When n is large, u 2 term will be much larger than the n term and will dominate 
the running time. 

We will say that the worst-case running time is 0(u 2 ). This is called the asymptotic growth rate of the 
function. We will discuss this 0-notation more formally later. 
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The analysis involved computing a summation. Summation should be familiar but let us review a bit 
here. Given a finite sequence of values cii , < 12 , . . . , a n , their sum ai + a 2 + • • • + a n is expressed in 
summation notation as 


5 > 


i=1 


If n = 0, then the sum is additive identity, 0. 
Some facts about summation: If c is a constant 


^ca t = c^ 


a. 


i=l 


i=1 


and 


Y (cn + bt) = Y aj+Y bj 

i=t i=l i=1 

Some important summations that should be committed to memory. 

Arithmetic series 


y~ 1 — 1 + 2 + . . . + 


n 


i=l 


nfn + 1 


= 0(ri 2 


Quadratic series 


Y_ i 2 = 1 + 4 + 9 + ...+ 
2n 3 + 3n 2 + n 


n 


i=l 


= 0(n d 


Geometric series 


x^ — 1 -)- x T x 2 T • . . T 
x (n+1) _ -| 


X 


i=1 


X — 1 


= 0(u 2 


If 0 < x < 1 then this is 0(1 ), and if x > 1 , then this is 0(x n ). 

Harmonic series For n > 0 

n 1 


i=l 


, 11 1 , 

= 1+ :r + - + ... + -^lnn 
2 3 n 

= 0(ln n) 



16 


CHAPTER 1. INTRODUCTION 


1.11 Analysis: A Harder Example 

Let us consider a harder example. 

NESTED-LOOPS () 

1 for i 1 to n 

2 do 

3 for j <— 1 to 2i 

4 do k = j . . . 

5 while (k > 0) 

6 dok = k-l ... 


How do we analyze the running time of an algorithm that has complex nested loop? The answer is we 
write out the loops as summations and then solve the summations. To convert loops into summations, we 
work from inside-out. 

Consider the inner most while loop. 


NESTED-LOOPS () 

1 for i <— 1 to n 

2 do for j <— 1 to 2i 

3 do k = j 

4 while (k > 0) JT| 

5 do k = k - 1 


It is executed for k = j, j — 1 , j — 2, . . . , 0. Time spent inside the while loop is constant. Let I() be the 
time spent in the while loop. Thus 

j 

I(i) = H 1 =3 + 1 

k=0 


Consider the middle for loop. 


NESTED-LOOPS () 

1 for i <— 1 to n 

2 do for j <— 1 to 2i |*T] 

3 do k = j 

4 while (k > 0) 

5 do k = k — 1 
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Its running time is determined by i. Let M.Q be the time spent in the for loop: 

2i 

j=i 

2i 

= £(j + i) 

j=1 

2i 2i 

= L> + L 1 

j=i j=i 
2i(2i + 1 ; 


2 

= 2i 2 + 3i 


2i 


Finally the outer-most for loop. 

NESTED-LOOPS () 

1 for i <— 1 to n 

2 do for ) <— 1 to 2i 

3 do k = j 

4 while (k > 0) 

5 do k = k - 1 


Let T( ) be running time of the entire algorithm: 

rt 

T(n) = Y_ M(i) 

i=1 

rt 

= ^(2i 2 + 3i) 

i=1 

n n 

= ^ 2i 2 + ^ 3i 

i=l i=l 

_ n 2n 3 + 3n 2 + n Tt(n+1) 
~ 2 6 + 3 2 
4n 3 + 1 5n 2 + 1 1 n 
6 

- 0(n 3 ) 


1.11.1 2-dimension Maxima Revisited 

Recall the 2-d maxima problem: Let a point p in 2-dimensional space be given by its integer coordinates, 
p = (p.x, p.p). A point p is said to dominated by point q if p.x < q.x and p.q < q.q. Given a set of n 
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points, P = {pi , P 2 , • • • , Pn} in 2-space a point is said to be maximal if it is not dominated by any other 
point in P. The problem is to output all the maximal points of P. We introduced a brute-force algorithm 
that ran in 0(n 2 ) time. It operated by comparing all pairs of points. Is there an approach that is 
significantly better? 

The problem with the brute-force algorithm is that it uses no intelligence in pruning out decisions. For 
example, once we know that a point pt is dominated by another point pj, we do not need to use pt for 
eliminating other points. This follows from the fact that dominance relation is transitive. If pj dominates 
Pi and pi dominates then pj also dominates p H ; Pi is not needed. 


1.11.2 Plane-sweep Algorithm 


The question is whether we can make an significant improvement in the running time? Here is an idea for 
how we might do it. We will sweep a vertical line across the plane from left to right. As we sweep this 
line, we will build a structure holding the maximal points lying to the left of the sweep line. When the 
sweep line reaches the rightmost point of P , then we will have constructed the complete set of maxima. 
This approach of solving geometric problems by sweeping a line across the plane is called plane sweep. 

Although we would like to think of this as a continuous process, we need some way to perform the plane 
sweep in discrete steps. To do this, we will begin by sorting the points in increasing order of their 
x-coordinates. For simplicity, let us assume that no two points have the same y-coordinate. (This limiting 
assumption is actually easy to overcome, but it is good to work with the simpler version, and save the 
messy details for the actual implementation.) Then we will advance the sweep-line from point to point in 
n discrete steps. As we encounter each new point, we will update the current list of maximal points. 

We will sweep a vertical line across the 2-d plane from left to right. As we sweep, we will build a 
structure holding the maximal points lying to the left of the sweep line. When the sweep line reaches the 
rightmost point of P, we will have the complete set of maximal points. We will store the existing 
maximal points in a list The points that p L dominates will appear at the end of the list because points are 
sorted by x-coordinate. We will scan the list left to right. Every maximal point with y -coordinate less 
than pi will be eliminated from computation. We will add maximal points onto the end of a list and 
delete from the end of the list. We can thus use a stack to store the maximal points. The point at the top 
of the stack will have the highest x-coordinate. 

Here are a series of figures that illustrate the plane sweep. The figure also show the content of the stack. 
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Figure 1.6: Sweep line at (2, 5) 
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Figure 1.7: Sweep line at (3, 13) 



sweep line 


i 

14- 


-► 


12 

(3,13)1 

O (12,12) 


10 - 

- (4,11) 

m O 0(14,10) 




(9,10) 


8 


(7,7)0 0(15,7) 


6 - 

(2,5) 



4 - 

- 

(10,5)° 


2 - 


0(13,3) 

(4,11) 



° 

(3,13) 

1 

2 ^ 

16 8 10 12 14 16 18 

stack 


Figure 1.8: Sweep line at (4,11) 
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Figure 1.9: Sweep line at (5,1) 
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Figure 1.10: Sweep line at (7, 7) 


Figure 1.11: Sweep line at (9, 10) 
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Figure 1.12: Sweep line at (10,5) 
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Figure 1.13: Sweep line at (12, 12) 
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Figure 1.14: The final maximal set 


Here is the algorithm. 

PL ANE- SWEEP- MAXIM A (n, P[l..ri]) 

1 sort P in increasing order by x; 

2 stack s; 

3 for i <— 1 to n 

4 do 

5 while (s.notEmpty () & s.topQ.p < m.y) 

6 do s.pop(); 

7 s.push(P[i]); 

8 output the contents of stack s; 
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1.11.3 Analysis of Plane-sweep Algorithm 

Sorting takes 0(nlogri); we will show this later when we discuss sorting. The for loop executes n 
times. The inner loop (seemingly) could be iterated (n — 1 ) times. It seems we still have an ri(u — 1 ) or 
0(n 2 ) algorithm. Got fooled by simple minded loop-counting. The while loop will not execute more n 
times over the entire course of the algorithm. Why is this? Observe that the total number of elements that 
can be pushed on the stack is n since we execute exactly one push each time during the outer for-loop. 

We pop an element off the stack each time we go through the inner while-loop. It is impossible to pop 
more elements than are ever pushed on the stack. Therefore, the inner while-loop cannot execute more 
than n times over the entire course of the algorithm. (Make sure that you understand this). 

The for-loop iterates n times and the inner while-loop also iterates n time for a total of 0(n). Combined 
with the sorting, the runtime of entire plane-sweep algorithm is 0(nlogrt). 

1.11.4 Comparison of Brute-force and Plane sweep algorithms 

How much of an improvement is plane-sweep over brute-force? Consider the ratio of running times: 

n 2 n 

nlogn log n 


n 

logn 

ri 

logn 

100 

7 

15 

1000 

10 

100 

10000 

13 

752 

100000 

17 

6021 

1000000 

20 

50171 


For n = 1 , 000, 000, if plane-sweep takes 1 second, the brute-force will take about 14 hours!. From this 
we get an idea about the importance of asymptotic analysis. It tells us which algorithm is better for large 
values of n. As we mentioned before, if n is not very large, then almost any algorithm will be fast. But 
efficient algorithm design is most important for large inputs, and the general rule of computing is that 
input sizes continue to grow until people can no longer tolerate the running times. Thus, by designing 
algorithms efficiently, you make it possible for the user to run large inputs in a reasonable amount of time. 
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Chapter 2 

Asymptotic Notation 


You may be asking that we continue to use the notation 0() but have never defined it. Let’s remedy this 
now. Given any function g(n), we define 0(g(n)) to be a set of functions that asymptotically equivalent 
to g(rt). Formally: 


0(g(n)) ={f(n) | there exist positive 

constants Ci , C 2 and n 0 such that 
0 < cig(n) < f(n) < c 2 g(n) 
for all n > n 0 } 

This is written as “f(n) G 0(g(n))” That is, f(n) and g(n) are asymptotically equivalent. This means 
that they have essentially the same growth rates for large n. For example, functions like 

• 4 n 2 , 

• (8n 2 + 2n — 3), 

• (u 2 /5 + s/n — lOlogn) 

• n(n — 3) 

are all asymptotically equivalent. As n becomes large, the dominant (fastest growing) term is some 
constant times n 2 . 

Consider the function 

f (n) = 8n 2 + 2n — 3 

Our informal rule of keeping the largest term and ignoring the constant suggests that f (n) G 0(n 2 ). Let’s 
see why this bears out formally. We need to show two things for 

f (n) = 8n 2 + 2n — 3 

Lower bound f (n) = 8n 2 + 2n — 3 grows asymptotically at least as fast as n 2 , 
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Upper bound f (n) grows no faster asymptotically than n 2 , 

Lower bound: f (n) grows asymptotically at least as fast as n 2 . For this, need to show that there exist 
positive constants Ci and n 0 , such that f (n) > c^n 2 for all n > n 0 . Consider the reasoning 

f (n) = 8n 2 + 2n — 3 > 8n 2 — 3 = 7n 2 + (n 2 — 3) > 7n 2 

Thus Ci = 7. We implicitly assumed that 2ri > 0 and n 2 — 3 > 0 These are not true for all n but if 
n > y/3, then both are true. So select rio > x/3. We then have f(n) > Ciri 2 for all n > no- 

Upper bound: f (n) grows asymptotically no faster than n 2 . For this, we need to show that there exist 
positive constants C 2 and rto, such that f (n) < C 2 n 2 for all n > rio. Consider the reasoning 

f (n) = 8n 2 + 2n — 3 < 8n 2 + 2 n < 8n 2 + 2n 2 = 1 0n 2 

Thus C 2 — 10. We implicitly made the assumption that 2n < 2n 2 . This is not true for all n but it is true 
for all n > 1 So select n 0 > 1 . We thus have f (n) < C 2 n 2 for all n > n 0 . 

From lower bound we have n 0 > i/3 From upper bound we have rto > 1 . Combining the two, we let rio 
be the larger of the two: n 0 > i/3. In conclusion, if we let Ci = 7, C 2 = 1 0 and n 0 > \/3, we have 

7n 2 < 8n 2 + 2n — 3 < 1 On 2 for all n > v/3 

We have thus established 

0 < Cig[n) < f(n) < c 2 g(n) foralln>n 0 
Here are plots of the three functions. Notice the bounds. 



Figure 2.1: Asymptotic Notation Example 


We have established that f (n) e n 2 . Let’s show why f (n) is not in some other asymptotic class. First, 
let’s show that f(n) ^ @(n). Show that f(n) ^ @(n). If this were true, we would have had to satisfy 



25 


both the upper and lower bounds. The lower bound is satisfied because f (n) = 8n 2 + 2n — 3 does grow 
at least as fast asymptotically as n. But the upper bound is false. Upper bounds requires that there exist 
positive constants C 2 and n 0 such that f (rt) < C 2 n for all n > n 0 . 

Informally we know that f (n) = 8n 2 + 2n — 3 will eventually exceed C 2 n no matter how large we make 
C 2 . To see this, suppose we assume that constants C 2 and n 0 did exist such that 8n 2 + 2n — 3 < C 2 n for 
all n > rto Since this is true for all sufficiently large n then it must be true in the limit as n tends to 
infinity. If we divide both sides by rt, we have 


lim ( 8n + 2 — — ] < C 2 . 

n— >00 y ny 

It is easy to see that in the limit, the left side tends to 00 . So, no matter how large C 2 is, the statement is 
violated. Thus f (n) ^ 0(n). 

Let’s show that f (rt) ^ 0(rt 3 ). The idea would be to show that the lower bound f (rt) > Ciri 3 for all 
n > no is violated, (ci and no are positive constants). Informally we know this to be true because any 
cubic function will overtake a quadratic. 

If we divide both sides by n 3 : 

(8 2 3 \ 

Inn — | > d 

n-^ooyn n z n A J 

The left side tends to 0. The only way to satisfy this is to set Ci = 0. But by hypothesis, Ci is positive. 
This means that f(n) ^ 0(n 3 ). 

The definition of 0-notation relies on proving both a lower and upper asymptotic bound. Sometimes we 
only interested in proving one bound or the other. The O-notation is used to state only the asymptotic 
upper bounds. 


0(g(n)) ={f(n) | there exist positive 

constants c and n 0 such that 
0 < f(n) < cg(n) 
for all n > rto} 

The 0-notation allows us to state only the asymptotic lower bounds. 

0(g(n)) = {f (n) | there exist positive 

constants c and n 0 such that 
0 < cg(n) < f (n) 
for all n > rt 0 } 
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The three notations: 


@(g(n)) : 0 < cig(n) < f(n) < c 2 g(n) 
0(g(n)) : 0 < f(n) < cg(n) 
a(g(n)) : 0 < cg(n) < f(n) 
for all n > n 0 


These definitions suggest an alternative way of showing the asymptotic behavior. We can use limits for 
define the asymptotic behavior. Limit rule for 0-notation: 


lim 

n— )oo 


f(n) 

g(n) 


= c, 


for some constant c > 0 (strictly positive but not infinity) then f (rt) G 0(g(n)). Similarly, the limit rule 
for O-notation is 


lim 

n— > oo 


f(rO 

g(n) 


= c, 


for some constant c > 0 (nonnegative but not infinite) then f(n) G 0(g(n)) and limit rule for 
0-notation: 


lim 

n-Aoo 


f(n) 

g(n) 


7 ^ 0 > 


(either a strictly positive constant or infinity) then f (n) G 0(g(n)). 


Here is a list of common asymptotic running times: 


• 0(1 ): Constant time; can’t beat it! 

• 0(log n): Inserting into a balanced binary tree; time to find an item in a sorted array of length ri 
using binary search. 

• 0(n): About the fastest that an algorithm can run. 

• 0(nlogn): Best sorting algorithms. 

• 0(n 2 ), 0(n 3 ): Polynomial time. These running times are acceptable when the exponent of n is 
small or n is not to large, e.g., n < 1 000. 

• 0(2 n ), 0(3 n ): Exponential time. Acceptable only if n is small, e.g., n < 50. 

• 0(n!), 0(ri n ): Acceptable only for really small ri, e.g. n < 20. 



Chapter 3 

Divide and Conquer Strategy 


The ancient Roman politicians understood an important principle of good algorithm design (although 
they were probably not thinking about algorithms at the time). You divide your enemies (by getting them 
to distrust each other) and then conquer them piece by piece. This is called divide-and-conquer. In 
algorithm design, the idea is to take a problem on a large input, break the input into smaller pieces, solve 
the problem on each of the small pieces, and then combine the piecewise solutions into a global solution. 
But once you have broken the problem into pieces, how do you solve these pieces? The answer is to 
apply divide-and-conquer to them, thus further breaking them down. The process ends when you are left 
with such tiny pieces remaining (e.g. one or two items) that it is trivial to solve them. Summarizing, the 
main elements to a divide-and-conquer solution are 

Divide: the problem into a small number of pieces 

Conquer: solve each piece by applying divide and conquer to it recursively 
Combine: the pieces together into a global solution. 


3.1 Merge Sort 

Divide and conquer strategy is applicable in a huge number of computational problems. The first 
example of divide and conquer algorithm we will discuss is a simple and efficient sorting procedure 
called We are given a sequence of n numbers A, which we will assume are stored in an array A[1 ..n]. 
The objective is to output a permutation of this sequence sorted in increasing order. This is normally 
done by permuting the elements within the array A. Here is how the merge sort algorithm works: 

• (Divide:) split A down the middle into two subsequences, each of size roughly n/2 

• (Conquer:) sort each subsequence by calling merge sort recursively on each. 

• (Combine:) merge the two sorted subsequences into a single sorted list. 
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Let’s design the algorithm top-down. We’ll assume that the procedure that merges two sorted list is 
available to us. We’ll implement it later. Because the algorithm is called recursively on sublists, in 
addition to passing in the array itself, we will pass in two indices, which indicate the first and last indices 
of the sub-array that we are to sort. The call MergeSort) A, p, r) will sort the sub-array A[p : r] and 
return the sorted result in the same sub-array.Here is the overview. If r = p, then this means that there is 
only one element to sort, and we may return immediately. Otherwise if (p ^ r) there are at least two 
elements, and we will invoke the divide-and-conquer. We find the index q, midway between p and r, 
namely q = (p + r)/2 (rounded down to the nearest integer). Then we split the array into sub-arrays 
A[p : q] and A[q + 1 : r]. Call MergeSort recursively to sort each sub-array. Finally, we invoke a 
procedure (which we have yet to write) which merges these two sub-arrays into a single sorted array. 

Here is the MergeSort algorithm. 

MERGE- SORT( array A, int p, int r) 

1 if (p < r) 

2 then 

3 q<-(p + r)/2 

4 MERGE-SORT(A,p, q) //sort A[p..q] 

5 MERGE-SORT(A, q + 1 , r) // sort A[q + 1 ..r] 

6 MERGE(A,p, q,r) // merge the two pieces 

The following figure illustrates the dividing (splitting) procedure. 


7 5 2 4 1 6 3 0 

7 5 2 4 

16 3 0 

7 5 2 4 

16 3 0 

7 5 2 4 

16 3 0 

split 


Figure 3.1: Merge sort divide phase 


Merging: All that is left is to describe the procedure that merges two sorted lists. Merge(A,p, q, r) 
assumes that the left sub-array, A[p : q], and the right sub-array, A[q + 1 : r], have already been sorted. 
We merge these two sub-arrays by copying the elements to a temporary working array called B. For 
convenience, we will assume that the array B has the same index range A, that is, B [p : r]. (One nice 
thing about pseudo-code, is that we can make these assumptions, and leave them up to the programmer to 
figure out how to implement it.) We have to indices i and j, that point to the current elements of each 
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sub-array. We move the smaller element into the next position of B (indicated by index k) and then 
increment the corresponding index (either i or j). When we run out of elements in one array, then we just 
copy the rest of the other array into B. Finally, we copy the entire contents of B back into A. (The use of 
the temporary array is a bit unpleasant, but this is impossible to overcome entirely. It is one of the 
shortcomings of MergeSort, compared to some of the other efficient sorting algorithms.) 

Here is the merge algorithm: 

merge( array A, int p, int q int r) 

1 int B[p..r]; int i <— k <— p; int j <— q + 1 

2 while (i < q) and (j < r) 

3 do if ( A[i] < A[j] ) 

4 then B [k++ ] <— A [i++ ] 

5 else B[k++ ] <— A[jf+ ] 

6 while (i < q ) 

7 do B [k++ ] <— A[b-+ ] 

8 while (j < r) 

9 do B [k++ ] <— A[jf+ ] 

10 for i <— p to r 

11 do A[i] B[i] 



Figure 3.2: Merge sort: combine phase 


3.1.1 Analysis of Merge Sort 

First let us consider the running time of the procedure Merge(A,p, q, r). Let n = r — p + 1 denote the 
total length of both the left and right sub-arrays. What is the running time of Merge as a function of n? 
The algorithm contains four loops (none nested in the other). It is easy to see that each loop can be 
executed at most n times. (If you are a bit more careful you can actually see that all the while-loops 
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together can only be executed n times in total, because each execution copies one new element to the 
array B, and B only has space for n elements.) Thus the running time to Merge n items is 0(n). Let us 
write this without the asymptotic notation, simply as n. (We’ll see later why we do this.) 

Now, how do we describe the running time of the entire MergeSort algorithm? We will do this through 
the use of a recurrence, that is, a function that is defined recursively in terms of itself. To avoid 
circularity, the recurrence for a given value of n is defined in terms of values that are strictly smaller than 
n. Finally, a recurrence has some basis values (e.g. for n = 1), which are defined explicitly. 

Let T(n) denote the worst case running time of MergeSort on an array of length n. If we call MergeSort 
with an array containing a single item (n = 1 ) then the running time is constant. We can just write 
T (n) = 1 , ignoring all constants. For n > 1 , MergeSort splits into two halves, sorts the two and then 
merges them together. The left half is of size pri/2] and the right half is |_n/2J . How long does it take to 
sort elements in sub array of size pri/2] ? We do not know this but because pn/2] < n for n > 1 , we can 
express this as T( |~n/2] ). Similarly the time taken to sort right sub array is expressed as T( |_n/2J ). In 
conclusion we have 

T(n) = / 1 if n = 1 , 

|T(|’n/2~|) +T(|_n/2J) +n otherwise 

This is called recurrence relation , i.e., a recursively defined function. Divide-and-conquer is an 
important design technique, and it naturally gives rise to recursive algorithms. It is thus important to 
develop mathematical techniques for solving recurrences, either exactly or asymptotically. 

Let’s expand the terms. 

T(1 ) = 1 

T(2) =T(1)+T(1)+2 = 1 +1 +2 = 4 
T(3) = T(2) + T(1 ) + 3 = 4 + 1 + 3 = 8 

T (4) = T(2) + T (2) + 4 =4 +4 +4 =1 2 

T(5) = T(3)+T(2) + 5 = 8 + 4 + 5 = 17 

T(8) = T(4) + T(4) + 8 = 12 + 12 + 8 = 32 

T(16) = T(8) + T(8) + 16 = 32 + 32 + 16 = 80 

T (32) =T(16) +T(16) +32 = 80 + 80 + 32 = 192 

What is the pattern here? Let’s consider the ratios T(n)/n for powers of 2: 

T(1 )/1 = 1 T(8)/8 = 4 

T(2)/2 = 2 T(16)/16 = 5 

T(4)/4 = 3 T(32)/32 = 6 


This suggests T(n)/n = logn + 1 Or, T(n) = nlogn + n which is 0(nlogn) (using the limit rule). 
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3.1.2 The Iteration Method for Solving Recurrence Relations 


Floor and ceilings are a pain to deal with. If n is assumed to be a power of 2 (2 k = n), this will simplify 
the recurrence to 


T(n) = 


2T(n/2) + n 


if n = 1 , 
otherwise 


The iteration method turns the recurrence into a summation. Let’s see how it works. Let’s expand the 
recurrence: 


T(n) = 2T(n/2) + n 

= 2(2T(n/4) + n/2) + n 
= 4T(n/4) + n + n 
= 4(2T(n/8) + n/4) + n + n 
= 8T(n/8) + n + n + n 
= 8(2T(n/16) + rt/8) + n + n + n 
= 16T(rt/16) + n + n + n + n 

If n is a power of 2 then let n = 2 k or k = log n. 

T(n) = 2 k T(n/(2 k )) + (n + n + nH hn) 

v y 

V* 

k times 

= 2 k T(n/(2 k )) + kn 
= 2 (logn) T(n/(2 ,logn) )) + (logn)n 
= 2 (logTV, T(n/n) + (logn)n 
= nT ( 1 ) + n log n = n + n log n 

3.1.3 Visualizing Recurrences Using the Recursion Tree 

Iteration is a very powerful technique for solving recurrences. But, it is easy to get lost in all the symbolic 
manipulations and lose sight of what is going on. Here is a nice way to visualize what is going on in 
iteration. We can describe any recurrence in terms of a tree, where each expansion of the recurrence takes 
us one level deeper in the tree. 

Recall that the recurrence for MergeSort (which we simplified by assuming that n is a power of 2, and 
hence could drop the floors and ceilings) 


T(n)= 1 ifn=1 ' 

I 2T(n/2) + n otherwise 

Suppose that we draw the recursion tree for MergeSort, but each time we merge two lists, we label that 
node of the tree with the time it takes to perform the associated (nonrecursive) merge. Recall that to 
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merge two lists of size m/2 to a list of size m takes 0(m) time, which we will just write as m. Below is 
an illustration of the resulting recursion tree. 


n 

time to merge 

= n 

+ 


n/2 n/2 

] 2(n/2) = n 


a! // 

+ 

^ 

log(n)+1 

n/4 n/4 n/4 

n/4 4(n/4) = n 

levels 

TV ZV A ZV 

n/8 n/8 n/8 n/8 n/8 n/8 n/8 

J n/8 8(n/8) = n 


i'i i'i i'i H i'M M 

1 r i + 


□ □ED] □□□□□ 

[T] n(n/n) = n 



n(log(n)+1) 



Figure 3.3: Merge sort Recurrence Tree 


3.1.4 A Messier Example 

The merge sort recurrence was too easy. Let’s try a messier recurrence. 

T(n) = {’ ifU=1 ' 

3T(rt/4) + n otherwise 


Assume n to be a power of 4, i.e., n — 4 k and k = log 4 n 

T [n] = 3T(n/4) + n 

= 3(3T (n/1 6) + (n/4) + n 
= 9T(n/16)+3(n/4)+n 
= 27T(n/64) + 9(n/l 6) + 3(n/4) + n 


= 3 k T(J)+3 k -’(n/4 k -') 

H F 9(n/16) + 3 (tl/ 4) + n 



+ L 
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With n — 4 k and T (1 ) = 1 


k-i 3i 

n , 3 


T (n)=3 t T(-)+X ? n 

i=0 


(log 4 n)-1 t 


3 ,og *"T(l) + Y. ^T n 
1=0 

(log 4 rt) — 1 t 

= n log 4 3 -|- y n 

Z_ Ax 


i=0 


We used the formula a logb n = n ,ogb a . rt remains constant throughout the sum and 3 L /4 L = (3/4 ) 1 
thus have 


T (n) = rt log4 3 + rt 

The sum is a geometric series; recall that for x ^ 1 


i=0 




i=0 


x m+1 _ -| 
X - 1 


In this case x = 3/4 and m = log 4 n — 1 . We get 

T(n) = n log4 3 + 


n- 


(3/4) log 4 n+1 - 1 
(3/4) - 1 

Applying the log identity once more 

(3 /4) lo g4 n — n log 4 (3/4) _ n log 4 3-log 4 4 


n 


log 4 3 


= n 


n 


If we plug this back, we get 


n log 4 3 


T(n) = n log4 3 + n-— v 


- 1 


(3/4) - 1 
n log 4 3 — rt 


— -p lo §4 3 _i_ 

-1/4 

= n log4 3 + 4(n — rt log4 3 ) 

= 4 n - 3rt log4 3 

With log 4 3 ~ 0.79, we finally have the result! 

T(n) =4 n — 3n log4 3 « 4 n — 3 n 0 ' 79 e 0(n) 


; we 
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3.2 Selection Problem 

Suppose we are given a set of n numbers. Define the rank of an element to be one plus the number of 
elements that are smaller. Thus, the rank of an element is its final position if the set is sorted. The 
minimum is of rank 1 and the maximum is of rank n. 

Consider the set: {5, 7, 2, 10, 8, 15, 21, 37, 41}. The rank of each number is its position in the sorted 
order. 


position 

1 

2 

3 

4 

5 

6 

7 

8 

9 

Number 

2 

5 

7 

8 

10 

15 

21 

37 

41 


For example, rank of 8 is 4, one plus the number of elements smaller than 8 which is 3. 

The selection problem is stated as follows: 

Given a set A of n distinct numbers and an integer k, 1 < k < n, output the element of A of rank 
k. 

Of particular interest in statistics is the median. If n is odd then the median is defined to be element of 
rank (n + 1 )/2. When n is even, there are two choices: n/2 and (n + 1 )/2. In statistics, it is common to 
return the average of the two elements. 

Medians are useful as a measures of the central tendency of a set especially when the distribution of 
values is highly skewed. For example, the median income in a community is a more meaningful measure 
than average. Suppose 7 households have monthly incomes 5000, 7000, 2000, 10000, 8000, 15000 and 
16000. In sorted order, the incomes are 2000, 5000, 7000, 8000, 10000, 15000, 16000. The median 
income is 8000; median is element with rank 4: (7 + 1 )/2 = 4. The average income is 9000. Suppose the 
income 16000 goes up to 450,000. The median is still 8000 but the average goes up to 71,000. Clearly, 
the average is not a good representative of the majority income levels. 

The selection problem can be easily solved by simply sorting the numbers of A and returning A[k]. 
Sorting, however, requires ©(ulog n) time. The question is: can we do better than that? In particular, is 
it possible to solve the selections problem in 0(n) time? The answer is yes. However, the solution is far 
from obvious. 

3.2.1 Sieve Technique 

The reason for introducing this algorithm is that it illustrates a very important special case of 
divide-and-conquer, which I call the sieve technique. We think of divide-and-conquer as breaking the 
problem into a small number of smaller subproblems, which are then solved recursively. The sieve 
technique is a special case, where the number of subproblems is just 1. 

The sieve technique works in phases as follows. It applies to problems where we are interested in finding 
a single item from a larger set of n items. We do not know which item is of interest, however after doing 
some amount of analysis of the data, taking say 0(nk) time, for some constant k, we find that we do not 
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know what the desired the item is, but we can identify a large enough number of elements that cannot be 
the desired value, and can be eliminated from further consideration. In particular “large enough” means 
that the number of items is at least some fixed constant fraction of n (e.g. n/2, n/3). Then we solve the 
problem recursively on whatever items remain. Each of the resulting recursive solutions then do the same 
thing, eliminating a constant fraction of the remaining set. 

3.2.2 Applying the Sieve to Selection 

To see more concretely how the sieve technique works, let us apply it to the selection problem. We will 
begin with the given array A[1 „n]. We will pick an item from the array, called the pivot element which 
we will denote by x. We will talk about how an item is chosen as the pivot later; for now just think of it as 
a random element of A. 

We then partition A into three parts. 

1. A[q] contains the pivot element x, 

2. A[1 ..q — 1] will contain all the elements that are less than x and 

3. A[q + 1 ..n] will contains all elements that are greater than x. 

Within each sub array, the items may appear in any order. The following figure shows a partitioned array: 



pivot 

p 

1 r 

5 9 2 

6 4 13 7 

P 

q r 

3 1 2 

4 6 9 5 7 

<x 

X > X 


Before partitioning 
After partitioning 


Figure 3.4: A[p..r] partitioned about the pivot x 


3.2.3 Selection Algorithm 

It is easy to see that the rank of the The rank of the pivot x is q — p + 1 in A[p..r] . Let 
rank.x = q— p + l.Ifk = rank.x then the pivot is k th smallest. If k < rank.x then search 
A[p..q — 1] recursively. If k > rank_x then search A[q + 1 ,.r\ recursively. Find element of rank (k — q) 
because we eliminated q smaller elements in A. 

SELECT( array A, int p, int r, int k) 
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1 if(p=r) 

2 then return A [p] 

3 else x <— choose_pivot(A,p,t) 

4 q <— partition(A,p,t,x) 

5 rarik.x <— q — p + 1 

6 if k = rarik.x 

7 then return x 

8 if k < rank_x 

9 then return select ( A, p, q — 1 , k) 

10 else return select( A, q + l,r,k— q) 

Example: select the 6 th smallest element of the set {5, 9, 2, 6, 4, 1, 3, 7} 


rank x = 4 



initial partition recurse partition recurse partition 

k=6 pivot=4 k=(6-4)=2 pivot=7 k=2 pivot=6 


DONE!! 


Figure 3.5: Sieve example: select 6 th smallest element 


3.2.4 Analysis of Selection 

We will discuss how to choose a pivot and the partitioning later. For the moment, we will assume that 
they both take 0(n) time. How many elements do we eliminate in each time? If x is the largest or the 
smallest then we may only succeed in eliminating one element. 

5, 9, 2, 6, 4, |T], 3,7 pivot is 1 
[T]> 5, 9, 2, 6, 4, 3, 7 after partition 

Ideally, x should have a rank that is neither too large or too small. 

Suppose we are able to choose a pivot that causes exactly half of the array to be eliminated in each phase. 
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This means that we recurse on the remaining n/2 elements. This leads to the following recurrence: 

T 


T(n) — 

If we expand this recurrence, we get 


if n = 1 , 
^T(n/2)+n otherwise 


, n n 

T(n) — rt+ — + — + ... 


< y — 

- Z— lx 

i =0 

OO 

=n>~i 

Z— 2 X 

i=0 


Recall the formula for infinite geometric series; for any |c| < 1 , 

CO 1 

Lc‘ = T^ 


i=0 


— C 


Using this we have 


T(n) < 2n e 0(n) 


Let’s think about how we ended up with a 0(n) algorithm for selection. Normally, a 0(n) time 
algorithm would make a single or perhaps a constant number of passes of the data set. In this algorithm, 
we make a number of passes. In fact it could be as many as log rt. 

However, because we eliminate a constant fraction of the array with each phase, we get the convergent 
geometric series in the analysis. This shows that the total running time is indeed linear in rt. This lesson 
is well worth remembering. It is often possible to achieve linear running times in ways that you would 
not expect. 



38 


CHAPTER 3. DIVIDE AND CONQUER STRATEGY 



Chapter 4 
Sorting 


For the next series of lectures, we will focus on sorting. There a number of reasons for sorting. Here are a 
few important ones. Procedures for sorting are parts of many large software systems. Design of efficient 
sorting algorithms is necessary to achieve overall efficiency of these systems. 

Sorting is well studied problem from the analysis point of view. Sorting is one of the few problems where 
provable lower bounds exist on how fast we can sort. In sorting, we are given an array A[1 ..n] of n 
numbers We are to reorder these elements into increasing (or decreasing) order. 

More generally, A is an array of objects and we sort them based on one of the attributes - the key value. 
The key value need not be a number. It can be any object from a totally ordered domain. Totally ordered 
domain means that for any two elements of the domain, x and y , either x < y , x = y or x > y . 


4.1 Slow Sorting Algorithms 

There are a number of well-known slow 0(n 2 ) sorting algorithms. These include the following: 


Bubble sort: Scan the array. Whenever two consecutive items are found that are out of order, swap 
them. Repeat until all consecutive items are in order. 


Insertion sort: Assume that A[1 ..i — 1] have already been sorted. Insert A[i] into its proper position in 
this sub array. Create this position by shifting all larger elements to the right. 

Selection sort: Assume that A[1 ..i — 1] contain the i — 1 smallest elements in sorted order. Find the 
smallest element in A[i..n] Swap it with A[i]. 


These algorithms are easy to implement. But they run in 0(rt 2 ) time in the worst case. 
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4.2 Sorting in O(nlogn) time 

We have already seen that Mergesort sorts an array of numbers in ©(nlogu) time. We will study two 
others: Heapsort and Quicksort. 


4.2.1 Heaps 

A heap is a left-complete binary tree that conforms to the heap order. The heap order property: in a 
(min) heap, for every node X, the key in the parent is smaller than or equal to the key in X. In other 
words, the parent node has key smaller than or equal to both of its children nodes. Similarly, in a max 
heap, the parent has a key larger than or equal both of its children Thus the smallest key is in the root in a 
min heap; in the max heap, the largest is in the root. 



Figure 4.1: A min-heap 


The number of nodes in a complete binary tree of height h. is 


H 

n = 2° + 2 1 + 2 2 + ■ ■ ■ + 2 h = Y_ 2i = 2 H+1 - 1 

i=0 


h. in terms of n is 


h. = (log(n + 1 )) — 1 ~ logn e 0(logn) 


One of the clever aspects of heaps is that they can be stored in arrays without using any pointers. This is 
due to the left-complete nature of the binary tree. We store the tree nodes in level-order traversal. Access 
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to nodes involves simple arithmetic operations: 


left(i) : returns 2i, index of left child of node i. 
righ.t(i) : returns 2i + 1 , the right child. 
parent(i) : returns |_i/2J , the parent of i. 


The root is at position 1 of the array. 


4.2.2 Heapsort Algorithm 


We build a max heap out of the given array of numbers A[1 „n]. We repeatedly extract the the maximum 
item from the heap. Once the max item is removed, we are left with a hole at the root. To fix this, we will 
replace it with the last leaf in tree. But now the heap order will very likely be destroyed. We will apply a 
heapify procedure to the root to restore the heap. Figure 4.2 shows an array being sorted. 


HEAPSORT( array A, int n) 

1 build-heap(A,tl) 

2 m <— n 

3 while (m > 2) 

4 do SWAP(A[1], A[m]) 

5 m <— m — 1 

6 heapify(A, 1 , m) 
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Figure 4.2: Example of heap sort 



















4.2. SORTING IN 0(N LOG N ) TIME 


43 


4.2.3 Heapify Procedure 

There is one principal operation for maintaining the heap property. It is called Heapify. (In other books it 
is sometimes called sifting down.) The idea is that we are given an element of the heap which we suspect 
may not be in valid heap order, but we assume that all of other the elements in the subtree rooted at this 
element are in heap order. In particular this root element may be too small. To fix this we “sift” it down 
the tree by swapping it with one of its children. Which child? We should take the larger of the two 
children to satisfy the heap ordering property. This continues recursively until the element is either larger 
than both its children or until its falls all the way to the leaf level. Here is the algorithm. It is given the 
heap in the array A, and the index i of the suspected element, and m the current active size of the heap. 
The element A [max] is set to the maximum of A[i] and it two children. If max ^ i then we swap A[t] 
and A [max] and then recurse on A [max]. 


HEAPIFY ( array A, int i, int m) 

1 l <— LEFT(i) 

2 r <— RIGHT(i) 

3 max <— i 

4 if (l < m)arid(A[l] > A[max]) 

5 then max <— l 

6 if (r < m)and(A[r] > A [max]) 

7 then max <— r 

8 if (max ^ i) 

9 then s wap ( A [i] , A [max] ) 

10 heapify(A, max, m) 


4.2.4 Analysis of Heapify 

We call heapify on the root of the tree. The maximum levels an element could move up is 0(log rt) levels. 
At each level, we do simple comparison which 0(1 ). The total time for heapify is thus O(logri). Notice 
that it is not 0(log n) since, for example, if we call heapify on a leaf, it will terminate in 0(1 ) time. 


4.2.5 BuildHeap 

We can use Heapify to build a heap as follows. First we start with a heap in which the elements are not in 
heap order. They are just in the same order that they were given to us in the array A. We build the heap by 
starting at the leaf level and then invoke Heapify on each node. (Note: We cannot start at the top of the 
tree. Why not? Because the precondition which Heapify assumes is that the entire tree rooted at node i is 
already in heap order, except for i.) Actually, we can be a bit more efficient. Since we know that each 
leaf is already in heap order, we may as well skip the leaves and start with the first non-leaf node. This 
will be in position n/2. (Can you see why?) 
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Here is the code. Since we will work with the entire array, the parameter m for Heapify, which indicates 
the current heap size will be equal to n, the size of array A, in all the calls. 


buildheap( array A, int n) 

1 for i <— rt/2 downto 1 

2 do 

3 heapify( A, i, rt) 


4.2.6 Analysis of BuildHeap 

For convenience, we will assume n = 2 H+1 — 1 where h is the height of tree. The heap is a left-complete 
binary tree. Thus at each level j , j < h, there are V nodes in the tree. At level H, there will be 2 h or less 
nodes. How much work does buildHeap carry out? Consider the heap in Figure 4.3: 



Figure 4.3: Total work performed in buildheap 


At the bottom most level, there are 2 h nodes but we do not heapify these. At the next level up, there are 
2 h 1 nodes and each might shift down 1. In general, at level j, there are 2 h ~i nodes and each may shift 
down j levels. 

So, if count from bottom to top, level-by-level, the total time is 

h h ~h 

T(n)=2j2 h -i=2j- 

j=0 j=0 

T(n)=2^i 

5=0 


We can factor out the 2 h term: 
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— x 


How do we solve this sum? Recall the geometric series, for any constant x < 1 

00 1 

£*' = t 

j=0 

Take the derivative with respect to x and multiply by x 

1 


i=o 


(1 -x) 


£jxi 


j=0 


(1 -x) 


We plug x = 1 /2 and we have the desired formula: 

1/2 


j=o 


1/2 

21 (1 — (1/2)) 2 1/4 




= 2 


In our case, we have a bounded sum, but since we the infinite series is bounded, we can use it as an easy 
approximation: 


H 

T(n)=2 h ^ 

1=0 


i_ 

21 


CO 



< 2 h • 2 = 2 h+1 


Recall that n = 2 H+1 — 1 . Therefore 


T(u) < n + 1 e O(n) 


The algorithm takes at least 12 (n) time since it must access every element at once. So the total time for 
BuildHeap is 0(n). 

BuildHeap is a relatively complex algorithm. Yet, the analysis yield that it takes 0(n) time. An intuitive 
way to describe why it is so is to observe an important fact about binary trees The fact is that the vast 
majority of the nodes are at the lowest level of the tree. For example, in a complete binary tree of height 
h, there is a total of n rs 2 h+1 nodes. 

The number of nodes at the bottom three levels alone is 

2 h + 2 H_1 + 2 h ~ 2 = t + T + ? = at = 0.875ri 
Z 4 o o 

Almost 90% of the nodes of a complete binary tree reside in the 3 lowest levels. Thus, algorithms that 
operate on trees should be efficient (as BuildHeap is) on the bottom-most levels since that is where most 
of the weight of the tree resides. 
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4.2.7 Analysis of Heapsort 

Heapsort calls BuildHeap once. This takes 0(n). Heapsort then extracts roughly n maximum elements 
from the heap. Each extract requires a constant amount of work (swap) and 0(log n) heapify. Heapsort 
is thus O(nlogn). 

Is HeapSort 0(nlogn)? The answer is yes. In fact, later we will show that comparison based sorting 
algorithms can not run faster than O(nlogn). Heapsort is such an algorithm and so is Mergesort that we 
saw ealier. 


4.3 Quicksort 


Our next sorting algorithm is Quicksort. It is one of the fastest sorting algorithms known and is the 
method of choice in most sorting libraries. Quicksort is based on the divide and conquer strategy. Here is 
the algorithm: 


QUICKSORT) array A, int p, int r) 

1 if (r > p) 

2 then 

3 i <— a random index from [p..r] 

4 swap A[i] with A[p] 

5 q <— partition(A,p,t) 

6 quicksort) A, p, q - 1) 

7 quicksort(A, q + 1 , r) 


4.3.1 Partition Algorithm 

Recall that the partition algorithm partitions the array A[p..r] into three sub arrays about a pivot element 
x. A[p..q — 1] whose elements are less than or equal to x, A[q] = x, A[q + 1 ..r] whose elements are 
greater than x 

We will choose the first element of the array as the pivot, i.e. x = A[p]. If a different rule is used for 
selecting the pivot, we can swap the chosen element with the first element. We will choose the pivot 
randomly. 

The algorithm works by maintaining the following invariant condition. A[p] = x is the pivot value. 
A[p..q — 1] contains elements that are less than x, A[q + 1 ,.s — 1] contains elements that are greater than 
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or equal to x A[s..r] contains elements whose values are currently unknown. 


PARTITION ( array A, int p, int r) 

1 x < — A[p] 

2 q <— p 

3 for s <— p + 1 to r 

4 do if (A[s] < x) 

5 then q <— q + 1 

6 swap A[q] with A[s] 

7 swap A[p] with A[q] 

8 return q 


Figure 4.4 shows the execution trace partition algorithm. 



Figure 4.4: Trace of partitioning algorithm 


4.3.2 Quick Sort Example 

The Figure 4.5 trace out the quick sort algorithm. The first partition is done using the last element, 1 0, of 
the array. The left portion are then partitioned about 5 while the right portion is partitioned about 1 3. 
Notice that 1 0 is now at its final position in the eventual sorted order. The process repeats as the 
algorithm recursively partitions the array eventually sorting it. 
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Figure 4.5: Example of quick sort 
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It is interesting to note (but not surprising) that the pivots form a binary search tree. This is illustrated in 
Figure 4.6. 



Figure 4.6: Quicksort BST 


4.3.3 Analysis of Quicksort 

The running time of quicksort depends heavily on the selection of the pivot. If the rank of the pivot is 
very large or very small then the partition (BST) will be unbalanced. Since the pivot is chosen randomly 
in our algorithm, the expected running time is O(nlogn). The worst case time, however, is 0(n 2 ). 
Luckily, this happens rarely. 


4.3.4 Worst Case Analysis of Quick Sort 

Let’s begin by considering the worst-case performance, because it is easier than the average case. Since 
this is a recursive program, it is natural to use a recurrence to describe its running time. But unlike 
MergeSort, where we had control over the sizes of the recursive calls, here we do not. It depends on how 
the pivot is chosen. Suppose that we are sorting an array of size n, A[1 : tl], and further suppose that the 
pivot that we select is of rank q, for some q in the range 1 to n. It takes 0(n) time to do the partitioning 
and other overhead, and we make two recursive calls. The first is to the subarray A[1 : q — 1] which has 
q — 1 elements, and the other is to the subarray A[q + 1 : n] which has n — q elements. So if we ignore 
the 0(n) (as usual) we get the recurrence: 


T(n) = T(q — 1 ) + T (n — q ) + rt 
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This depends on the value of q. To get the worst case, we maximize over all possible values of q. Putting 
is together, we get the recurrence 

f 1 if n < 1 

~ 1 max (T(q — 1 ) + T(n — q) + n) otherwise 

Recurrences that have max’s and min’s embedded in them are very messy to solve. The key is 
determining which value of q gives the maximum. (A rule of thumb of algorithm analysis is that the 
worst cases tends to happen either at the extremes or in the middle. So I would plug in the value q = 1 , 
q = Tt, and q = n/2 and work each out.) In this case, the worst case happens at either of the extremes. If 
we expand the recurrence for q = 1 , we get: 

T(n) < T(0) +T(n— 1) +n 
= 1 + T(n — 1 ) + n 
= T(n — 1) + (n+ 1) 

= T(n - 2) + n + (n + 1 ) 

= T(n — 3) + (n- 1) +n + (n + 1) 

= T(n — 4) + (n - 2) + (n — 1 ) + n + (n + 1 ) 

k—2 

= T(n — k) + ^ (n — i) 

i=— 1 


For the basis T ( 1 


= 1 we set k = n — 1 and get 

n3 

T(n) <T(1) + Jjn-i 


i =— 1 


— 1 + (3 + 4 + 5 + -- - + (n — 1 ) + n + (n + 1 )) 


n+t 




< > t = 

i=1 


(n+ l)(n + 2) 


G 0(n 2 


4.3.5 Average-case Analysis of Quicksort 

We will now show that in the average case, quicksort runs in 0(rilog n) time. Recall that when we talked 
about average case at the beginning of the semester, we said that it depends on some assumption about 
the distribution of inputs. However, in the case of quicksort, the analysis does not depend on the 
distribution of input at all. It only depends upon the random choices of pivots that the algorithm makes. 
This is good, because it means that the analysis of the algorithm’s performance is the same for all inputs. 
In this case the average is computed over all possible random choices that the algorithm might make for 
the choice of the pivot index in the second step of the Quicksort procedure above. 

To analyze the average running time, we let T (n) denote the average running time of Quicksort on a list 
of size n. It will simplify the analysis to assume that all of the elements are distinct. The algorithm has n 
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random choices for the pivot element, and each choice has an equal probability of 1 /n of occuring. So 
we can modify the above recurrence to compute an average rather than a max, giving: 


T(n) = 


1 if n < 1 

n^q=i(T(q - 1) +T(n- q) +n) otherwise 


The time T (n) is the weighted sum of the times taken for various choices of q. I.e., 


T(n) 


— (T(0) + T(n — 1 ) + n) 

n 

+ — (T(1 ) + T(n — 2) + n) 

n 

+ -(T(2)+T(n-3)+n) 

n 

+ •■■ 4 — (T(n — 1) + T(0) + n) 

TL 


We have not seen such a recurrence before. To solve it, expansion is possible but it is rather tricky. We 
will attempt a constructive induction to solve it. We know that we want a 0(nlogn). Let us assume that 
T (n) < cnlog n for n > 2 where c is a constant. 

For the base case n = 2we have 


T(n) = i£(T(c-l) + T( 2 -c,) + 2 ) 

Z q=1 

= 1 [fT(0J + T(T) + 2) + (T(l) + T(0) +2) 


We want this to be at most c2 log 2, i.e., 

T(2) < c2 log 2 


or 

4 < c2 log 2 


therefore 

c >4/(2 log 2) « 2.88. 

For the induction step, we assume that n > 3 and The induction hypothesis is that for any n' < n, we 
have T (n') < cn' log n'. We want to prove that it is true for T (n). By expanding T (rt) and moving the 
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factor of n outside the sum we have 


T(n) = - Y (T(q - 1 ) + T(n - q) + n) 
q=l 

n 

Y (T(q — 1)+T(n — q))+n 


1 

n 


1 

n 


q=1 

n 


^T(q-1 

q=l 


n 


Y_ T(n - q) + n 


q=l 


T(n) = - Y T(q - 1) + - Y T(n- q) +n 
n * — n — 


q=l 


q=l 


Observe that the two sums add up the same values T(0) + T(1 ) + • • • + T(n — 1 ). One counts up and 
other counts down. Thus we can replace them with 2 X.q=o T(q). We will extract T(0) and T(1 ) and treat 
them specially. These two do not follow the formula. 


n-1 


T (n) = -(^ T (q)J + n 


q=0 


= ^(T(0)+T(l) + ^T(q)) +n 

' q=2 


n-1 


We will apply the induction hypothesis for q < n we have 


n-1 


T(rt) = — ^T(O) + T(1) + Y_ T(q)j + n 


q=2 


n-1 


< ~(l + 1 + Y (cqlogq)j + n 

q=2 

2° i A 4 

= — ( _>_(cqln q) ) +n+ - 
n \ z — / n 

q=2 

We have never seen this sum before: 

n-1 

S(n) = ^(cqlnq) 

q=2 

Recall from calculus that for any monotonically increasing function f (x) 


b-l j-b 

Y f (i) < f(x)dx 
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The function f(x) = xlnx is monotonically increasing, and so 


We can integrate this by parts: 


n— 1 

S(n) = ^(cqln q) < 

q=2 


'IT 

xln x dx 

2 


'TT 

xln x dx 

2 



X 2 n 

4 x=2 


(4.1) 


xln x dx = —In x — — 


n 

x—2 

2 . 


We thus have 


= (ylnn- y ) - (21n 2 — 1 

2 2 

nr- nr 

-Y lnn -T 


n-1 y 2 

S(n) = ^(cqlnq) < -Inn- — 

q=2 


Plug this back into the expression for T (n) to get 


2c / n 2 tv* 

T( n ) = — (—Inn - — ) + n + 
n V 2 4 


4 

n 


. 2c /n 2 , n 2 \ 

T (n) = _(_ in n -_) +n+ 

, cn 4 

= cnln n — + n H — 

2 n 

c 4 

= cnlnn + n(l — 4 — 

v 2 n 


4 

n 


c 4 

T(n) = cnlnn + n(l — 4 — 

v V n 

To finish the proof, we want all of this to be at most cnln n. For this to happen, we will need to select c 
such that 


n(l — 

If we select c = 3, and use the fact that n > 3 we 

n (l ~j) + 


< 



< 0 


get 


4 3 n 


n 

1 


n 2 

3 1 

2~~2 


< 0. 


From the basis case we had c > 2.88. Choosing c = 3 satisfies all the constraints. Thus 
T(n) = 3nlnn e 0(nlogn). 
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4.4 In-place, Stable Sorting 

An in-place sorting algorithm is one that uses no additional array for storage. A sorting algorithm is 
stable if duplicate elements remain in the same relative position after sorting. 


9 

3 

3' 

5 

6 

5' 

2 

1 

3" 

1 

2 

3 

3' 

3" 

5 

5' 

6 

9 

1 

2 

3' 

3 

3" 

5' 

5 

6 

9 


unsorted 
stable sort 
unstable 


Bubble sort, insertion sort and selection sort are in-place sorting algorithms. Bubble sort and insertion 
sort can be implemented as stable algorithms but selection sort cannot (without significant modifications). 
Mergesort is a stable algorithm but not an in-place algorithm. It requires extra array storage. Quicksort is 
not stable but is an in-place algorithm. Heapsort is an in-place algorithm but is not stable. 


4.5 Lower Bounds for Sorting 

The best we have seen so far is O(nlogn) algorithms for sorting. Is it possible to do better than 
O(rilogn)? If a sorting algorithm is solely based on comparison of keys in the array then it is impossible 
to sort more efficiently than £1 ( n log n ) time. All algorithms we have seen so far are comparison-based 
sorting algorithms. 

Consider sorting three numbers ai , a 2 , < 13 . There are 3! = 6 possible combinations: 


(ai, a. 2 , a 3 ), (ai, a 3 , a 2 ) , (a 3 , a 2 , aQ 
(a 3 , ai, a 2 ), (a 2 , ai, a 3 ) , (a 2 , a 3 , ai) 


One of these permutations leads to the numbers in sorted order. 

The comparison based algorithm defines a decision tree. Here is the tree for the three numbers. 
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Figure 4.7: Decision Tree 


For n elements, there will be n! possible permutations. The height of the tree is exactly equal to T(rt), 
the running time of the algorithm. The height is T (tl) because any path from the root to a leaf 
corresponds to a sequence of comparisons made by the algorithm. 

Any binary tree of height T(n) has at most 2~^ nl leaves. Thus a comparison based sorting algorithm can 
distinguish between at most 2 T(n) different final outcomes. So we have 

2 T ' :m > n! and therefore 
T(n) > log(n!) 


We can use Stirling’s approximation for n!: 

n! > 


Thereofore 

T(rt) > log^v/ 27 m^j j 

= \og{V2Ttn) + nlogn — nlog e 
G fl(nlogn) 


We thus have the following theorem. 

Theorem 1 

Any comparison-based sorting algorithm has worst-case running time jQ(nlog n). 
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Chapter 5 

Linear Time Sorting 


The lower bound implies that if we hope to sort numbers faster than 0(nlog n), we cannot do it by 
making comparisons alone. Is it possible to sort without making comparisons? The answer is yes, but 
only under very restrictive circumstances. Many applications involve sorting small integers (e.g. sorting 
characters, exam scores, etc.). We present three algorithms based on the theme of speeding up sorting in 
special cases, by not making comparisons. 


5.1 Counting Sort 


We will consider three algorithms that are faster and work by not making comparisons. Counting sort 
assumes that the numbers to be sorted are in the range 1 to k where k is small. The basic idea is to 
determine the rank of each number in final sorted array. 

Recall that the rank of an item is the number of elements that are less than or equal to it. Once we know 
the ranks, we simply copy numbers to their final position in an output array. 

The question is how to find the rank of an element without comparing it to the other elements of the 
array?. The algorithm uses three arrays. As usual, A[1 ..n] holds the initial input, B[1 ..n] holds the sorted 
output and C [1 ..k] is an array of integers. C[x] is the rank of x in A, where x 6 [1 . .k] . The algorithm is 
remarkably simple, but deceptively clever. The algorithm operates by first constructing C. This is done in 
two steps. First we set C[x] to be the number of elements of A[j] that are equal to x. We can do this 
initializing C to zero, and then for each j, from 1 to n, we increment C[A[j]] by 1 . Thus, if A[j] = 5, then 
the 5th element of C is incremented, indicating that we have seen one more 5. To determine the number 
of elements that are less than or equal to x, we replace C[x] with the sum of elements in the sub array 
R[1 : x]. This is done by just keeping a running total of the elements of C. 

C[x] now contains the rank of x. This means that if x = A[j] then the final position of A[j] should be at 
position C[x] in the final sorted array. Thus, we set B[C[x]] = A[j]. Notice We need to be careful if there 
are duplicates, since we do not want them to overwrite the same location of B. To do this, we decrement 


57 



58 


CHAPTER 5. LINEAR TIME SORTING 


C[i] after copying. 


C0UNTING-S0RT( array A, array B, int k) 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 


1 


to 

0 


k times 


for i 

do C[i] <— 

for j <- 1 to length[A] 
do C[A[j]] <— C[A[j]] + 1 n times 
// C[i] now contains the number of elements = i 

for i <— 2 to k 

do C [i] 


C[i] + C[i-1] 


k times 


// C[i] now contains the number of elements < i 
for j <— length[A] downto 1 
do B[C[A[j]]] <— A[j] 

C[A[j]] <- C [A [j]] - 1 


n times 


There are four (unnested) loops, executed k times, n times, k — 1 times, and rt times, respectively, so the 
total running time is 0(n + k) time. If k = 0(n), then the total running time is 0(rt). 


Figure 5.1 through 5.19 shows an example of the algorithm. You should trace through the example to 
convince yourself how it works. 



Figure 5.1: Initial A and C arrays. 
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Figure 5.2: A[l] = 7 processed 



Figure 5.3: A [2] = 1 processed 



Figure 5.4: A [3] =3 processed 
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Figure 5.7: C now contains count of elements of A 
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Figure 5.8: C set to rank each number of A 


Input 

A[L.n] 

Output 

B[l..n\ 



Figure 5.9: A[1 1] placed in output array B 
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Figure 5.10: A[1 0] placed in output array B 


Input 

A[L.n\ 

Output 



Figure 5.11: A [9] placed in output array B 
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Input 

A[l..n\ 

Output 

B[l..n\ 



Figure 5.12: A [8] placed in output array B 



Figure 5.13: A [7] placed in output array B 
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Figure 5.14: A [6] placed in output array B 



Figure 5.15: A [5] placed in output array B 
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Figure 5.16: A [4] placed in output array B 


Input 

A[l..n] 

Output 

/?[!.. ft] I 
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Figure 5.17: A [3] placed in output array B 
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Figure 5.18: A [2] placed in output array B 



Figure 5.19: B now contains the final sorted data. 


Counting sort is not an in-place sorting algorithm but it is stable. Stability is important because data are 
often carried with the keys being sorted, radix sort (which uses counting sort as a subroutine) relies on it 
to work correctly. Stability achieved by running the loop down from n to 1 and not the other way around: 
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COUNTING-sort( array A, array B, int k) 
1 : 

2 for j <— length[A] downto 1 

3 do B[C[A[j]]] <— A[j] 

4 C [A [j]] <- C[A[j]] - 1 


Figure 5.20 illustrates the stability. The numbers 1, 2, 3, 4, and 7, each appear twice. The two 4’s have 
been given the superscript Numbers are placed in the output B array starting from the right. The two 
4’s maintain their relative position in the B array. If the sorting algorithm had caused 4** to end up on the 
left of 4*, the algorithm would be termed unstable. 




68 


CHAPTER 5. LINEAR TIME SORTING 





Figure 5.20: Stability of counting sort 
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5.2 Bucket or Bin Sort 

Assume that the keys of the items that we wish to sort lie in a small fixed range and that there is only one 
item with each value of the key. Then we can sort with the following procedure: 

1. Set up an array of “bins” - one for each value of the key - in order, 

2. Examine each item and use the value of the key to place it in the appropriate bin. 


Now our collection is sorted and it only took n operations, so this is an O(n) operation. However, note 
that it will only work under very restricted conditions. To understand these restrictions, let’s be a little 
more precise about the specification of the problem and assume that there are m values of the key. To 
recover our sorted collection, we need to examine each bin. This adds a third step to the algorithm above, 


3. Examine each bin to see whether there’s an item in it. 


which requires m operations. So the algorithm’s time becomes: 

T(rt) = Ci n + c 2 m 

and it is strictly 0(n + m). If m < n, this is clearly 0(n). However if m >> n, then it is O(m). An 
implementation of bin sort might look like: 


BucektSort( array A, int n, int M) 

1 // Pre-condition: for 1 < i < n, 0 < a[i] < M 

2 // Mark all the bins empty 

3 for i <— 1 to M 

4 do bin[i] <— Empty 

5 for i <— 1 to n 

6 do biu[A[i]] <— A[i] 


If there are duplicates , then each bin can be replaced by a linked list. The third step then becomes: 


3. Link all the lists into one list. 


We can add an item to a linked list in 0(1 ) time. There are n items requiring 0(n) time. Linking a list to 
another list simply involves making the tail of one list point to the other, so it is 0(1). Linking m such 
lists obviously takes 0(m) time, so the algorithm is still 0(n + m). Ligures 5.21 through 5.23 show the 
algorithm in action using linked lists. 
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Figure 5.21: Bucket sort: step 1, placing keys in bins in sorted order 



Figure 5.22: Bucket sort: step 2, concatenate the lists 
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Figure 5.23: Bucket sort: the final sorted sequence 


5.3 Radix Sort 


The main shortcoming of counting sort is that it is useful for small integers, i.e., 1 ..k where k is small. If 
k were a million or more, the size of the rank array would also be a million. Radix sort provides a nice 
work around this limitation by sorting numbers one digit at a time. 


576 49 [4] 9[5]4 [1]76 176 

494 19 [4] 5 [7] 6 [1]94 194 

194 95 [4] 1[7]6 [2]78 278 

296 =$■ 57[6] =>■ 2[7]8 =>■ [2]96 =$■ 296 

278 29[6] 4[9]4 [4]94 494 

176 17[6] 1[9]4 [5]76 576 

954 27[8] 2[9]6 [9]54 954 

Here is the algorithm that sorts A[1 ..it] where each number is d digits long. 


radix-sort( array A, int n, int d) 

1 for i <— 1 to d 

2 do stably sort A w.r.t i th lowest order digit 
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Chapter 6 

Dynamic Programming 

6.1 Fibonacci Sequence 


Suppose we put a pair of rabbits in a place surrounded on all sides by a wall. How many pairs of rabbits 
can be produced from that pair in a year if it is supposed that every month each pair begets a new pair 
which from the second month on becomes productive? Resulting sequence is 
1,1,2,3,5,8,13,21, 34, 55, . . . where each number is the sum of the two preceding numbers. 

This problem was posed by Leonardo Pisano, better known by his nickname Fibonacci (son of Bonacci, 
born 1 170, died 1250). This problem and many others were in posed in his book Liber abaci, published 
in 1202. The book was based on the arithmetic and algebra that Fibonacci had accumulated during his 
travels. The book, which went on to be widely copied and imitated, introduced the Hindu- Arabic 
place-valued decimal system and the use of Arabic numerals into Europe. The rabbits problem in the 
third section of Liber abaci led to the introduction of the Fibonacci numbers and the Fibonacci sequence 
for which Fibonacci is best remembered today. 

This sequence, in which each number is the sum of the two preceding numbers, has proved extremely 
fruitful and appears in many different areas of mathematics and science. The Fibonacci Quarterly is a 
modern journal devoted to studying mathematics related to this sequence. The Fibonacci numbers F n are 
defined as follows: 


F o = 0 

Ft = 1 

Fn = F n— 1 + F n— 2 


73 



74 


CHAPTER 6. DYNAMIC PROGRAMMING 


The recursive definition of Fibonacci numbers gives us a recursive algorithm for computing them: 


fib(tl) 

1 if (n < 2) 

2 then return n 

3 else return fib (n — 1 ) + fib (n — 2) 


Figure ?? shows four levels of recursion for the call fib (8): 



fib(4) fib(3) fib(3) fib(2) fib(3) fib(2) fib(2) fib(l) fib(3) fib(2) fib(2) fib(l) fib(2) fib(l) fib(l) fib(O) 
Figure 6.1: Recursive calls during computation of Fibonacci number 


A single recursive call to fib(n) results in one recursive call to fib(n — 1), two recursive calls to 
fib(n — 2), three recursive calls to fib(u — 3), five recursive calls to fib(n — 4) and, in general, F k _i 
recursive calls to fib(n — k) For each call, we’re recomputing the same fibonacci number from scratch. 

We can avoid this unnecessary repetitions by writing down the results of recursive calls and looking them 
up again if we need them later. This process is called memoization. Here is the algorithm with 
memoization. 


MemoFib(ti) 

1 if (n < 2) 

2 then return n 

3 if (F[n] is undefined) 

4 then F[n] <— MemoFib (n - 1 ) + MemoFib (n — 2) 

5 return F[n] 


If we trace through the recursive calls to MemoFib, we find that array FO gets filled from bottom up. I.e., 
first F[2], then F[3], and so on, up to F[n]. We can replace recursion with a simple for-loop that just fills 
up the array FO in that order. 
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This gives us our first explicit dynamic programming algorithm. 


IterFib(ti) 

1 F[0] i — 0 

2 F[1] <- 1 

3 for i <— 2 to n 

4 do 

5 F[i] <— F[i — 1] + F[i — 2] 

6 return F[n] 


This algorithm clearly takes only O(n) time to compute F n . By contrast, the original recursive algorithm 
takes 0(4) n ), <j) = 1 1 « 1 .61 8. IterFib achieves an exponential speedup over the original recursive 

algorithm. 


6.2 Dynamic Programming 

Dynamic programming is essentially recursion without repetition. Developing a dynamic programming 
algorithm generally involves two separate steps: 

• Formulate problem recursively. Write down a formula for the whole problem as a simple 
combination of answers to smaller subproblems. 

• Build solution to recurrence from bottom up. Write an algorithm that starts with base cases and 
works its way up to the final solution. 

Dynamic programming algorithms need to store the results of intermediate subproblems. This is often 
but not always done with some kind of table. We will now cover a number of examples of problems in 
which the solution is based on dynamic programming strategy. 


6.3 Edit Distance 

The words “computer” and “commuter” are very similar, and a change of just one letter, p-^m, will 
change the first word into the second. The word “sport” can be changed into “sort” by the deletion of the 
‘p’, or equivalently, ‘sort’ can be changed into ‘sport’ by the insertion of ‘p’. The edit distance of two 
strings, si and s2, is defined as the minimum number of point mutations required to change si into s2, 
where a point mutation is one of: 

• change a letter, 


• insert a letter or 
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• delete a letter 

For example, the edit distance between FOOD and MONEY is at most four: 

FOOD — > MOOD — > MON . D 

A 

— > MONED — > MONEY 


6.3.1 Edit Distance: Applications 

There are numerous applications of the Edit Distance algorithm. Here are some examples: 

Spelling Correction 

If a text contains a word that is not in the dictionary, a ‘close’ word, i.e. one with a small edit distance, 
may be suggested as a correction. Most word processing applications, such as Microsoft Word, have 
spelling checking and correction facility. When Word, for example, finds an incorrectly spelled word, it 
makes suggestions of possible replacements. 

Plagiarism Detection 

If someone copies, say, a C program and makes a few changes here and there, for example, change 
variable names, add a comment of two, the edit distance between the source and copy may be small. The 
edit distance provides an indication of similarity that might be too close in some situations. 

Computational Molecular Biology DNA is a polymer. The monomer units of DNA are nucleotides, and 
the polymer is known as a “polynucleotide.” Each nucleotide consists of a 5-carbon sugar (deoxyribose), 
a nitrogen containing base attached to the sugar, and a phosphate group. There are four different types of 
nucleotides found in DNA, differing only in the nitrogenous base. The four nucleotides are given one 
letter abbreviations as shorthand for the four bases. 

• A- adenine 

• G-guanine 

• C-cytosine 

• T-thymine 

Double-helix of DNA molecule with nucleotides Figure of Double-helix of DNA molecule with 
nucleotides goes here 

The edit distance like algorithms are used to compute a distance between DNA sequences (strings over 
A,C,G,T, or protein sequences (over an alphabet of 20 amino acids), for various purposes, e.g.: 

• to find genes or proteins that may have shared functions or properties 

• to infer family relationships and evolutionary trees over different organisms. 
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Speech Recognition 

Algorithms similar to those for the edit-distance problem are used in some speech recognition systems. 
Find a close match between a new utterance and one in a library of classified utterances. 

6.3.2 Edit Distance Algorithm 

A better way to display this editing process is to place the words above the other: 

S D 1 M D M 

M A I T H 

A _ R T _ S 

The first word has a gap for every insertion (I) and the second word has a gap for every deletion (D). 
Columns with two different characters correspond to substitutions (S). Matches (M) do not count. The 
Edit transcript is defined as a string over the alphabet M, S, I, D that describes a transformation of one 
string into another. For example 

S D I M D M 

1+ 1+ 1+ 0+ 1+ 0+ = 4 

In general, it is not easy to determine the optimal edit distance. For example, the distance between 
ALGORITHM and ALTRUISTIC is at most 6. 

ALGOR I THM 
AL TRUISTIC 


Is this optimal? 

6.3.3 Edit Distance: Dynamic Programming Algorithm 

Suppose we have an m-character string A and an ri-character string B. Define E(i, j) to be the edit 
distance between the first i characters of A and the first j characters of B. For example, 

1 

A L G O R I T H M 
A L _ TRUISTIC 

' v ^ 


The edit distance between entire strings A and B is E(m, n). The gap representation for the edit 
sequences has a crucial “ optimal substructure ’ property. If we remove the last column, the remaining 
columns must represent the shortest edit sequence for the remaining substrings. The edit distance is 6 for 
the following two words. 



78 


CHAPTER 6. DYNAMIC PROGRAMMING 


ALGOR I THM 
AL TRUISTIC 

If we remove the last column, the edit distance reduces to 5. 

ALGOR _I_TH 
A L T R U I S T I 

We can use the optimal substructure property to devise a recursive formulation of the edit distance 
problem. There are a couple of obvious base cases: 

• The only way to convert an empty string into a string of j characters is by doing j insertions. Thus 

E(0,j)=j 

• The only way to convert a string of i characters into the empty string is with i deletions: 

E(i, 0) = i 

There are four possibilities for the last column in the shortest possible edit sequence: 

Deletion: Last entry in bottom row is empty. 

i=3 

A L G O R I T H M 

A L TRUISTIC 

j=2 

In this case 

E(i,j)=E(i-l,j) + 1 

Insertion: The last entry in the top row is empty. 

i—5 

A L O R I T H M 

AL _ TRUISTIC 

V v y 

j=5 


In this case 

E(i, j) = E(i, j — 1 ) + 1 
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Substitution: Both rows have characters in the last column. 


i=4 

A L G O' R I T H M 
A L _ TRUISTIC 

v * v ^ 

i=3 

If the characters are different, then 

E(r,j)=E(i-1J-1) + l 

1=5 

'algor' I t h m 
AL _ TRUISTIC 

v v 

v 

j=4 


If characters are same, no substitution is needed: 


E(i, j) = E(i — I , j — 1 ) 


Thus the edit distance E(i, j) is the smallest of the four possibilities: 

( E(i — 1 , j) + 1 \ 

p n _ • E(i, j — 1) + 1 

1 MJ mm E(i — 1 , j — 1 ) + 1 ifA[i]^B[j] 

\ E(i-lJ-l) if A[i] = B[j]y 

Consider the example of edit between the words “ARTS” and “MATHS”: 


ARTS 

MATHS 


The edit distance would be in E(4, 5). If we recursion to compute, we will have 


E(4, 5) = min 


/E(3, 5) + 1 \ 

E (4, 4) + 1 

E(3,4) + 1 if A [4] ^ B [5] 

V E(3,4) if A [4] = B[5]/ 


Recursion clearly leads to the same repetitive call pattern that we saw in Fibonnaci sequence. To avoid 
this, we will use the DP approach. We will build the solution bottom-up. We will use the base case 
E(0, j) to fill first row and the base case E(i, 0) to fill first column. We will fill the remaining E matrix 
row by row. 
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A 

R 

T 

S 


0 

->1 

— >2 

— >3 

— >4 

M 

1 





A 

2 





T 

3 





H 

4 





S 

5 








A 

R 

T 

S 


0 

->1 

— >2 

— >3 

— >4 

M 






A 






T 






H 






S 







Table 6.1: First row and first column entries using the base cases 


We can now fill the second row. The table not only shows the values of the cells E[i, j] but also arrows 
that indicate how it was computed using values in E[i — 1 , j], E[i, j — 1] and E[i — 1 , j — 1]. Thus, if a cell 
E [1, j] has a down arrow from E[i — 1 , j] then the minimum was found using E[i — 1 , j]. For a right arrow, 
the minimum was found using E[i, j — 1]. For a diagonal down right arrow, the minimum was found 
using E [i — 1 , j — 1], There are certain cells that have two arrows pointed to it. In such a case, the 
minimum could be obtained from the diagonal E[i — 1 , j — 1] and either of E[i — 1 , j] and E[i, j — 1], We 
will use these arrows later to determine the edit script. 
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A 

R 

T 

S 


0 

^1 

— >2 

— >3 

— >4 

M 

1 

\ 

1 

\ 

— >2 



A 

2 





T 

3 





H 

4 





S 

5 








A 

R 

T 

S 


0 


— >2 

— >3 

— >4 

M 

1 

\ 

1 




A 

2 





T 

3 





H 

4 





S 

5 






Table 6.2: Computing E[1 , 1] and E[1 , 2] 




A 

R 

T 

S 


0 


— >2 

— >3 

— >4 

M 

1 

\ 

1 

\ 

— >2 

\ 

— >3 

\ 

— >4 

A 

2 





T 

3 





H 

4 





S 

5 








A 

R 

T 

S 


0 

->1 

— >2 

— >3 

— >4 

M 

1 

\ 

1 

\ 

— >2 

\ 

— >3 


A 

2 





T 

3 





H 

4 





S 

5 






Table 6.3: Computing E[1 , 3] and E[1 ,4] 


An edit script can be extracted by following a unique path from E[0, 0] to E [4, 5]. There are three possible 
paths in the current example. Let us follow these paths and compute the edit script. In an actual 
implementation of the dynamic programming version of the edit distance algorithm, the arrows would be 
recorded using an appropriate data structure. For example, each cell in the matrix could be a record with 
fields for the value (numeric) and flags for the three incoming arrows. 
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A 

R 

T 

S 


0 

->1 

— >2 

— >3 

— >4 

M 


\ 

\ 

\ 

\ 

1 

1 

— >2 

— >3 

— >4 

A 


\ 

\ 

\ 

\ 


2 

1 

— >2 

— >3 

— >4 

T 

3 

2 

\ 

2 

\ 

2 

— >3 

H 



V 

V 

\ 

4 

3 

3 

3 

3 

<5 



V 

V 

\ 


5 

4 

4 

4 

3 


Table 6.4: The final table with all T[i, j] entries computed 




A 

R 

T 

S 


0 

->1 

— >2 

— >3 

— >4 

M 

1 

\ 

1 

\ 

— >2 

\ 

— >3 

\ 

— >4 

A 

2 

\ 

1 

\ 

— >2 

\ 

^3 

\ 

— >4 

T 

3 

2 

\ 

2 

\ 

2 

— >3 

H 

4 

3 

V 

3 

V 

3 

\ 

3 

S 

5 

4 

V 

4 

V 

4 

\ 

3 


Solution path 1: 

1 + 0 + 1 + 1 + 
D M S S 


0 =3 

M 


MATH 
_ ART 


S 

s 


6.3. EDIT DISTANCE 


83 




A 

R 

T 

S 


0 

->1 

— >2 

^3 

— >4 

M 


\ 

\ 

\ 

\ 

1 

1 

— >2 

^3 

— >4 

A 


\ 

\ 

\ 

\ 


2 

1 

— >2 

^3 

— >4 

T 

3 

2 

\ 

2 

\ 

2 

— >3 

H 



V 

V 

\ 

4 

3 

3 

3 

3 

<5 



V 

V 

\ 


5 

4 

4 

4 

3 


Table 6.5: Possible edit scripts. The red arrows from E[0, 0] to E[4, 5] show the paths that can be followed 
to extract edit scripts. 




A 

R 

T 

S 


0 

->1 

— >2 

— >3 

— >4 

M 

1 

\ 

1 

\ 

— >2 

\ 

^3 

\ 

— >4 

A 

2 

\ 

1 

\ 

— >2 

\ 

^3 

\ 

— >4 

T 

3 

2 

\ 

2 

\ 

2 

^3 

H 

4 

3 

V 

3 

V 

3 

\ 

3 

S 

5 

4 

V 

4 

V 

4 

\ 

3 


Solution path 2: 

1+ 1+0+ 1+ 0 =3 

S S M D M 

~"M A T H S 

A R T _ S 
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A 

R 

T 

S 


0 

->1 

— >2 

— >3 

— >4 

M 

7 

1 

\ 

1 

\ 

— >2 

\ 

^3 

\ 

— >4 

A 

7 

2 

\ 

1 

\ 

— >2 

\ 

^3 

\ 

— >4 

T 

7 

3 

7 

2 

\ 

2 

\ 

2 

^3 

H 

7 

4 

7 

3 

V 

3 

V 

3 

\ 

3 

S 

7 

5 

7 

4 

V 

4 

V 

4 

\ 

3 


Solution path 3: 

1 + 0 + 1 + 0 + 1 + 0 =3 

D M 1 M D M 

~"M A 7 T H S 

ARTS 


6.3.4 Analysis of DP Edit Distance 

There are 0(ri 2 ) entries in the matrix. Each entry E(i, j) takes 0(1 ) time to compute. The total running 
time is 0(n 2 ). 


6.4 Chain Matrix Multiply 

Suppose we wish to multiply a series of matrices: 

Ai A 2 . . . A n 

In what order should the multiplication be done? A p x q matrix A can be multiplied with a q x r matrix 
B. The result will be a p x r matrix C. In particular, for 1 < i < p and 1 < j < r, 

q 

C[i,j] = ^A[i,k]B[k,j] 

k=1 

There are (p ■ r) total entries in C and each takes 0(q) to compute. 

q 

C[i,j]=^A[i,k]B[k,j] 

k=1 

Thus the total number of multiplications is p ■ q ■ r. Consider the case of 3 matrices: A-| is 5 x 4, A 2 is 
4x6 and A3 is 6 x 2 The multiplication can be carried out either as ((A-| A 2 )A3) or (Ai (A 2 A3)). The 
cost of the two is 


((A 1 A 2 )A 3 ) = (5 ■ 4 • 6) + (5 ■ 6 • 2)= 180 
(At (A 2 A 3 )) = (4 • 6 ■ 2) + (5 ■ 4 ■ 2) =88 
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There is considerable savings achieved even for this simple example. In general, however, in what order 
should we multiply a series of matrices Ai A2 . . . A n . Matrix multiplication is an associative but not 
commutative operation. We are free to add parenthesis the above multiplication but the order of matrices 
can not be changed. The Chain Matrix Multiplication Problem is stated as follows: 

Given a sequence A] , A 2 , . . . , A n and dimensions p 0 , Pi , . . . , p n where A^ is of dimension 
Pi_i x pt, determine the order of multiplication that minimizes the number of operations. 

We could write a procedure that tries all possible parenthesizations. Unfortunately, the number of ways 
of parenthesizing an expression is very large. If there are n items, there are n — 1 ways in which outer 
most pair of parentheses can placed. 



(Ai)(A 2 A 3 A 4 . 

..A n ) 

or 

(Ai A 2 )(A3A 4 . 

..A n ) 

or 

(A 1 A 2 A 3 )(A 4 . 

..A n ) 

or 

(Ai A 2 A3A 4 . . 

• A n _i ) (A n 


Once we split just after the k th matrix, we create two sublists to be parethesized, one with k and other 
with n — k matrices. 


(Ai A 2 . . . Ak) (Ak+i . . . A n ) 


We could consider all the ways of parenthesizing these two. Since these are independent choices, if there 
are L ways of parenthesizing the left sublist and R ways to parenthesize the right sublist, then the total is 
L ■ R. This suggests the following recurrence for P(n), the number of different ways of parenthesizing n 
items: 

1 if n = 1 , 

^k=i P(k)P(n — k) ifn>2 

This is related to a famous function in combinatorics called the Catalan numbers. Catalan numbers are 
related the number of different binary trees on n nodes. Catalan number is given by the formula: 


P(n) = 


C(n) = 


1 


n + 1 


2 n 


n 


In particular, P(n) = C(n — 1 ) C(n) € Ll[4 n /n 3 ^ 2 ) The dominating term is the exponential 4 n thus 
P(n) will grow large very quickly. So this approach is not practical. 


6.4.1 Chain Matrix Multiplication-Dynamic Programming Formulation 

The dynamic programming solution involves breaking up the problem into subproblems whose solutions 
can be combined to solve the global problem. Let At..j be the result of multiplying matrices i through j. It 
is easy to see that At..j is a x pj matrix. 

A3 A 4 A 5 Ag = A3 g 

4x5 5x2 2x8 8x7 4x7 
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At the highest level of parenthesization, we multiply two matrices 

Ai.. n = Ai.j< ■ A k+1 .. n 1 < k < n — 1 . 

The question now is: what is the optimum value of k for the split and how do we parenthesis the 
sub-chains Ai.. k and A k+1 n . We can not use divide and conquer because we do not know what is the 
optimum k. We will have to consider all possible values of k and take the best of them. We will apply 
this strategy to solve the subproblems optimally. 

We will store the solutions to the subproblem in a table and build the table bottom-up (why)? For 
1 < i < j < n, let m[i, j] denote the minimum number of multiplications needed to compute A k .j. The 
optimum can be described by the following recursive formulation. 

• If i = j, there is only one matrix and thus m[i, I] = 0 (the diagonal entries). 

• If i < j, the we are asking for the product At..j. 

• This can be split by considering each k, i < k < j, as A k . k times A k+k .j. 

The optimum time to compute A k k is m[i, k] and optimum time for A k+ i 0 - is in m[k + 1 , j]. Since A k . k 
is a pt_i x p k matrix and A k+ i j is p k x p, matrix, the time to multiply them is xp k xpj. This 
suggests the following recursive rule: 


m[i, i] = 0 

m[i,j] = min (m[i,k] +m[k+ 1,j] +Pt-ip k Pj) 

i< k<j 


We do not want to calculate m entries recursively. So how should we proceed? We will fill m along the 
diagonals. Here is how. Set all m[i, i] = 0 using the base condition. Compute cost for multiplication of a 
sequence of 2 matrices. These are m[l ,2], m[2,3], m[3,4], . . . , m[n — 1 ,n]. m[l ,2], for example is 


m[1,2] =m[l,l] + m[2, 2] + p 0 -pi • p 2 
For example, for m for product of 5 matrices at this stage would be: 


m[l, 1] 

<— m[l , 2] 





m[2,2] 

<— m[2, 3] 





m[3, 3] 

<— m[3,4] 





m[4,4] 

<— m[4, 5] 





m[5,5] 
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Next, we compute cost of multiplication for sequences of three matrices. These are 
m[1 , 3], m[2,4], m[3,5], . . . , m[n — 2, n] . m[l ,3], for example is 


m[l , 3] = min 


m[l,l] + m[2, 3] + p 0 -pi • P 3 
m[l , 2] + m[3, 3] + p 0 • P 2 • P 3 


We repeat the process for sequences of four, five and higher number of matrices. The final result will end 
up in m[1 , n]. 

Example: Let us go through an example. We want to find the optimal multiplication order for 


Ai • A 2 ■ A3 ■ A4 ■ A5 

(5x4) (4x6) (6x2) (2x7) (7x3) 

We will compute the entries of the m matrix starting with the base condition. We first fill that main 
diagonal: 


0 






0 






0 






0 






0 


Next, we compute the entries in the first super diagonal, i.e., the diagonal above the main diagonal: 

m[l ,2] = m[l , 1] + m[2, 2} + p 0 ■ pi • P 2 = 0 + 0 + 5 ■ 4 ■ 6 = 120 

m[2, 3] = m[2, 2] + m[3, 3] + pi ■ p 2 • P 3 = 0 + 0 + 4 ■ 6 • 2 = 48 

m[3,4] = m[3, 3] + m[4,4] + p 2 • P 3 • P 4 = 0 + 0 + 6 ■ 2 • 7 = 84 

m[4, 5] = m[4, 4] + m[5, 5] + p 3 ■ p 4 ■ p 5 = 0 + 0 + 2 • 7 ■ 3 = 42 

The matrix m now looks as follows: 


0 

120 





0 

48 





0 

84 





0 

42 





0 


We now proceed to the second super diagonal. This time, however, we will need to try two possible 
values for k. For example, there are two possible splits for computing m[l , 3]; we will choose the split 
that yields the minimum: 


m[1 , 3] = m[1 , 1] + m[2,3] + po • Pi - Vs == 0 + 48 + 5- 4- 2 = 88 
m[1 , 3] = m[l ,2] + m[3,3] + Po ■ P2 • P3 = 120 + 0 + 5 • 6 ■ 2 = 180 
the minimum m[l ,3] =88 occurs with k = 1 
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Similarly, for m[2,4] and m[3,5]: 

m[2, 4] = m[2, 2] + m[3, 4] + p! • p 2 ■ P4 = 0 + 84 + 4 • 6 ■ 7 = 252 
m[2,4] = m[2, 3] + m[4,4] + pi • p 3 ■ p 4 = 48 + 0 + 4 ■ 2 ■ 7 = 104 
minimum m[2, 4] = 1 04 at k = 3 


m[3, 5] = m[3, 3] + m[4, 5] + p 2 • P3 ■ ps = 0 + 42 + 6 ■ 2 • 3 = 78 
m[3, 5] = m[3,4] + m[5,5] + p 2 • P4 • Ps = 84 + 0 + 6 ■ 7 ■ 3 = 210 
minimum m[3, 5] = 78 at k = 3 


With the second super diagonal computed, the m matrix looks as follow: 


0 

120 

88 




0 

48 

104 




0 

84 

78 




0 

42 





0 


We repeat the process for the remaining diagonals. However, the number of possible splits (values of k) 
increases: 


m[l ,4] = m[1 , 1] + m[2,4] + po • pi • P4 = 0 + 1 04 + 5- 4- 7 = 244 
m[l ,4] = m[l ,2] + m[3,4] + p 0 • p 2 • P4 = 120 + 84 + 5 • 6 • 7 = 414 
m[1 ,4] = m[l , 3] + m[4,4] + p 0 • P3 • P4 = 88 + 0 + 5 • 2 • 7 = 158 
minimum m[1 , 4] = 1 58 at k = 3 


m[2,5] = m[2,2] + m[3,5] +p] ■ p 2 • p 5 = 0 + 78 + 4 ■ 6 ■ 3 = 150 

m[2,5] = m[2,3] + m[4,5] + pi ■ p 3 • Ps = 48 + 42 + 4 ■ 2 ■ 3 = 1 14 

m[2,5] = m[2,4] + m[5,5] + pi ■ P4 • Ps = 104 + 0 + 4- 7- 3 = 188 

minimum m[2, 5] = 1 14 at k = 3 


The matrix m at this stage is: 
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0 

120 

88 

158 



0 

48 

104 

114 



0 

84 

78 




0 

42 





0 


That leaves the m[1 , 5] which can now be computed: 


m[l 

,5] 

= m[l,l] 

+ m[2, 5] 

+ Po ' 

Pi 

■P5 = 

= 0 + 114 + 5- 

m[l 

,5] 

= m[l , 2] 

+ m[3,5] 

+ Po • 

P2 

■P5 = 

= 120 + 78 + 5 

m[l 

,5] 

= m[l , 3] 

+ m[4,5] 

T Po • 

P3 

■P5 = 

= 88 + 42 + 5- 

m[l 

,5] 

= m[l , 4] 

+ m[5,5] 

+ Po ' 

' P4 

■P5 = 

= 158 + 0 + 5- 


minimum : 

m[l , 5] = 

160 at 

k = 

= 3 



We thus have the final cost matrix. 


0 

120 

88 

158 

160 

0 

0 

48 

104 

114 

0 

0 

0 

84 

78 

0 

0 

0 

0 

42 

0 

0 

0 

0 

0 


Here is the order in which m entries are calculated 


0 

1 

5 

8 

10 

0 

0 

2 

6 

9 

0 

0 

0 

3 

7 

0 

0 

0 

0 

4 

0 

0 

0 

0 

0 


and the split k values that led to a minimum m[i, j] value 


0 

1 

1 

3 

3 


0 

2 

3 

3 



0 

3 

3 




0 

4 





0 


. = 174 
3 =288 
. = 160 
. = 263 
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Based on the computation, the minimum cost for multiplying the five matrices is 160 and the optimal 
order for multiplication is 


( (At ( A 2 A 3 ) ) ( A 4 A 5 ) ) 


This can be represented as a binary tree 



Ai 



A 4 a 5 


a 2 a 3 


Figure 6.2: Optimum matrix multiplication order for the five matrices example 


Here is the dynamic programming based algorithm for computing the minimum cost of chain matrix 
multiplication. 

MATRIX-CHAIN(p, N) 

1 for i = 1 , N 

2 do m[i, i] <— 0 

3 for L = 2, N 

4 do 

5 for i = 1 , ri — L + 1 

6 do j i — i T L — 1 

7 m[i, j] 00 

8 for k = 1 , j — 1 

9 do t <— m[i, k] + m[k + 1 , j] + Vi-i • Pk • Vj 

10 if (t < m[i, j]) 

11 then m[i, j] t; s[i, j] <— k 


Analysis: There are three nested loops. Each loop executes a maximum n times. Total time is thus 
0(n 3 ). 

The s matrix stores the values k. The s matrix can be used to extracting the order in which matrices are to 
be multiplied. Here is the algorithm that caries out the matrix multiplication to compute At p 
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MULTIPLY 


1 

if (i = 

)) 

2 

then 

return A[i] 

3 

else 

k < s[i, j] 

4 


X <— MULTIPLY (i, k) 

5 


Y <— MULTIPLY (k+ l,j) 

6 


return X ■ Y 


6.5 0/1 Knapsack Problem 

A thief goes into a jewelry store to steal jewelry items. He has a knapsack (a bag) that he would like to 
fill up. The bag has a limit on the total weight of the objects placed in it. If the total weight exceeds the 
limit, the bag would tear open. The value of of the jewelry items varies for cheap to expensive. The 
thief’s goal is to put items in the bag such that the value of the items is maximized and the weight of the 
items does not exceed the weight limit of the bag. Another limitation is that an item can either be put in 
the bag or not - fractional items are not allowed. The problem is: what jewelry should the thief choose 
that satisfy the constraints? 

Formally, the problem can be stated as follows: Given a knapsack with maximum capacity W, and a set S 
consisting of ri items Each item i has some weight w; and value value Vi (all wi , Vi and W are integer 
values) How to pack the knapsack to achieve maximum total value of packed items? For example, 
consider the following scenario: 



Figure 6.3: Knapsack can hold W — 20 

The knapsack problem belongs to the domain of optimization problems. Mathematically, the problem is 

maximize v t 

ieT 

subject to wt < W 

ieT 
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The problem is called a “0-1 ” problem, because each item must be entirely accepted or rejected. How do 
we solve the problem. We could try the brute-force solution: 


• Since there are ri items, there are 2 n possible combinations of the items (an item either chosen or 
not). 

• We go through all combinations and find the one with the most total value and with total weight 
less or equal to W 


Clearly, the running time of such a brute-force algorithm will be 0(2 n ). Can we do better? The answer is 
“yes”, with an algorithm based on dynamic programming Let us recap the steps in the dynamic 
programming strategy: 


1. Simple Subproblems: We should be able to break the original problem to smaller subproblems 
that have the same structure 

2. Principle of Optimality: Recursively define the value of an optimal solution. Express the solution 
of the original problem in terms of optimal solutions for smaller problems. 

3. Bottom-up computation: Compute the value of an optimal solution in a bottom-up fashion by 
using a table structure. 

4. Construction of optimal solution: Construct an optimal solution from computed information. 


Let us try this: If items are labelled 1 , 2, . . . , n, then a subproblem would be to find an optimal solution 
for 

S k = items labelled 1 , 2, . . . , k 

This is a valid subproblem definition. The question is: can we describe the final solution S n in terms of 
subproblems Sk? Unfortunately, we cannot do that. Here is why. Consider the optimal solution if we can 
choose items 1 through 4 only. 


Solution S 4 

Item 

Wt 

Vi 


1 

2 

3 

• Items chosen are 1 , 2, 3, 4 

2 

3 

4 

• Total weight: 2 + 3 + 4 + 5 = 14 

3 

4 

4 

5 

5 

8 

• Total value: 3 + 4 + 5 + 8 = 20 

5 

9 

10 


Now consider the optimal solution when items 1 through 5 are available. 
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Solution S5 

• Items chosen are 1 , 3, 4, 5 

• Total weight: 2 + 4 + 5 + 9 = 20 

• Total value: 3 + 5 + 8 + 10 = 26 

S 4 is not part of S 5 !! 


Item 

Wt 

Vi 

1 

2 

3 

2 

3 

4 

3 

4 

5 

4 

5 

8 

5 

9 

10 


The solution for S4 is not part of the solution for S5. So our definition of a subproblem is flawed and we 
need another one. 


6.5.1 0/1 Knapsack Problem: Dynamic Programming Approach 

For each i < n and each w < W, solve the knapsack problem for the first i objects when the capacity is 
w. Why will this work? Because solutions to larger subproblems can be built up easily from solutions to 
smaller ones. We construct a matrix V[0 . . . n, 0 . . . W]. For 1 < i < n, and 0 < j < W, V[i, j] will store 
the maximum value of any set of objects {1 , 2, . . . , i} that can fit into a knapsack of weight j. V[n, W] 
will contain the maximum value of all n objects that can fit into the entire knapsack of weight W. 

To compute entries of V we will imply an inductive approach. As a basis, V[0, j] =0 for 0 < j < W 
since if we have no items then we have no value. We consider two cases: 


Leave object i: If we choose to not take object i, then the optimal value will come about by considering 
how to fill a knapsack of size j with the remaining objects {1 , 2, . . . , i — 1 }. This is just V[i — 1 , j] . 

Take object i: If we take object i, then we gain a value of Vt. But we use up w ( of our capacity. With the 
remaining j — Wj capacity in the knapsack, we can fill it in the best possible way with objects 
{1 , 2, . . . , i — 1}. This is V| + V[i — 1 , j — wj. This is only possible if w ( < j. 

This leads to the following recursive formulation: 

V[i, j] = —00 if j < 0 
V[ 0 J ]=0 if j > 0 

VR -i = / _ 1 « A if Wi > j 

|max{V[i- l,j], v t + V[i - 1 , j - wj} ifw t <j 

A naive evaluation of this recursive definition is exponential. So, as usual, we avoid re-computation by 
making a table. 

Example: The maximum weight the knapsack can hold is W is 11. There are five items to choose from. 
Their weights and values are presented in the following table: 
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Weight limit (j): 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

Wi = 1 Vi = 1 

w 2 = 2 V2 = 6 
w 3 = 5 v 3 = 1 8 
w 4 = 6 v 4 = 22 
W5 = 7 V5 = 28 














The [i, j] entry here will be V[i, j], the best value obtainable using the first i rows of items if the 
maximum capacity were j . We begin by initializating and first row. 


Weight limit: 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

W| = 1 Vi = 1 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

w 2 = 2 V2 = 6 

0 












w 3 = 5 v 3 = 1 8 

0 












W4 = 6 V4 = 22 

0 












w 5 = 7 V5 = 28 

0 













Recall that we take V[i, j] to be 0 if either i or j is < 0 . We then proceed to fill in top-down, left-to-right 
always using 


V[i, j] = max{V[i- l,j], + V[i - 1,j -wj} 


Weight limit: 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

Wi = 1 Vi = 1 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

w 2 = 2 v 2 = 6 

0 

1 

6 

7 

7 

7 

7 

7 

7 

7 

7 

7 

w 3 = 5 v 3 = 1 8 

0 












w 4 = 6 v 4 = 22 

0 












W5 = 7 V5 = 28 

0 













Weight limit: 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

Wi = 1 Vi = 1 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

w 2 = 2 v 2 = 6 

0 

1 

6 

7 

7 

7 

7 

7 

7 

7 

7 

7 

w 3 = 5 v 3 = 1 8 

0 

1 

6 

7 

7 

18 

19 

24 

25 

25 

25 

25 

w 4 — 6 v 4 — 22 

0 












W5 = 7 v 5 = 28 

0 













As an illustration, the value of V[ 3 , 7 ] was computed as follows: 

V[ 3 , 7 ] =max{V [3 - 1 , 7 ], v 3 + V [3 - 1 , 7 - w 3 ]} 
=max{V[ 2 , 7 ], 18 + V[ 2 , 7 — 5 ]} 

=max{ 7 , 18 + 6} 

=24 
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Weight limit: 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

Wi = 1 Vi = 1 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

w 2 = 2 V 2 = 6 

0 

1 

6 

7 

7 

7 

7 

7 

7 

7 

7 

7 

W 3 = 5 V 3 = 18 

0 

1 

6 

7 

7 

18 

19 

24 

25 

25 

25 

25 

w 4 = 6 v 4 = 22 
W 5 = 7 V 5 = 28 

0 

0 

1 

6 

7 

7 

18 

22 

24 

28 

29 

29 

40 


Finally, we have 


Weight limit: 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

Wi = 1 Vi = 1 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

W 2 = 2 V2 = 6 

0 

1 

6 

7 

7 

7 

7 

7 

7 

7 

7 

7 

W 3 = 5 V 3 = 18 

0 

1 

6 

7 

7 

18 

19 

24 

25 

25 

25 

25 

w 4 = 6 v 4 = 22 

0 

1 

6 

7 

7 

18 

22 

24 

28 

29 

29 

40 

W 5 = 7 V 5 = 28 

0 

1 

6 

7 

7 

18 

22 

28 

29 

34 

35 

40 


The maximum value of items in the knapsack is 40, the bottom-right entry). The dynamic programming 
approach can now be coded as the following algorithm: 


KNAPSACK(n, W) 

1 forw = 0,W 

2 do V[0, w] <— 0 

3 fori = 0,ri 

4 do V[i, 0] «- 0 

5 forw = 0,W 

6 do if (wi < w & Vi + V[i — 1 , w — wj > V[i — 1 , w] ) 

7 then V[i, w] <— Vt + V[i — 1 , w — wj 

8 else V[i, w] <— V[i — 1 , w] 


The time complexity is clearly 0(n • W). It must be cautioned that as n and W get large, both time and 
space complexity become significant. 

Constructing the Optimal Solution 

The algorithm for computing V[i, j] does not keep record of which subset of items gives the optimal 
solution. To compute the actual subset, we can add an auxiliary boolean array keep[i, j] which is 1 if we 
decide to take the i tH item and 0 otherwise. We will use all the values keep[i, j] to determine the optimal 
subset T of items to put in the knapsack as follows: 

• If keep[n, W] is 1, then n e T. We can now repeat this argument for keep[n — 1 , W — wj. 

• If kee[n, W] is 0, the n ^ T and we repeat the argument for keep[n — 1 , W]. 
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We will add this to the knapsack algorithm: 


KNAPSACK(n, W) 

1 forw = 0,W 

2 do V[0, w] <— 0 

3 fori = 0,ri 

4 do V[i, 0] <- 0 

5 forw = 0,W 

6 do if (wi < w & Vt + V[i — 1 , w — wj > V[i — 1 , w] ) 

7 then V[i, w] <— Vt + V[i — 1 , w — wj; keep [i, w] <— 1 

8 else V[i,w] <— V[i— l,w]; keep[i, w] <— 0 

9 // output the selected items 

10 k <— W 

1 1 for i = n downto 1 

12 do if (keep[i, k] = 1 ) 

13 then output i 

14 k <— k — wt 


Here is the keep matrix for the example problem. 


Weight limit: 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

Wi = 1 Vi = 1 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

W2 = 2 V2 = 6 

0 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

w 3 = 5 v 3 = 18 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

W4 = 6 V4 = 22 

0 

0 

0 

0 

0 

0 

1 

0 

1 

1 

1 

E 

w 5 = 7 v 5 = 28 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

0 


When the item selection algorithm is applied, the selected items are 4 and 3. This is indicated by the 
boxed entries in the table above. 




Chapter 7 

Greedy Algorithms 


An optimization problem is one in which you want to find, not just a solution, but the best solution. 
Search techniques look at many possible solutions. E.g. dynamic programming or backtrack search. A “ 
greedy algorithm” sometimes works well for optimization problems 

A greedy algorithm works in phases. At each phase: 

• You take the best you can get right now, without regard for future consequences. 

• You hope that by choosing a local optimum at each step, you will end up at a global optimum. 

For some problems, greedy approach always gets optimum. For others, greedy finds good, but not always 
best. If so, it is called a greedy heuristic, or approximation. For still others, greedy approach can do very 
poorly. 


7.1 Example: Counting Money 

Suppose you want to count out a certain amount of money, using the fewest possible bills (notes) and 
coins. A greedy algorithm to do this would be: at each step, take the largest possible note or coin that 
does not overshoot. 


while (N > 0){ 

give largest denomination coin < N 
reduce N by value of that coin 

} 

Consider the currency in U.S.A. There are paper notes for one dollar, five dollars, ten dollars, twenty 
dollars, fifty dollars and hundred dollars. The notes are also called “bills”. The coins are one cent, five 
cents (called a “nickle”), ten cents (called a “dime”) and twenty five cents (a “quarter”). In Pakistan, the 
currency notes are five rupees, ten rupees, fifty rupees, hundred rupees, five hundred rupees and thousand 
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rupees. The coins are one rupee and two rupees. Suppose you are asked to give change of $6.39 (six 
dollars and thirty nine cents), you can choose: 

• a $5 note 

• a $1 note to make $6 

• a 25 cents coin (quarter), to make $6.25 

• a 10 cents coin (dime), to make $6.35 

• four 1 cents coins, to make $6.39 


Notice how we started with the highest note, $5, before moving to the next lower denomination. 

Formally, the Coin Change problem is: Given k denominations di , d 2 , . . . , dk and given N, find a way of 
writing 

N = ii di + i 2 d 2 + • • • + ikdk 


such that 


i i + i 2 + • • • + ik is minimized. 


The “size” of problem is k. 

The greedy strategy works for the coin change problem but not always. Here is an example where it fails. 
Suppose, in some (fictional) monetary system, “ krons” come in 1 kron, 7 kron, and 10 kron coins Using 
a greedy algorithm to count out 15 krons, you would get A 10 kron piece Five 1 kron pieces, for a total of 
15 krons This requires six coins. A better solution, however, would be to use two 7 kron pieces and one 1 
kron piece This only requires three coins The greedy algorithm results in a solution, but not in an optimal 
solution 


The greedy approach gives us an optimal solution when the coins are all powers of a fixed denomination. 


N = i 0 D° + GD 1 + i 2 D 2 + • • • + i k D k 


Note that this is N represented in based D. U.S.A coins are multiples of 5: 5 cents, 10 cents and 25 cents. 


7.1.1 Making Change: Dynamic Programming Solution 

The general coin change problem can be solved using Dynamic Programming. Set up a Table, 

C [1 ..k, 0..N] in which the rows denote available denominations, di; (1 < i < k) and columns denote the 
amount from 0 . . . N units, (0 < j < N). C [i, j] denotes the minimum number of coins, required to pay 
an amount j using only coins of denominations 1 to i. C[k, N] is the solution required. 

To pay an amount j units, using coins of denominations 1 to i, we have two choices: 

1. either chose NOT to use any coins of denomination i, 

2. or chose at least one coin of denomination i, and also pay the amount (j — d|). 
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To pay (j — di) units it takes C[i, j — dj coins. Thus, 

C[i, j] = 1 + C[i, j - dj 

Since we want to minimize the number of coins used, 

C[i, j] = min(C[i — 1, j], 1 + C[i, j - dj) 

Here is the dynamic programming based algorithm for the coin change problem. 


coins(N) 

1 d[l ..rx] ={1,4,6} // ( coinage, for example ) 

2 for i = 1 to k 

3 do c[i, 0] < — 0 

4 for i = 1 to k 

5 do for j = 1 to N 

6 do if (i = 1 & j < d[i]) 

7 then c[i, j] <— oo 

8 else if (i = 1 ) 


9 

then c[i, j] <— 1 + 

c[l > j d[l]] 

10 

else if (j < d[i] ) 


11 

then c[i , )] 

<— c[i — 1 , j] 

12 

else c[i, j] 

<— min (c[i 

13 

return c[k, N] 



,)],!+ c[i, j - d[i]]) 


7.1.2 Complexity of Coin Change Algorithm 

Greedy algorithm (non-optimal) takes O(k) time. Dynamic Programming takes O(kN) time. Note that 
N can be as large as 2 k so the dynamic programming algorithm is really exponential in k. 


7.2 Greedy Algorithm: Huffman Encoding 

The Huffman codes provide a method of encoding data efficiently. Normally, when characters are coded 
using standard codes like ASCII. Each character is represented by a fixed-length codeword of bits, e.g., 8 
bits per character. Fixed-length codes are popular because it is very easy to break up a string into its 
individual characters, and to access individual characters and substrings by direct indexing. However, 
fixed-length codes may not be he most efficient from the perspective of minimizing the total quantity of 
data. 

Consider the string “ abacdaacac”. if the string is coded with ASCII codes, the message length would be 
1 0 x 8 = 80 bits. We will see shortly that the same string encoded with a variable length Huffman 
encoding scheme will produce a shorter message. 
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7.2.1 Huffman Encoding Algorithm 


Here is how the Huffman encoding algorithm works. Given a message string, determine the frequency of 
occurrence (relative probability) of each character in the message. This can be done by parsing the 
message and counting how many time each character (or symbol) appears. The probability is the number 
of occurrence of a character divided by the total characters in the message. The frequencies and 
probabilities for the example string “ abacdaacac” are 


character 

a 

b 

c 

d 

frequency 

5 

1 

3 

1 

probability 

0.5 

0.1 

0.3 

0.1 


Next, create binary tree (leaf) node for each symbol (character) that occurs with nonzero frequency Set 
node weight equal to the frequency of the symbol. Now comes the greedy part: Find two nodes with 
smallest frequency. Create a new node with these two nodes as children, and with weight equal to the 
sum of the weights of the two children. Continue until we have a single tree. 

Finding two nodes with the smallest frequency can be done efficiently by placing the nodes in a 
heap-based priority queue. The min-heap is maintained using the frequencies. When a new node is 
created by combining two nodes, the new node is placed in the priority queue. Here is the Huffman tree 
building algorithm. 


huffman(N, symbol[l ..N], freq[l ..N]) 

1 for i = 1 to N 

2 do t <— TreeNode(spmbol[i], freq[i]) 

3 pq.insert(t, freq [i] ) // priority queue 

4 for i = 1 to N — 1 

5 do x = pq.removeQ; y = pq.removeQ 

6 z <— new TreeNode 

7 z.left <— x; z.riglvt <— y 

8 z.freq <— x. freq + y. freq 

9 pq.insert(z, z.freq); 

10 return pq.removeQ; //root 


Figure 7.1 shows the tree built for the example message “abacdaacac” 
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Prefix Property: 

The codewords assigned to characters by the Huffman algorithm have the property that no codeword is a 
prefix of any other: 


character 

a 

b 

c 

d 

frequency 

5 

1 

3 

1 

probability 

0.5 

0.1 

0.3 

0.1 

codeword 

0 

110 

10 

111 


The prefix property is evident by the fact that codewords are leaves of the binary tree. Decoding a prefix 
code is simple. We traverse the root to the leaf letting the input 0 or 1 tell us which branch to take. 

Expected encoding length: 

If a string of rt characters over the alphabet C = {a, b, c, d} is encoded using 8-bit ASCII, the length of 
encoded string is 8n. For example, the string “abacdaacac” will require 8 x 1 0 = 80 bits. The same 
string encoded with Huffman codes will yield 


a 

b 

a 

c 

d 

a 

a 

c 

a 

c 

0 

110 

0 

10 

111 

0 

0 

10 

0 

10 


This is just 17 bits, a significant saving!. For a string of n characters over this alphabet, the expected 
encoded string length is 

n(0.5- 1 +0.1 -3 + 0.3 -2 + 0.1 -3) = 1.7n 

In general, let p(x) be the probability of occurrence of a character, and let d T ( x ) denote the length of the 
codeword relative to some prefix tree T. The expected number of bits needed to encode a text with n 
characters is given by 

B (T) = n y~ p(x)d T (x) 

xgC 
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7.2.2 Huffman Encoding: Correctness 

Huffman algorithm uses a greedy approach to generate a prefix code T that minimizes the expected 
length B(T) of the encoded string. In other words, Huffman algorithm generates an optimum prefix code. 
The question that remains is that why is the algorithm correct? 

Recall that the cost of any encoding tree T is 

B(T) = n y~ p(x)d T (x) 

xSC 

Our approach to prove the correctness of Huffman Encoding will be to show that any tree that differs 
from the one constructed by Huffman algorithm can be converted into one that is equal to Huffman’s tree 
without increasing its costs. Note that the binary tree constructed by Huffman algorithm is a full binary 
tree. 

Claim: 

Consider two characters x and y with the smallest probabilities. Then there is optimal code tree in which 
these two characters are siblings at the maximum depth in the tree. 

Proof: 

Let T be any optimal prefix code tree with two siblings b and c at the maximum depth of the tree. Such a 
tree is shown in Figure 7.2Assume without loss of generality that 

p(b)<p(c) and p(x) < p(y) 



Figure 7.2: Optimal prefix code tree T 

Since x and y have the two smallest probabilities (we claimed this), it follows that 


p(x) < p(b) and p(y) < p(c) 
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Since b and c are at the deepest level of the tree, we know that 


d(b) > d(x) and d(c) > d(y) (d is the depth) 


Thus we have 


p(b) — p(x) > 0 


and 


d(b) — d(x) > 0 


Hence their product is non-negative. That is, 


(p(b) — p(x)) • (d(b) - d(x)) > 0 


Now swap the positions of x and b in the tree 



Figure 7.3: Swap x and b in tree prefix tree T 


This results in a new tree T ; 
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Figure 7.4: Prefix tree T' after x and b are swapped 

Let’s see how the cost changes. The cost of T' is 

B(T') = B(T) -p(x)d(x) + p(x)d(b) — p(b)d(b) + p(b)d(x) 

= B(T) + p(x)[d(b) - d(x)] -p(b)[d(b) - d(x)] 

= B(T) - (p(b) — p(x))(d(b) - d(x)) 

< B(T) because (p(b) — p(x))(d(b) — d(x)) > 0 

Thus the cost does not increase, implying that T' is an optimal tree. 

By switching y with c we get the tree T". Using a similar argument, we can show that T" is also optimal. 



The final tree T" satisfies the claim we made earlier, i.e., consider two characters x and y with the 
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smallest probabilities. Then there is optimal code tree in which these two characters are siblings at the 
maximum depth in the tree. 

The claim we just proved asserts that the first step of Huffman algorithm is the proper one to perform (the 
greedy step). The complete proof of correctness for Huffman algorithm follows by induction on n. 

Claim: Huffman algorithm produces the optimal prefix code tree. 

Proof: The proof is by induction on n, the number of characters. For the basis case, n = 1 , the tree 
consists of a single leaf node, which is obviously optimal. We want to show it is true with exactly n 
characters. 

Suppose we have exactly n characters. The previous claim states that two characters x and y with the 
lowest probability will be siblings at the lowest level of the tree. Remove x and y and replace them with a 
new character z whose probability is p(z) = p(x) + p(y). Thus n — 1 characters remain. 

Consider any prefix code tree T made with this new set of n — 1 characters. We can convert T into prefix 
code tree T' for the original set of n characters by replacing z with nodes x and y. This is essentially 
undoing the operation where x and y were removed an replaced by z. The cost of the new tree T' is 

B(T') = B(T) — p(z)d(z) +p(x)[d(z) + 1] +p(y)[d(z) + 1] 

= B(T) - (p(x) +p(y))d(z) + (p(x) +p(y))[d(z) + 1] 

= B(T) + (p(x) +p(y))[d(z) + 1 - d(z)] 

= B(T) T p (x) T p(tj) 


The cost changes but the change depends in no way on the structure of the tree T (T is for n — 1 
characters). Therefore, to minimize the cost of the final tree T', we need to build the tree T on n — 1 
characters optimally. By induction, this is exactly what Huffman algorithm does. Thus the final tree is 
optimal. 


7.3 Activity Selection 

The activity scheduling is a simple scheduling problem for which the greedy algorithm approach provides 
an optimal solution. We are given a set S = {ai , a 2 , . . . , a n } of n activities that are to be scheduled to use 
some resource. Each activity a L must be started at a given start time s L and ends at a given finish time f t . 

An example is that a number of lectures are to be given in a single lecture hall. The start and end times 
have be set up in advance. The lectures are to be scheduled. There is only one resource (e.g., lecture hall) 
Some start and finish times may overlap. Therefore, not all requests can be honored. We say that two 
activities at and cq are non-interfering if their start-finish intervals do not overlap. I.e, 

(si, ft) fl (sj, fj) = 0. The activity selection problem is to select a maximum-size set of mutually 
non-interfering activities for use of the resource. 

So how do we schedule the largest number of activities on the resource? Intuitively, we do not like long 
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activities Because they occupy the resource and keep us from honoring other requests. This suggests the 
greedy strategy: Repeatedly select the activity with the smallest duration (ft — s0 and schedule it, 
provided that it does not interfere with any previously scheduled activities. Unfortunately, this turns out 
to be non-optimal 

Here is a simple greedy algorithm that works: Sort the activities by their finish times. Select the activity 
that finishes first and schedule it. Then, among all activities that do not interfere with this first job, 
schedule the one that finishes first, and so on. 


SCHEDULE(a[l . . M ] ) 

1 sort a[l ..N] by finish times 

2 A <— {a[l ]}; // schedule activity 1 first 

3 prev < — 1; II most recently scheduled 

4 for i = 2 to N 

5 do if (a[i]. start > a[prev], finish.) 

6 then A <— A U a[i]; prev <— i 


Figure 7.5 shows an example of the activity scheduling algorithm. There are eight activities to be 
scheduled. Each is represented by a rectangle. The width of a rectangle indicates the duration of an 
activity. The eight activities are sorted by their finish times. The eight rectangles are arranged to show the 
sorted order. Activity a i is scheduled first. Activities a 2 and a 3 interfere with so they ar not selected. 
The next to be selected is <14. Activities a 5 and a 6 interfere with a 4 so are not chosen. The last one to be 
chosen is 0.7. Eventually, only three out of the eight are scheduled. 

Timing analysis: Time is dominated by sorting of the activities by finish times. Thus the complexity is 
O(NlogN). 
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7.3.1 Correctness of Greedy Activity Selection 


Our proof of correctness is based on showing that the first choice made by the algorithm is the best 
possible. And then using induction to show that the algorithm is globally optimal. The proof structure is 
noteworthy because many greedy correctness proofs are based on the same idea: Show that any other 
solution can be converted into the greedy solution without increasing the cost. 

Claim: 

Let S = {di , a. 2 , . . . , a n } of n activities, sorted by increasing finish times, that are to be scheduled to use 
some resource. Then there is an optimal schedule in which activity cii is scheduled first. 

Proof: 

Let A be an optimal schedule. Let x be the activity in A with the smallest finish time. If x = ai then we 
are done. Otherwise, we form a new schedule A' by replacing x with activity ai . 
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Figure 7.6: Activity X = ai 


We claim that A' = A — {x} U {a-|} is a feasible schedule, i.e., it has no interfering activities. This 
because A — {x} cannot have any other activities that start before x finishes, since otherwise, these 
activities will interfere with x. 



Since ai is by definition the first activity to finish, it has an earlier finish time than x. Thus ai cannot 
interfere with any of the activities in A — {x}. Thus, A' is a feasible schedule. Clearly A and A' contain 
the same number of activities implying that A' is also optimal. 


Claim: 
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The greedy algorithm gives an optimal solution to the activity scheduling problem. 

Proof: 

The proof is by induction on the number of activities. For the basis case, if there are no activities, then the 
greedy algorithm is trivially optimal. For the induction step, let us assume that the greedy algorithm is 
optimal on any set of activities of size strictly smaller than |S| and we prove the result for S. Let S' be the 
set of activities that do not interfere with activity ai , That is 


S' = {a t e S\si > f i} 


Any solution for S' can be made into a solution for S by simply adding activity ai, and vice versa. 
Activity ai is in the optimal schedule (by the above previous claim). It follows that to produce an optimal 
schedule for the overall problem, we should first schedule ai and then append the optimal schedule for 
S'. But by induction (since |S'| < |S|), this exactly what the greedy algorithm does. 


7.4 Fractional Knapsack Problem 


Earlier we saw the 0-1 knapsack problem. A knapsack can only carry W total weight. There are n items; 
the i th item is worth vt and weighs w L . Items can either be put in the knapsack or not. The goal was to 
maximize the value of items without exceeding the total weight limit of W. In contrast, in the fractional 
knapsack problem, the setup is exactly the same. But, one is allowed to tak e, fraction of an item for a 
fraction of the weight and fraction of value. The 0-1 knapsack problem is hard to solve. However, there is 
a simple and efficient greedy algorithm for the fractional knapsack problem. 

Let pi = vt/wi denote the value per unit weight ratio for item i. Sort the items in decreasing order of Pi. 
Add items in decreasing order of P|. If the item fits, we take it all. At some point there is an item that 
does not fit in the remaining space. We take as much of this item as possible thus filling the knapsack 
completely. 
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Figure 7.8: Greedy solution to the fractional knapsack problem 


It is easy to see that the greedy algorithm is optimal for the fractional knapsack problem. Given a room 
with sacks of gold, silver and bronze, one (thief?) would probably take as much gold as possible. Then 
take as much silver as possible and finally as much bronze as possible. It would never benefit to take a 
little less gold so that one could replace it with an equal weight of bronze. 

We can also observe that the greedy algorithm is not optimal for the 0-1 knapsack problem. Consider 
the example shown in the Figure 7.9. If you were to sort the items by p t , then you would first take the 
items of weight 5, then 20, and then (since the item of weight 40 does not fit) you would settle for the 
item of weight 30, for a total value of $30 + $100 + $90 = $220. On the other hand, if you had been less 
greedy, and ignored the item of weight 5, then you could take the items of weights 20 and 40 for a total 
value of $100+$160 = $260. This is shown in Figure 7.10. 
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knapsack $30 $20 $100 $90 $160 



$220 

Greedy solution 
to 0-1 knapsack 


Figure 7.9: Greedy solution for the 0-1 knapsack problem (non-optimal) 



knapsack $30 $20 $100 $90 $160 

6 2 5 3 4 



$260 

Optimal solution 
to 0-1 knapsack 


Figure 7.10: Optimal solution for the 0-1 knapsack problem 
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Chapter 8 
Graphs 


We begin a major new topic: Graphs. Graphs are important discrete structures because they are a flexible 
mathematical model for many application problems. Any time there is a set of objects and there is some 
sort of “connection” or “relationship” or “interaction” between pairs of objects, a graph is a good way to 
model this. Examples of this can be found in computer and communication networks transportation 
networks, e.g., roads VLSI, logic circuits surface meshes for shape description in computer-aided design 
and GIS precedence constraints in scheduling systems. 

A graph G = (V, E) consists of a finite set of vertices V (or nodes) and E, a binary relation on V called 
edges. E is a set of pairs from V. If a pair is ordered , we have a directed graph. For unordered pair, we 
have an undirected graph. 



Figure 8.1: Types of graphs 


A vertex w is adjacent to vertex v if there is an edge from v to w. 
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In an undirected graph, we say that an edge is incident on a vertex if the vertex is an endpoint of the 
edge, of the edge 



e2 



( 1 r 


{ 3 ) 

el incident on vertices 1 & 2 

el 

\e3 


e2 incident on vertices 1 & 3 

( o \— 


I A 

e3 incident on vertices 1 & 4 

2 

e4 

1 T- 

e4 incident on vertices 2 & 4 


Figure 8.3: Incidence of edges on vertices 


In a digraph, the number of edges coming out of a vertex is called the out-degree of that vertex. Number 
of edges coming in is the in-degree. In an undirected graph, we just talk of degree of a vertex. It is the 
number of edges incident on the vertex. 
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Figure 8.4: In and out degrees of vertices of a graph 


For a digraph G = (V, t), 


y~ in-degree{v) = out-degree[v ) = |E 

vev vev 


where |E| means the cardinality of the set E, i.e., the number of edges. 
For an undirected graph G = (V, E), 


y~ degree{v) = 2|E| 
vev 


where |E| means the cardinality of the set E, i.e., the number of edges. 

A path in a directed graphs is a sequence of vertices (v 0 , Vi , . . . , v k ) such that , v-t) is an edge for 
i = 1 , 2, . . . , k. The length of the paths is the number of edges, k. A vertex w is reachable from vertex 
u is there is a path from u to w. A path is simple if all vertices (except possibly the fist and last) are 
distinct. 

A cycle in a digraph is a path containing at least one edge and for which Vo = V| C . A Hamiltonian cycle 
is a cycle that visits every vertex in a graph exactly once. A Eulerian cycle is a cycle that visits every 
edge of the graph exactly once. There are also “path” versions in which you do not need return to the 
starting vertex. 
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Figure 8.5: Cycles in a directed graph 


A graph is said to be acyclic if it contains no cycles. A graph is connected if every vertex can reach 
every other vertex. A directed graph that is acyclic is called a directed acyclic graph (DAG). 

There are two ways of representing graphs: using an adjacency matrix and using an adjacency list. Let 
G = (V, E) be a digraph with n = |V| and let e = |E|. We will assume that the vertices of G are indexed 
{1,2,..., n}. 


An adjacency matrix is a n x n matrix defined for 1 < v, w < n. 


A[v, w] 


1 if (v, w) £ E 
0 otherwise 


An adjacency list is an array Adj [1 . .n] of pointers where for 1 < v < n, Adj [v] points to a linked list 
containing the vertices which are adjacent to v 

Adjacency matrix requires 0(n 2 ) storage and adjacency list requires 0(n + e) storage. 



Figure 8.6: Graph Representations 


8.1 Graph Traversal 

To motivate our first algorithm on graphs, consider the following problem. We are given an undirected 
graph G = (V, E) and a source vertex s £ V. The length of a path in a graph is the number of edges on 




8.1. GRAPH TRAVERSAL 


117 


the path. We would like to find the shortest path from s to each other vertex in the graph. The final result 
will be represented in the following way. For each vertex v e V, we will store d[v] which is the distance 
(length of the shortest path) from s to v. Note that d[s] = 0. We will also store a predecessor (or parent) 
pointer 7t[v] which is the first vertex along the shortest path if we walk from v backwards to s. We will set 
7t[s] = Nil. 

There is a simple brute-force strategy for computing shortest paths. We could simply start enumerating 
all simple paths starting at s, and keep track of the shortest path arriving at each vertex. However, there 
can be as many as rt! simple paths in a graph. To see this, consider a fully connected graph shown in 
Figure 8.7 



Figure 8.7: Fully connected graph 


There n choices for source node s, (n — 1 ) choices for destination node, (n — 2) for first hop (edge) in 
the path, (n — 3) for second, (n — 4) for third down to (n — (n — 1 )) for last leg. This leads to n! simple 
paths. Clearly this is not feasible. 


8.1.1 Breadth-first Search 

Here is a more efficient algorithm called the breadth-first search (BFS) Start with s and visit its adjacent 
nodes. Label them with distance 1. Now consider the neighbors of neighbors of s. These would be at 
distance 2. Now consider the neighbors of neighbors of neighbors of s. These would be at distance 3. 
Repeat this until no more unvisited neighbors left to visit. The algorithm can be visualized as a wave 
front propagating outwards from s visiting the vertices in bands at ever increasing distances from s. 
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Figure 8.8: Source vertex for breadth-first-search (BFS) 



Figure 8.9: Wave reaching distance 1 vertices during BFS 



Figure 8.10: Wave reaching distance 2 vertices during BFS 
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Figure 8.11: Wave reaching distance 3 vertices during BFS 


8.1.2 Depth-first Search 

Breadth-first search is one instance of a general family of graph traversal algorithms. Traversing a graph 
means visiting every node in the graph. Another traversal strategy is depth-first search (DFS). DFS 
procedure can be written recursively or non-recursively. Both versions are passed s initially. 


recursiveDFS(v) 

1 if (v is unmarked ) 

2 then mark v 

3 for each edge (v,w) 

4 do recursiveDFS(w) 


iterativeDFS(s) 

1 push(s) 

2 while stack not empty 

3 dov<— pop() 

4 if v is unmarked 

5 then mark v 

6 for each edge (v,w) 

7 do push(w) 


8.1.3 Generic Graph Traversal Algorithm 

The generic graph traversal algorithm stores a set of candidate edges in some data structures we’ll call a 
“bag”. The only important properties of the “bag” are that we can put stuff into it and then later take stuff 
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back out. Here is the generic traversal algorithm. 


Traverse(s) 

1 put (0,s) in bag 

2 while bag not empty 

3 do take (p,v) from bag 

4 if (v is unmarked ) 

5 then mark v 

6 parent (v) p 

7 for each edge (v,w) 

8 do put (v,w) in bag 


Notice that we are keeping edges in the bag instead of vertices. This is because we want to remember, 
whenever we visit v for the first time, which previously-visited vertex p put v into the bag. The vertex p 
is call the parent ofv. 

The running time of the traversal algorithm depends on how the graph is represented and what data 
structure is used for the bag. But we can make a few general observations. 

• Since each vertex is visited at most once, the for loop in line 7 is executed at most V times. 

• Each edge is put into the bag exactly twice; once as (u, v) and once as (v, u), so line 8 is executed 
at most 2E times. 

• Finally, since we can’t take out more things out of the bag than we put in, line 3 is executed at most 
2E + 1 times. 

• Assume that the graph is represented by an adjacency list so the overhead of the for loop in line 7 is 
constant per edge. 

If we implement the bag by using a stack , we have depth-first search (DFS) or traversal. 

Traverse(s) 

1 push(0,s) 

2 while stack not empty 

3 dopop(p,v) 

4 if (v is unmarked ) 

5 then mark v 

6 parent (v) <— p 

7 for each edge (v,w) 

8 dopush(v,w) 


Figures 8.12 to 8.20 show a trace of the DFS algorithm applied to a graph. The figures show the content 
of the stack during the execution of the algorithm. 
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Figure 8.12: Trace of Depth-first-search algorithm: source vertex ‘s’ 



Figure 8.13: Trace of DFS algorithm: vertex ‘a’ popped 



Figure 8.14: Trace of DFS algorithm: vertex ‘c’ popped 
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Figure 8.15: Trace of DFS algorithm: vertex ‘f’ popped 



Figure 8.16: Trace of DFS algorithm: vertex ‘g’ popped 
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Figure 8.17: Trace of DFS algorithm: vertex ‘e’ popped 



Figure 8.18: Trace of DFS algorithm: vertex ‘b’ popped 
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Figure 8.19: Trace of DFS algorithm: vertex ‘d’ popped 



Figure 8.20: Trace of DFS algorithm: the final DFS tree 


Each execution of line 3 or line 8 in the Traverse-DFS algorithm takes constant time. So the overall 
running time is 0(V + E). Since the graph is connected, V < E + 1 , this is 0(E). 

If we implement the bag by using a queue , we have breadth-first search (BFS). Each execution of line 3 
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or line 8 still takes constant time. So overall running time is still 0(E). 


Traverse(s) 

1 enqueue (0,s) 

2 while queue not empty 

3 do dequeue (p,v) 

4 if (v is unmarked ) 

5 then mark v 

6 parent (v) p 

7 for each edge (v,w) 

8 do enqueue (v,w) 


If the graph is represented using an adjacency matrix, the finding of all the neighbors of vertex in line 7 
takes 0(V) time. Thus depth-first and breadth-first take 0(V 2 ) time overall. 

Either DFS or BFS yields a spanning tree of the graph. The tree visits every vertex in the graph. This fact 
is established by the following lemma: 

Lemma: 

The generic Traverse(S) marks every vertex in any connected graph exactly once and the set of edges 
(v,parent(v)) with parent(v) / 0 form a spanning tree of the graph. 

Proof: 

First, it should be obvious that no vertex is marked more than once. Clearly, the algorithm marks s. Fet 
v ^ s be a vertex and let s —>•••—> u — > v be a path from s to v with the minimum number of edges. 

Since the graph is connected, such a path always exists. If the algorithm marks u, then it must put (u, v) 
into the bag, so it must take (u, v) out of the bag at which point v must be marked. Thus, by induction on 
the shortest-path distance from s, the algorithm marks every vertex in the graph. 

Call an edge (v, parent(v) ) with parent(v) ^ 0, a parent edge. For any node v, the path of parent 
edges v — > parent(v) — > parent(parent(v)) — > . . . eventually leads back to s. So the set of parent 
edges form a connected graph. 

Clearly, both end points of every parent edge are marked, and the number of edges is exactly one less 
than the number of vertices. Thus, the parent edges form a spanning tree. 


8.1.4 DFS - Timestamp Structure 

As we traverse the graph in DFS order, we will associate two numbers with each vertex. When we first 
discover a vertex u, store a counter in d[u]. When we are finished processing a vertex, we store a counter 
in f[u]. These two numbers are time stamps. 

Consider the recursive version of depth-first traversal 
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DFS(G) 

1 for (each u e V) 

2 do color[u] <— white 

3 pred[u] nil 

4 time <— 0 

5 for each u e V 

6 do if ( color [u] = white) 

7 then DFSvisit(u) 


The DFSvisit routine is as follows: 


DFSvisit(u) 

1 color[u] <— gray; // mark u visited 

2 d[u] <— ++ time 

3 for (each v e Adj[u]) 

4 do if ( color [v] = white) 

5 then pred[v] <— u 

6 DFSvisit(v) 

7 color[u] <— black; // we are done with u 

8 f [u] <— ++ time; 


Figures 8.21 through 8.25 present a trace of the execution of the time stamping algorithm. Terms like 
“2/5” indicate the value of the counter (time). The number before the is the time when a vertex was 
discovered (colored gray) and the number after the is the time when the processing of the vertex 
finished (colored black). 



Figure 8.21: DFS with time stamps: recursive calls initiated at vertex ‘a’ 
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Figure 8.22: DFS with time stamps: processing of ‘b’ and ‘c’ completed 



Figure 8.23: DFS with time stamps: recursive processing of T and ‘g’ 
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Figure 8.24: DFS with time stamps: processing of ‘f’ and ‘g’ completed 


DFS(d) 
DFS(e) 
return e 
return f 


d e 

(11/144 — 412/T3 




Figure 8.25: DFS with time stamps: processing of ‘d’ and 


Notice that the DFS tree structure (actually a collection of trees, or a forest) on the structure of the graph 
is just the recursion tree, where the edge (u, v) arises when processing vertex u we call DFS VlSlT(v) for 
some neighbor v. For directed graphs the edges that are not part of the tree (indicated as dashed edges in 
Figures 8.21 through 8.25) edges of the graph can be classified as follows: 

Back edge: (u, v) where v is an ancestor of u in the tree. 
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Forward edge: (u, v) where v is a proper descendent of u in the tree. 


Cross edge: (u, v) where u and v are not ancestor or descendent of one another. In fact, the edge may 
go between different trees of the forest. 


The ancestor and descendent relation can be nicely inferred by the parenthesis lemma, u is a descendent 
of v if and only if [d[u], f[u]] C [d[v],f[v]]. u is a ancestor of v if and only if [d[u], f[u]] D [d[v],f[v]]. u 
is unrelated to v if and only if [d[u] , f [u]] and [d[v] , f [v]] are disjoint. The is shown in Figure 8.26. The 
width of the rectangle associated with a vertex is equal to the time the vertex was discovered till the time 
the vertex was completely processed (colored black). Imagine an opening parenthesis ‘(’ at the start of 
the rectangle and and closing parenthesis ‘)’ at the end of the rectangle. The rectangle (parentheses) for 
vertex ‘b’ is completely enclosed by the rectangle for ‘a’. Rectangle for ‘c’ is completely enclosed by 
vertex ‘b’ rectangle. 



Figure 8.26: Parenthesis lemma 


Figure 8.27 shows the classification of the non-tree edges based on the parenthesis lemma. Edges are 
labelled ‘F’, ‘B’ and ‘C’ for forward, back and cross edge respectively. 
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Figure 8.27: Classfication of non-tree edges in the DFS tree for a graph 


For undirected graphs, there is no distinction between forward and back edges. By convention they are 
all called back edges. Furthermore, there are no cross edges (can you see why not?) 

8.1.5 DFS - Cycles 

The time stamps given by DFS allow us to determine a number of things about a graph or digraph. For 
example, we can determine whether the graph contains any cycles. We do this with the help of the 
following two lemmas. 

Lemma: Given a digraph G = (V, E), consider any DFS forest of G and consider any edge (u, v) e E. 
If this edge is a tree, forward or cross edge, then f [u] > f [v]. If this edge is a back edge, then 
f[u] < f[v]. 

Proof: For the non-tree forward and back edges the proof follows directly from the parenthesis lemma. 
For example, for a forward edge (u, v), v is a descendent of u and so v’s start-finish interval is 
contained within u’s implying that v has an earlier finish time. For a cross edge (u, v) we know 
that the two time intervals are disjoint. When we were processing u, v was not white (otherwise 
(u, v) would be a tree edge), implying that v was started before u. Because the intervals are 
disjoint, v must have also finished before u. 
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Lemma: Consider a digraph G = (V, E) and any DFS forest for G. G has a cycle if and only if the DFS 
forest has a back edge. 


Proof: If there is a back edge (u, v) then v is an ancestor of u and by following tree edge from v to u, 
we get a cycle. 


We show the contrapositive: suppose there are no back edges. By the lemma above, each of the 
remaining types of edges, tree, forward, and cross all have the property that they go from vertices 
with higher finishing time to vertices with lower finishing time. Thus along any path, finish times 
decrease monotonically, implying there can be no cycle. 


The DFS forest in Figure 8.27 has a back edge from vertex ‘g’ to vertex ‘a’. The cycle is ‘a-g-f’. 

Beware : No back edges means no cycles. But you should not infer that there is some simple relationship 
between the number of back edges and the number of cycles. For example, a DFS tree may only have a 
single back edge, and there may anywhere from one up to an exponential number of simple cycles in the 
graph. 

A similar theorem applies to undirected graphs, and is not hard to prove. 


8.2 Precedence Constraint Graph 


A directed acyclic graph (DAG) arise in many applications where there are precedence or ordering 
constraints. There are a series of tasks to be performed and certain tasks must precede other tasks. For 
example, in construction, you have to build the first floor before the second floor but you can do electrical 
work while doors and windows are being installed. In general, a precedence constraint graph is a DAG in 
which vertices are tasks and the edge (u, v) means that task u must be completed before task v begins. 

For example, consider the sequence followed when one wants to dress up in a suit. One possible order 
and its DAG are shown in Figure 8.28. Figure 8.29 shows the DFS with time stamps of the DAG. 
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Figure 8.28: Order of dressing up in a suit 



Figure 8.29: DFS of dressing up DAG with time stamps 


Another example of precedence constraint graph is the sets of prerequisites for CS courses in a typical 
undergraduate program. 
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Cl 

C2 

C3 

C4 

Introduction to Computers 

Introduction to Computer Programming 
Discrete Mathematics 

Data Structures 

C2 

C5 

C6 

C7 

C8 

Digital Logic Design 

Automata Theory 

Analysis of Algorithms 

Computer Organization and Assembly 

C2 

C3 

C3,C4 

C2 

C9 

CIO 

Cll 

C12 

Data Base Systems 

Computer Architecture 

Computer Graphics 

Software Engineering 

C4, Cl 

C4, C5,C8 

C4,C7 

C7.C11 

C13 

C14 

C15 

Operating System 

Compiler Construction 

Computer Networks 

C4,C7,C1 1 

C4,C6,C8 

C4,C7,C10 


Table 8.1: Prerequisites for CS courses 


The prerequisites can be represented with a precedence constraint graph which is shown in Figure 8.30 



Figure 8.30: Precedence constraint graph for CS courses 


8.3 Topological Sort 

A topological sort of a DAG is a linear ordering of the vertices of the DAG such that for each edge (u, v), 
u appears before v in the ordering. 

Computing a topological ordering is actually quite easy, given a DFS of the DAG. For every edge (u, v) 
in a DAG, the finish time of u is greater than the finish time of v (by the lemma). Thus, it suffices to 
output the vertices in the reverse order of finish times. 
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We run DFS on the DAG and when each vertex is finished, we add it to the front of a linked. Note that in 
general, there may be many legal topological orders for a given DAG. 


TopologicalSort(G) 

1 for (each, u e V) 

2 do color[u] <— 'white 

3 L <— new LinkedListQ 

4 for each u e V 

5 do if ( color [u] = white) 

6 then TopVisit(u) 

7 return L 


TopVisit(u) 

1 color[u] <— gray; // mark u visited 

2 for (each v e Adj[u]) 

3 do if ( color [v] = white) 

4 then TopVisit(v) 

5 Append u to the front of L 


Figure 8.31 shows the linear order obtained by the topological sort of the sequence of putting on a suit. 
The DAG is still the same; it is only that the order in which the vertices of the graph have been laid out is 
special. As a result, all directed edges go from left to right. 



Figure 8.31: Topological sort of the dressing up sequence 


This is a typical example of how DFS used in applications. The running time is 0(V + E). 
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8.4 Strong Components 


We consider an important connectivity problem with digraphs. When diagraphs are used in 
communication and transportation networks, people want to know that their networks are complete. 
Complete in the sense that that it is possible from any location in the network to reach any other location 
in the digraph. 

A digraph is strongly connected if for every pair of vertices u, v e V, u can reach v and vice versa. We 
would like to write an algorithm that determines whether a digraph is strongly connected. In fact, we will 
solve a generalization of this problem, of computing the strongly connected components of a digraph. 

We partition the vertices of the digraph into subsets such that the induced subgraph of each subset is 
strongly connected. We say that two vertices u and v are mutually reachable if u can reach v and vice 
versa. Consider the directed graph in Figure 8.32. The strong components are illustrated in Figure 8.33. 



Figure 8.32: A directed graph 




136 


CHAPTER 8. GRAPHS 



Figure 8.33: Digraph with strong components 


It is easy to see that mutual reachability is an equivalence relation. This equivalence relation partitions 
the vertices into equivalence classes of mutually reachable vertices and these are the strong components. 

If we merge the vertices in each strong component into a single super vertex , and join two super vertices 
(A, B) if and only if there are vertices u G A and v G B such that (u, v) G E, then the resulting digraph is 
called the component digraph. The component digraph is necessarily acyclic. The is illustrated in Figure 
8.34. 



Figure 8.34: Component DAG of super vertices 
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8.4.1 Strong Components and DFS 

Consider DFS of a digraph given in Figure ??. Once you enter a strong component, every vertex in the 
component is reachable. So the DFS does not terminate until all the vertices in the component have been 
visited. Thus all vertices in a strong component must appear in the same tree of the DFS forest. 



Figure 8.35: DFS of a digraph 


fig :dfsofdi graph 

Observe that each strong component is a subtree in the DFS forest. Is it always true for any DFS? The 
answer is “no”. In general, many strong components may appear in the same DFS tree as illustrated in 
Figure 8.36 



Figure 8.36: Another DFS tree of the digraph 
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Is there a way to order the DFS such that it true? Fortunately, the answer is “yes”. Suppose that you knew 
the component DAG in advance. (This is ridiculous, because you would need to know the strong 
components and this is the problem we are trying to solve.) Further, suppose that you computed a 
reversed topological order on the component DAG. That is, for edge (u, v) in the component DAG, then 
v comes before u. This is presented in Figure 8.37. Recall that the component DAG consists of super 
vertices. 



Topological order of component DAG 



Reversed topological order 


Figure 8.37: Reversed topological sort of component DAG 


Now, run DFS, but every time you need a new vertex to start the search from, select the next available 
vertex according to this reverse topological order of the component digraph. Here is an informal 
justification. Clearly once the DFS starts within a given strong component, it must visit every vertex 
within the component (and possibly some others) before finishing. If we do not start in reverse 
topological, then the search may “leak out” into other strong components, and put them in the same DFS 
tree. For example, in the Figure 8.36, when the search is started at vertex ‘a’, not only does it visit its 
component with ‘b’ and ‘c’, but it also visits the other components as well. However, by visiting 
components in reverse topological order of the component tree, each search cannot “leak out” into other 
components, because other components would have already have been visited earlier in the search. 

This leaves us with the intuition that if we could somehow order the DFS, so that it hits the strong 
components according to a reverse topological order, then we would have an easy algorithm for 
computing strong components. However, we do not know what the component DAG looks like. (After 
all, we are trying to solve the strong component problem in the first place). The trick behind the strong 
component algorithm is that we can find an ordering of the vertices that has essentially the necessary 
property, without actually computing the component DAG. 

We will discuss the algorithm without proof. Define G T to be the digraph with the same vertex set at G 
but in which all edges have been reversed in direction. This is shown in Figure 8.38. Given an adjacency 
list for G, it is possible to compute G T in 0(V + E) time. 
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Figure 8.38: The digraph G T 


Observe that the strongly connected components are not affected by reversal of all edges. If u and v are 
mutually reachable in G, then this is certainly true in G T . All that changes is that the component DAG is 
completely reversed. The ordering trick is to order the vertices of G according to their finish times in a 
DFS. Then visit the nodes of G T in decreasing order of finish times. All the steps of the algorithm are 
quite easy to implement, and all operate in 0(V + E) time. Here is the algorithm: 


StrongComponents ( G ) 

1 Run DFS(G) computing finish times f [u] 

2 Compute G T 

3 Sort vertices of G T in decreasing f [u] 

4 Run DFS(G t ) using this order 

5 Each DFS tree is a strong component 


The execution of the algorithm is illustrated in Figures 8.39, 8.40 and 8.41. 
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Figure 8.39: DFS of digraph with vertices in descending order by finish times 



Figure 8.40: Digraph G T and the vertex order for DFS 
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Final DFS with Components 


Figure 8.41: Final DFS with strong components of G T 


The complete proof for why this algorithm works is in CLR. We will discuss the intuition behind why the 
algorithm visits vertices in decreasing order of finish times and why the graph is revered. Recall that the 
main intent is to visit the strong components in a reverse topological order. The problem is how to order 
the vertices so that this is true. Recall from the topological sorting algorithm, that in a DAG, finish times 
occur in reverse topological order (i.e., the first vertex in the topological order is the one with the highest 
finish time). So, if we wanted to visit the components in reverse topological order, this suggests that we 
should visit the vertices in increasing order of finish time, starting with the lowest finishing time. 

This is a good starting idea, but it turns out that it doesn’t work. The reason is that there are many vertices 
in each strong component, and they all have different finish times. For example, in Figure 8.36, observe 
that in the first DFS, the lowest finish time (of 4) is achieved by vertex ‘c\ and its strong component is 
first, not last, in topological order. 

However, there is something to notice about the finish times. If we consider the maximum finish time in 
each component, then these are related to the topological order of the component graph. In fact it is 
possible to prove the following (but we won’t). 


Lemma: Consider a digraph G on which DFS has been run. Label each component with the maximum 
finish time of all the vertices in the component, and sort these in decreasing order. Then this order 
is a topological order for the component digraph. 


For example, in Figure 8.36, the maximum finish times for each component are 18 (for {a, b, c}), 17 (for 
(i d,e }), and 12 (for {fig,h,i}). The order (18, 17, 12) is a valid topological order for the component 
digraph. The problem is that this is not what we wanted. We wanted a reverse topological order for the 
component digraph. So, the final trick is to reverse the digraph. This does not change the component 
graph, but it reverses the topological order, as desired. 
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8.5 Minimum Spanning Trees 


A common problem is communications networks and circuit design is that of connecting together a set of 
nodes by a network of total minimum length. The length is the sum of lengths of connecting wires. 
Consider, for example, laying cable in a city for cable t.v. 

The computational problem is called the minimum spanning tree (MST) problem. Formally, we are given 
a connected, undirected graph G = (V, E) Each edge (u, v) has numeric weight of cost. We define the 
cost of a spanning tree T to be the sum of the costs of edges in the spanning tree 


w(T) = ^ w(u,v) 

(u,v)eT 


A minimum spanning tree is a tree of minimum weight. 

Figures ??, ?? and ?? show three spanning trees for the same graph. The first is a spanning tree but is not 
a MST; the other two are. 



Figure 8.42: A spanning tree Figure 8.43: A minimum Figure 8.44: Another mini- 

that is not MST spanning tree mum spanning tree 

We will present two greedy algorithms (Kruskal’s and Prim’s) for computing MST. Recall that a greedy 
algorithm is one that builds a solution by repeatedly selecting the cheapest among all options at each 
stage. Once the choice is made, it is never undone. 

Before presenting the two algorithms, let us review facts about free trees. A free tree is a tree with no 
vertex designated as the root vertex. A free tree with n vertices has exactly n — 1 edges. There exists a 
unique path between any two vertices of a free tree. Adding any edge to a free tree creates a unique cycle. 
Breaking any edge on this cycle restores the free tree. This is illustrated in Figure 8.45. When the edges 
(b, e) or (b, d) are added to the free tree, the result is a cycle. 
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Figure 8.45: Free tree facts 


8.5.1 Computing MST: Generic Approach 

Let G = (V, E) be an undirected, connected graph whose edges have numeric weights. The intuition 
behind greedy MST algorithm is simple: we maintain a subset of edges E of the graph . Call this subset 
A. Initially, A is empty. We will add edges one at a time until A equals the MST. 

A subset A C E is viable if A is a subset of edges of some MST. An edge (u,v) e E — A is safe if 
A U {(u, v)} is viable. In other words, the choice (u, v) is a safe choice to add so that A can still be 
extended to form a MST. 

Note that if A is viable, it cannot contain a cycle. A generic greedy algorithm operates by repeatedly 
adding any safe edge to the current spanning tree. 

When is an edge safe? Consider the theoretical issues behind determining whether an edge is safe or not. 
Let S be a subset of vertices S C V. A cut (S, V — S) is just a partition of vertices into two disjoint 
subsets. An edge (u, v) crosses the cut if one endpoint is in S and the other is in V — S. 

Given a subset of edges A, a cut respects A if no edge in A crosses the cut. It is not hard to see why 
respecting cuts are important to this problem. If we have computed a partial MST and we wish to know 
which edges can be added that do not induce a cycle in the current MST, any edge that crosses a 
respecting cut is a possible candidate. 
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8.5.2 Greedy MST 


An edge of E is a light edge crossing a cut if among all edges crossing the cut, it has the minimum 
weight. Intuition says that since all the edges that cross a respecting cut do not induce a cycle, then the 
lightest edge crossing a cut is a natural choice. The main theorem which drives both algorithms is the 
following: 


MST Lemma: Let G = (V, E) be a connected, undirected graph with real-valued weights on the edges. 
Let A be a viable subset of E (i.e., a subset of some MST). Let (S, V — S) be any cut that respects 
A and let (u, v) be a light edge crossing the cut. Then the edge (u, v) is safe for A. This is 
illustrated in Figure 8.46. 



Figure 8.46: Subset A with a cut (wavy line) that respects A 


MST Proof: It would simplify the proof if we assume that all edge weights are distinct. Let T be any 
MST for G. If T contains (u, v) then we are done. This is shown in Figure 8.47 where the lightest 
edge (u, v) with cost 4 has been chosen. 
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Figure 8.47: MST T which contains light edge (u, v) 


Suppose no MST contains (u, v). Such a tree is shown in Figure 8.48. We will derive a 
contradiction. 



Figure 8.48: MST T which does not contains light edge (u, v) 


Add (u, v) to T thus creating a cycle as illustrated in Figure 8.49. 





146 


CHAPTER 8. GRAPHS 



Figure 8.49: Cycle created due to T + (u,v) 


Since u and v are on opposite sides of the cut, and any cycle must cross the cut an even number of 
times, there must be at least one other edge (x, y ) in T that crosses the cut. The edge (x, y ) is not in 
A because the cut respects A. By removing (x,y ) we restore a spanning tree, call it T ; . This is 
shown in Figure 8.50 



Figure 8.50: Tree T' = T — (x,y) + (u,v) 


Wehavew(T') = w(T) —w[x,y) +w(u,v). Since (u,v) is the lightest edge crossing the cut we 
have w(u, v) < w(x,y). Thus w(T') < w(T) which contradicts the assumption that T was an 
MST. 
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8.5.3 KruskaPs Algorithm 


Kruskal’s algorithm works by adding edges in increasing order of weight (lightest edge first). If the next 
edge does not induce a cycle among the current set of edges, then it is added to A. If it does, we skip it 
and consider the next in order. As the algorithm runs, the edges in A induce a forest on the vertices. The 
trees of this forest are eventually merged until a single tree forms containing all vertices. 

The tricky part of the algorithm is how to detect whether the addition of an edge will create a cycle in A. 
Suppose the edge being considered has vertices (u, v). We want a fast test that tells us whether u and v 
are in the same tree of A. This can be done using the Union-Find data structure which supports the 
following O(logn) operations: 

Create-set(u): Create a set containing a single item u. 

Find-set(u): Find the set that contains u 

Union(u,v): merge the set containing u and set containing v into a common set. 

In Kruskal’s algorithm, the vertices will be stored in sets. The vertices in each tree of A will be a set. The 
edges in A can be stored as a simple list. Here is the algorithm: Figures 8.51 through ?? demonstrate the 
algorithm applied to a graph. 


kruskal(G = (V, E)) 

1 A <— {} 

2 for ( each u G V) 

3 do create_set(u) 

4 sort E in increasing order by weight w 

5 for ( each (u, v) in sorted edge list) 

6 do if (find(u) ^ find(v)) 

7 then add (u,v) to A 

8 union(u,v) 

9 return A 




148 


CHAPTER 8. GRAPHS 



Figure 8.51: Kruskal algorithm: (b, d) and (d, e) added 



Figure 8.52: Kruskal algorithm: (c, g) and (a, e) added 



Figure 8.53: Kruskal algorithm: unsafe edges 
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Figure 8.54: Kruskal algorithm: (e,f)added 



Figure 8.55: Kruskal algorithm: more unsafe edges and final MST 


Analysis: 

Since the graph is connected, we may assume that E > V — 1 . Sorting edges ( line 4) takes 0(E log E). 
The for loop {line 5) performs 0(E) find and 0(V) union operations. Total time for union — find is 
0(Ea(V)) where a(V) is the inverse Ackerman function. a(V) < 4 for V less the number of atoms in 
the entire universe. Thus the time is dominated by sorting. Overall time for Kruskal is 
0(ElogE) = 0(Elog V) if the graph is sparse. 

8.5.4 Prim’s Algorithm 

Kruskal’s algorithm worked by ordering the edges, and inserting them one by one into the spanning tree, 
taking care never to introduce a cycle. Intuitively Kruskal’s works by merging or splicing two trees 
together, until all the vertices are in the same tree. 

In contrast, Prim’s algorithm builds the MST by adding leaves one at a time to the current tree. We start 
with a root vertex r; it can be any vertex. At any time, the subset of edges A forms a single tree (in 
Kruskal’s, it formed a forest). We look to add a single vertex as a leaf to the tree. 
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Figure 8.56: Prim’s algorithm: a cut of the graph 


Consider the set of vertices S currently part of the tree and its complement (V — S) as shown in Figure 
8.56. We have cut of the graph. Which edge should be added next? The greedy strategy would be to add 
the lightest edge which in the figure is edge to ’u’. Once u is added, Some edges that crossed the cut are 
no longer crossing it and others that were not crossing the cut are as shown in Figure 8.57 



Figure 8.57: Prim’s algorithm: u selected 


We need an efficient way to update the cut and determine the light edge quickly. To do this, we will make 
use of a priority queue. The question is what do we store in the priority queue? It may seem logical that 
edges that cross the cut should be stored since we choose light edges from these. Although possible, there 
is more elegant solution which leads to a simpler algorithm. 
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For each vertex u G (V — S) (not part of the current spanning tree), we associate a key key [u]. The 
key [u] is the weight of the lightest edge going from u to any vertex in S. If there is no edge from u to a 
vertex in S, we set the key value to oo. We also store in pred[u] the end vertex of this edge in S. We will 
also need to know which vertices are in S and which are not. To do this, we will assign a color to each 
vertex. If the color of a vertex is black then it is in S otherwise not. Here is the algorithm: 


Prim((G,w,t)) 


1 

for ( each u G V) 


2 

do key [u] <— oo; pq. insert (u, key [u] ) 


3 

color [u] <— white 


4 

key [r] <— 0; pred [r] <— nil; pq. decrease key 

(r, key [r] ) ; 

5 

while ( pq. nonempty ()) 


6 

douf- pq.extract_min () 


7 

for ( each u G adj [u] ) 


8 

do if ( color [v] == white )and( w (u,v) < 

key [v]) 

9 

then key [v] = w (u,v) 


10 

pq. decrease key (v, key [v]) 


11 

pred [v] = u 


12 

color [u] = black 



Figures 8.58 through 8.60 illustrate the algorithm applied to a graph. The contents of the priority queue 
are shown as the algorithm progresses. The arrows indicate the predecessor pointers and the numeric 
labels in each vertex is its key value. 



Figure 8.58: Prim’s algorithm: edge with weight 4 selected 
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Figure 8.59: Prim’s algorithm: edges with weights 8 and 1 selected 
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Figure 8.60: Prim’s algorithm: the final MST 


Analysis: 

It takes 0(log V) to extract a vertex from the priority queue. For each incident edge, we spend potentially 
0(log V) time decreasing the key of the neighboring vertex. Thus the total time is 
0(log V + deg(u) log V). The other steps of update are constant time. 

So the overall running time is 

T(V,E) = ^( lo S V + de g( u ) lo g V ) 

uev 

= logV y~ (1 + deg(u)) 

uev 

= (log V) (V + 2E) = ©( ( V + E) log V) 
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Since G is connected, V is asymptotically no greater than E so this is @(E log V), same as Kruskal’s 
algorithm. 


8.6 Shortest Paths 

A motorist wishes to find the shortest possible route between Peshawar and Karachi. Given a road map of 
Pakistan on which the distance between each pair of adjacent cities is marked Can the motorist determine 
the shortest route? 

In the shortest-paths problem We are given a weighted, directed graph G = (V, E) The weight of path 
p =< v 0 , V] , . . . , v k > is the sum of the constituent edges: 

k 

W(p) = Y wtVj-T.Vj) 

i=1 

We define the shortest-path weight from u to v by 

f min{w (p ) : u v} if there is a path from u to v 

S(u,v) = ? 

I oo otherwise 

Problems such as shortest route between cities can be solved efficiently by modelling the road map as a 
graph. The nodes or vertices represent cities and edges represent roads. Edge weights can be interpreted 
as distances. Other metrics can also be used, e.g., time, cost, penalties and loss. 

Similar scenarios occur in computer networks like the Internet where data packets have to be routed. The 
vertices are routers. Edges are communication links which may be be wire or wireless. Edge weights can 
be distance, link speed, link capacity link delays, and link utilization. 

The breadth-first-search algorithm we discussed earlier is a shortest-path algorithm that works on 

un- weighted graphs. An un- weighted graph can be considered as a graph in which every edge has weight 

one unit. 

There are a few variants of the shortest path problem. We will cover their definitions and then discuss 
algorithms for some. 

Single-source shortest-path problem: Find shortest paths from a given (single) source vertex s G V to 
every other vertex v e V in the graph G . 

Single-destination shortest-paths problem: Find a shortest path to a given destination vertex t from 

each vertex v. We can reduce the this problem to a single-source problem by reversing the direction 
of each edge in the graph. 

Single-pair shortest-path problem: Find a shortest path from u to v for given vertices u and v. If we 
solve the single-source problem with source vertex u, we solve this problem also. No algorithms 
for this problem are known to run asymptotically faster than the best single-source algorithms in 
the worst case. 
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All-pairs shortest-paths problem: Find a shortest path from u to v for every pair of vertices u and v. 
Although this problem can be solved by running a single-source algorithm once from each vertex, 
it can usually be solved faster. 

8.6.1 Dijkstra’s Algorithm 

Dijkstra’s algorithm is a simple greedy algorithm for computing the single-source shortest-paths to all 
other vertices. Dijkstra’s algorithm works on a weighted directed graph G = (V, E) in which all edge 
weights are non-negative, i.e., w(u, v) > 0 for each edge (u, v) e E. 

Negative edges weights maybe counter to intuition but this can occur in real life problems. However, we 
will not allow negative cycles because then there is no shortest path. If there is a negative cycle between, 
say, s and t, then we can always find a shorter path by going around the cycle one more time. 



Figure 8.61: Negative weight cycle 


The basic structure of Dijkstra’s algorithm is to maintain an estimate of the shortest path from the source 
vertex to each vertex in the graph. Call this estimate d[v]. Intuitively, d[v] will the length of the shortest 
path that the algorithm knows of from s to v. This value will always be greater than or equal to the true 
shortest path distance from s to v. I.e., d[v] > 6(u, v). Initially, we know of no paths, so d[v] = oo. 
Moreover, d[s] = 0 for the source vertex. 

As the algorithm goes on and sees more and more vertices, it attempts to update d[v] for each vertex in 
the graph. The process of updating estimates is called relaxation. Here is how relaxation works. 

Intuitively, if you can see that your solution is not yet reached an optimum value, then push it a little 
closer to the optimum. In particular, if you discover a path from s to v shorter than d[v], then you need to 
update d[v]. This notion is common to many optimization algorithms. 

Consider an edge from a vertex u to v whose weight is w(u, v). Suppose that we have already computed 
current estimates on d[u] and d[v]. We know that there is a path from s to u of weight d[u]. By taking 



8.6. SHORTEST PATHS 


155 


this path and following it with the edge (u,v) we get a path to v of length d[u] + w(u, v). If this path is 
better than the existing path of length d[v] to v, we should take it. The relaxation process is illustrated in 
the following figure. We should also remember that the shortest way back to the source is through u by 
updating the predecessor pointer. 



Figure 8.62: Vertex u relaxed Figure 8.63: Vertex v relaxed 


RELAX( (u, v) ) 

1 if (d[u] +w(u,v) < d[v]) 

2 then d[v] <— d[u] + w(u,v) 

3 pred[v] = u 


Observe that whenever we set d[v] to a finite value, there is always evidence of a path of that length. 
Therefore d[v] > 6(s, v). If d[v] = 6(s,v), then further relaxations cannot change its value. 

It is not hard to see that if we perform Relax(U,V) repeatedly over all edges of the graph, the d[v] 
values will eventually converge to the final true distance value from s. The cleverness of any shortest path 
algorithm is to perform the updates in a judicious manner, so the convergence is as fast as possible. 

Dijkstra’s algorithm is based on the notion of performing repeated relaxations. The algorithm operates by 
maintaining a subset of vertices, S C V, for which we claim we know the true distance, d[v] = 6(s, v). 

Initially S = 0, the empty set. We set d[u] = 0 and all others to oo. One by one we select vertices from 
V — S to add to S. 

How do we select which vertex among the vertices of V — S to add next to S? Here is greediness comes 
in. For each vertex u G ( V — S), we have computed a distance estimate d[u]. 

The greedy thing to do is to take the vertex for which d[u] is minimum, i.e., take the unprocessed vertex 
that is closest by our estimate to s. Later, we justify why this is the proper choice. In order to perform 
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this selection efficiently, we store the vertices of V — S in a priority queue. 


Dijkstra((G,w, s)) 

1 for ( each u e V) 

2 do d[u] <— oo 

3 pq.insert (u, d[u]) 

4 d[s] <— 0; pred [s] <— nil; pq. decrease Jkey (s, d [s] ) ; 

5 while ( pq. nonempty ()) 

6 dou<— pq.extract_min () 

7 for ( each v e adj[u]) 

8 do if (d[u] +w(u,v) < d[v] ) 

9 then d[v] = d[u] + w(u, v) 

10 pq.decrease_key (v, d[v]) 

1 1 pred [v] = u 


Note the similarity with Prim’s algorithm, although a different key is used here. Therefore the running 
time is the same, i.e., 0(Elog V). 

Figures 8.64 through ?? demonstrate the algorithm applied to a directed graph with no negative weight 
edges. 



Figure 8.64: Dijkstra’s algorithm: select 0 
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Figure 8.65: Dijkstra’s algorithm: select 2 



Figure 8.66: Dijkstra’s algorithm: select 5 


select 6 
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select 7 



fig:dijlast 


8.6.2 Correctness of Dijkstra’s Algorithm 

We will prove the correctness of Dijkstr’s algorithm by Induction. We will use the definition that 6(s, v) 
denotes the minimal distance from s to v. 

For the base case 

1. S ={s} 

2. d(s) = 0, which is 6(s, s) 

Assume that d(v) = S(s, v) for every v 6 S and all neighbors of v have been relaxed. If d(u) < d(u') 
for every u'6 V then d(u) = 6(s,u), and we can transfer u from V to S, after which d(v) = 5 ( s , v) for 
every v e S. 

We do this as a proof by contradiction. Suppose that d(u) > 6(s,u). The shortest path from s to u, 
p(s, u), must pass through one or more vertices exterior to S. Let x be the last vertex inside S and y be 
the first vertex outside S on this path to u. Then p(s,u) = p(s, x) U {(x,p)} U p(y,u). 
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Figure 8.69: Correctness of Dijkstra’s algorithm 


The length of p(s,u) is S(s,u) = 6(s,y) + 6(y,u). Because y was relaxed when x was put into S, 
d(y) = 6(s,y) by the convergence property. Thus d(y) < 6(s,u) < d(u). But, because d(u) is the 
smallest among vertices not in S, d(u) < d(y ) . The only possibility is d(u) = d(y ), which requires 
d(u) = 6(s,u) contradicting the assumption. 

By the upper bound property, d(u) > 6(s,u). Since d(u) > 6(s,u) is false, d(u) = 6(s,u), which is 
what we wanted to prove. Thus, if we follow the algorithm’s procedure, transferring from V to S, the 
vertex in V with the smallest value of d(u) then all vertices in S have d(v) = 6(s, v) 


8.6.3 Bellman-Ford Algorithm 

Dijkstra’s single-source shortest path algorithm works if all edges weights are non-negative and there are 
no negative cost cycles. Bellman-Ford allows negative weights edges and no negative cost cycles. The 
algorithm is slower than Dijkstra’s, running in 0(VE) time. 

Like Dijkstra’s algorithm, Bellman-Ford is based on performing repeated relaxations. Bellman-Ford 
applies relaxation to every edge of the graph and repeats this V — 1 times. Here is the algorithm; its is 
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illustrated in Figure 8.70. 


Bellman-Ford(G, w, s) 

1 for ( each u e V) 

2 do d[u] <— oo 

3 pred[u] = nil 

4 

5 d[s] <- 0; 

6 for i = 1 to V — 1 

7 do for ( each (u,v) in E ) 

8 do relax(u, v) 



Figure 8.70: The Bellman-Ford algorithm 


8.6.4 Correctness of Bellman-Ford 

Think of Bellman-Ford as a sort of bubble-sort analog for shortest path. The shortest path information is 
propagated sequentially along each shortest path in the graph. Consider any shortest path from s to some 
other vertex u: (v 0 , Vi , . . . , v k ) where v 0 = s and v k = u. 

Since a shortest path will never visit the same vertex twice, we know that k < V — 1 . Hence the path 
consists of at most V — 1 edges. Since this a shortest path, it is 6 ( s , Vi), the true shortest path cost from s 
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to Vt that satisfies the equation: 


6(s,vt) = 6(s, Vi t ) +w(Vi_i,Vi) 

Claim: We assert that after the i th pass of the “for-i” loop, d[vj = 8(s, v L ). 

Proof: The proof is by induction on i. Observe that after the initialization (pass 0), d[v-|] = d[s] = 0. 

In general, prior to the i th pass through the loop, the induction hypothesis tells us that 

d[Vi i] = 6(s, Vi_i). After the i th pass, we have done relaxation on the edge (Vi_i, Vi) (since we do 

relaxation along all edges). Thus after the i th pass we have 


d[vj < d[vi_i] +w(v i _ 1 ,v i ) 

= 6(s,Vi_i) +w(Vi_i,Vi) 
= 6(s,vt) 


Recall from Dijkstra’s algorithm that d[vj is never less than 5(s, V|). Thus, d[vj is in fact equal to 
6(s, Vi). This completes the induction proof. 

In summary, after i passes through the for loop, all vertices that are i edges away along the shortest path 
tree from the source have the correct values stored in d[u]. Thus, after the (V — 1 ) st iteration of the for 
loop, all vertices u have correct distance values stored in d[u]. 

8.6.5 Floyd- Warshall Algorithm 

We consider the generalization of the shortest path problem: to compute the shortest paths between all 
pairs of vertices. This is called the all-pairs shortest paths problem. 

Let G = (V, E) be a directed graph with edge weights. If (u, v) e E is an edge then w(u, v) denotes its 

weight. 6(u, v) is the distance of the minimum cost path between u and v. We will allow G to have 

negative edges weights but will not allow G to have negative cost cycles. We will present an 0(n 3 ) 
algorithm for the all pairs shortest path. The algorithm is called the Floyd-War shall algorithm and is 
based on dynamic programming . 

We will use an adjacency matrix to represent the digraph. Because the algorithm is matrix based, we will 
employ the common matrix notation, using i, j and k to denote vertices rather than u, v and w. 

The input is an n x n matrix of edge weights: 

f 0 if i = j 

Wij = < w(i, j) if i ^ ) and (i, j) e E 

[oo if i ^ j and (i, j) ^ E 

The output will be an n x n distance matrix D = dy, where = 6(i, j), the shortest path cost from 
vertex i to j . 
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The algorithm dates back to the early 60’s. As with other dynamic programming algorithms, the genius 
of the algorithm is in the clever recursive formulation of the shortest path problem. For a path 
p = (vi,v 2 , . . . , Vi, we say that the vertices v 2 , v 3 , . . . , v L i are the intermediate vertices of this path. 

Formulation: Define djj^ to be the shortest path from i to j such that any intermediate vertices on the 
path are chosen from the set {1 , 2, ... , k}. The path is free to visit any subset of these vertices and in any 
order. How do we compute d-^ 1 assuming we already have the previous matrix d !k 1 '? There are two 
basic cases: 

1. Don’t go through vertex k at all. 

2. Do go through vertex k. 



Don’t go through k at all 

Then the shortest path from i to j uses only intermediate vertices {1 , 2, . . . , k — 1 }. Hence the length of 
the shortest is d-|" 1 1 

Do go through k 

First observe that a shortest path does not go through the same vertex twice, so we can assume that we 
pass through k exactly once. That is, we go from i to k and then from k to j. In order for the overall path 
to be as short as possible, we should take the shortest path from i to k and the shortest path from k to j. 
Since each of these paths uses intermediate vertices {1 , 2, . . . , k — 1 }, the length of the path is 

d (k-p d (k-t) 
u ik ' u kj 

The following illustrate the process in which the value of d^ 2 is updated as k goes from 0 to 4. 
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Figure 8.76: k = 4, dj, 4 ] = 7 (3 -> 1 -> 4 -> 2) 

This suggests the following recursive (DP) formulation: 

< = 

d; ) k) = mm(d;r l Ur 5 + d^ 1 >) 

The final answer is djj 1 * because this allows all possible vertices as intermediate vertices. 

As is the case with DP algorithms, we will avoid recursive evaluation by generating a table for d-^. The 
algorithm also includes mid-vertex pointers stored in mid[i, j] for extracting the final path. 


Floyd-Warshall(u, w[1 ..n, 1 ..n]) 

1 for (i = 1 , n) 

2 do for (j = 1 , n) 

3 do d[i, j] w[i, j]; mid[i, j] <— null 

4 for (k = 1 ,n) 

5 do for (i = 1 , n) 

6 do for (j = 1 , n) 

7 do if ( d [i, k] T- d[k, j] < d[i, j]) 

8 then d[i, j] = d[i, k] + d[k, j] 

9 mid[i, j] = k 


Clearly, the running time is 0(n 3 ). The space used by the algorithm is 0 (n 2 ) . 

Figure 8.77 through 8.81 demonstrate the algorithm when applied to a graph. The matrix to left of the 
graph contains the matrix d entries. A circle around an entry dt^ indicates that it was updated in the 
current k iteration. 
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Figure 8.77: Floyd- Warshall Algorithm example: d ( 0) 



Figure 8.78: Floyd- Warshall Algorithm example: d ( l ) 
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Figure 8.79: Floyd- Warshall Algorithm example: d ( 2) 



Figure 8.80: Floyd- Warshall Algorithm example: d ( 3) 
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Figure 8.81: Floyd- Warshall Algorithm example: d ( 4) 


Extracting Shortest Path: 

The matrix d holds the final shortest distance between pairs of vertices. In order to compute the shortest 
path, the mid-vertex pointers mid[i, j] can be used to extract the final path. Whenever we discovered that 
the shortest path from i to j passed through vertex k, we set mid[i, j] = k. If the shortest path did not 
pass through k then mid[i, j] = null. 

To find the shortest path from i to j, we consult mid[i, j]. If it is null, then the shortest path is just the 
edge (i, j). Otherwise we recursively compute the shortest path from i to mid[i, j] and the shortest path 
from mid[i, j] to j. 


PATH(i,j) 

1 if (mid[i, j] == null) 

2 then output(i, j) 

3 else PATH(i, mid[i, j]) 

4 PATH(mld[i, j], j) 
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Chapter 9 

Complexity Theory 


So far in the course, we have been building up a “bag of tricks” for solving algorithmic problems. 
Hopefully you have a better idea of how to go about solving such problems. What sort of design 
paradigm should be used: divide-and-conquer, greedy, dynamic programming. 

What sort of data structures might be relevant: trees, heaps, graphs. What is the running time of the 
algorithm. All of this is fine if it helps you discover an acceptably efficient algorithm to solve your 
problem. 

The question that often arises in practice is that you have tried every trick in the book and nothing seems 
to work. Although your algorithm can solve small problems reasonably efficiently (e.g., n < 20), for the 
really large problems you want to solve, your algorithm never terminates. When you analyze its running 
time, you realize that it is running in exponential time, perhaps n^, or 2 n , or 2 2n , or n! or worse!. 

By the end of 60’s, there was great success in finding efficient solutions to many combinatorial problems. 
But there was also a growing list of problems for which there seemed to be no known efficient 
algorithmic solutions. 

People began to wonder whether there was some unknown paradigm that would lead to a solution to 
there problems. Or perhaps some proof that these problems are inherently hard to solve and no 
algorithmic solutions exist that run under exponential time. 

Near the end of the 1960’s, a remarkable discovery was made. Many of these hard problems were 
interrelated in the sense that if you could solve any one of them in polynomial time, then you could solve 
all of them in polynomial time, this discovery gave rise to the notion of NP-completeness. 

This area is a radical departure from what we have been doing because the emphasis will change. The 
goal is no longer to prove that a problem can be solved efficiently by presenting an algorithm for it. 
Instead we will be trying to show that a problem cannot be solved efficiently. 

Up until now all algorithms we have seen had the property that their worst-case running time are bounded 
above by some polynomial in n. A polynomial time algorithm is any algorithm that runs in Ofn k ) time. 
A problem is solvable in polynomial time if there is a polynomial time algorithm for it. 

Some functions that do not look like polynomials (such as 0(nlog rt) are bounded above by polynomials 
(such as 0(n 2 )). Some functions that do look like polynomials are not. For example, suppose you have 
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an algorithm that takes as input a graph of size n and an integer k and run in 0(n k ) time. 

Is this a polynomial time algorithm? No, because k is an input to the problem so the user is allowed to 
choose k = rt, implying that the running time would be 0(ri n ). 0(n n ) is surely not a polynomial in n. 
The important aspect is that the exponent must be a constant independent of n. 


9.1 Decision Problems 

Most of the problems we have discussed involve optimization of one form of another. Find the shortest 
path, find the minimum cost spanning tree, maximize the knapsack value. For rather technical reasons, 
the NP-complete problems we will discuss will be phrased as decision problems. 

A problem is called a decision problem if its output is a simple “yes” or “no” (or you may this of this as 
true/false, 0/1, accept/reject.) We will phrase may optimization problems as decision problems. For 
example, the MST decision problem would be: Given a weighted graph G and an integer k, does G have 
a spanning tree whose weight is at most k? 

This may seem like a less interesting formulation of the problem. It does not ask for the weight of the 
minimum spanning tree, and it does not even ask for the edges of the spanning tree that achieves this 
weight. However, our job will be to show that certain problems cannot be solved efficiently. If we show 
that the simple decision problem cannot be solved ejciently, then the more general optimization problem 
certainly cannot be solved efficiently either. 


9.2 Complexity Classes 

Before giving all the technical definitions, let us say a bit about what the general classes look like at an 

intuitive level. 

Class P: This is the set of all decision problems that can be solved in polynomial time. We will 
generally refer to these problems as being “easy” or “efficiently solvable”. 

Class NP: This is the set of all decision problems that can be verified in polynomial time. This class 

contains P as a subset. It also contains a number of problems that are believed to be very “ hard” to 
solve. 

Class NP: The term “NP” does not mean “not polynomial”. Originally, the term meant “ 

non-deterministic polynomial” but it is a bit more intuitive to explain the concept from the 
perspective of verification. 

Class NP-hard: In spite of its name, to say that a problem is NP-hard does not mean that it is hard to 

solve. Rather, it means that if we could solve this problem in polynomial time, then we could solve 
all NP problems in polynomial time. Note that for a problem to NP-hard, it does not have to be in 
the class NP. 
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Class NP-complete: A problem is NP-complete if (1) it is in NP and (2) it is NP-hard. 


The Figure 9.1 illustrates one way that the sets P, NP, NP-hard, and NP-complete (NPC) might look. We 
say might because we do not know whether all of these complexity classes are distinct or whether they 
are all solvable in polynomial time. The Graph Isomorphism, which asks whether two graphs are 
identical up to a renaming of their vertices. It is known that this problem is in NP, but it is not known to 
be in P. The other is QBF, which stands for Quantified Boolean Formulas. In this problem you are given a 
boolean formula and you want to know whether the formula is true or false. 


O' 

O- 


NP-hard 


NP o 




Quantified Boolean 
Formulas 


No Hamiltonian cycle 



0/1 knapsack 
Hamiltonian cycle 
Satisfiability 


Graph Isomorphism 


MST 

Strong connectivity 


Figure 9.1: Complexity Classes 


9.3 Polynomial Time Verification 


Before talking about the class of NP-complete problems, it is important to introduce the notion of a 
verification algorithm. Many problems are hard to solve but they have the property that it easy to verify 
the solution if one is provided. Consider the Hamiltonian cycle problem. 

Given an undirected graph G, does G have a cycle that visits every vertex exactly once? There is no 
known polynomial time algorithm for this problem. 
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Figure 9.2: Hamiltonian Cycle 


However, suppose that a graph did have a Hamiltonian cycle. It would be easy for someone to convince 
of this. They would simply say: “the cycle is (v 3 ,v 7 ,vi, . . . ,Vi3) We could then inspect the graph and 
check that this is indeed a legal cycle and that it visits all of the vertices of the graph exactly once. Thus, 
even though we know of no efficient way to solve the Hamiltonian cycle problem, there is a very efficient 
way to verify that a a given cycle is indeed a Hamiltonian cycle. 

The piece of information that allows verification is called a certificate. Note that not all problems have 
the property that they are easy to verify. For example, consider the following two: 

1. UHC = {(G) |G has a unique Hamiltonian cycle} 

2. HC = {(G)|G has no Hamiltonian cycle} 

Suppose that a graph G is in UHC. What information would someone give us that would allow us to 
verify this? They could give us an example of the unique Hamiltonian cycle and we could verify that it is 
a Hamiltonian cycle. But what sort of certificate could they give us to convince us that this is the only 
one? 

They could give another cycle that is not Hamiltonian. But this does not mean that there is not another 
cycle somewhere that is Hamiltonian. They could try to list every other cycle of length n, but this is not 
efficient at all since there are n! possible cycles in general. Thus it is hard to imagine that someone could 
give us some information that would allow us to efficiently verify that the graph is in UHC. 


9.4 The Class NP 

The class NP is a set of all problems that can be verified by a polynomial time algorithm. Why is the set 
called “NP” and not “VP”? The original term NP stood for non-deterministic polynomial time. This 
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referred to a program running on a non-deterministic computer that can make guesses. Such a computer 
could non-deterministically guess the value of the certificate, and then verify it in polynomial time. We 
have avoided introducing non-determinism here; it is covered in other courses such as automata or 
complexity theory. 

Observe that P C NP. In other words, if we can solve a problem in polynomial time, we can certainly 
verify the solution in polynomial time. More formally, we do not need to see a certificate to solve the 
problem; we can solve it in polynomial time anyway. 

However, it is not known whether P = NP. It seems unreasonable to think that this should be so. Being 
able to verify that you have a correct solution does not help you in finding the actual solution. The belief 
is that P 7 ^ N P but no one has a proof for this. 


9.5 Reductions 


The class NP-complete (NPC) problems consists of a set of decision problems (a subset of class NP) that 
no one knows how to solve efficiently. But if there were a polynomial solution for even a single 
NP-complete problem, then ever problem in NPC will be solvable in polynomial time. For this, we need 
the concept of reductions. 

Consider the question: Suppose there are two problems, A and B. You know (or you strongly believe at 
least) that it is impossible to solve problem A in polynomial time. You want to prove that B cannot be 
solved in polynomial time. We want to show that 

(A^P)^(B^P) 


How would you do this? Consider an example to illustrate reduction: The following problem is 
well-known to be NPC: 


3-color: Given a graph G, can each of its vertices be labelled with one of 3 different colors such that two 
adjacent vertices have the same label (color). 


Coloring arises in various partitioning problems where there is a constraint that two objects cannot be 
assigned to the same set of partitions. The term “coloring” comes from the original application which 
was in map drawing. Two countries that share a common border should be colored with different colors 

It is well known that planar graphs can be colored (maps) with four colors. There exists a polynomial 
time algorithm for this. But determining whether this can be done with 3 colors is hard and there is no 
polynomial time algorithm for it. In Figure 9.3, the graph on the left can be colored with 3 colors while 
the graph on the right cannot be colored. 
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Figure 9.3: Examples of 3-colorable and non-3 -colorable graphs 


Examplel: Fish tank problem 

Consider the following problem than can be solved with the graph coloring approach. A tropical fish 
hobbyist has six different types of fish designated by A, B, C, D, E, and F, respectively. Because of 
predator-prey relationships, water conditions and size, some fish can be kept in the same tank. The 
following table shows which fish cannot be together: 


Type 

Cannot be with 

A 

B, C 

B 

A ,C, E 

C 

A, B, D, E 

D 

C, F 

E 

B, C, F 

F 

D, E 


These constraints can be displayed as a graph where an edge between two vertices exists if the two 
species cannot be together. This is shown in Figure 9.4. For example, A cannot be with B and C; there is 
an edge between A and B and between A and C. 

Given these constraints, What is the smallest number of tanks needed to keep all the fish? The answer 
can be found by coloring the vertices in the graph such that no two adjacent vertices have the same color. 
This particular graph is 3-colorable and therefore, 3 fish tanks are enough. This is depicted in Figure 9.5. 
The 3 fish tanks will hold fish as follows: 


Tank 1 

Tank 2 

Tank 3 

A, D 

F,C 

B, E 
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Figure 9.4: Graph representing constraints be- Figure 9.5: Fish tank graph colored with 3 colors 
tween fish species 

The 3-color (3Col) problem will the play the role of A, which we strongly suspect to not be solvable in 
polynomial time. For our problem B, consider the following problem: Given a graph G = (V, E), we say 
that a subset of vertices V'CV forms a clique if for every pair of vertices u, v e V', the edge 
(u, v) e V' That is, the subgraph induced by V' is a complete graph. 

Clique Cover: Given a graph G and an integer k, can we find k subsets of vertices Vi , V 2 , . . . , V k , such 
that (J t Vi = V, and that each Vt is a clique of G . 

The following figure shows a graph that has a clique cover of size 3. There are three subgraphs that are 
complete. 



Figure 9.6: Graph with clique cover of size 3 
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The clique cover problem arises in applications of clustering. We put an edge between two nodes if they 
are similar enough to be clustered in the same group. We want to know whether it is possible to cluster all 
the vertices into k groups. 

Suppose that you want to solve the CCov problem. But after a while of fruitless effort, you still cannot 
find a polynomial time algorithm for the CCov problem. How can you prove that CCov is likely to not 
have a polynomial time solution? 

You know that 3Col is NP-complete and hence, experts believe that 3Col ^ P. You feel that there is some 
connection between the CCov problem and the 3Col problem. Thus, you want to show that 


(3Col £ P) =» (CCov £ P) 


Both problems involve partitioning the vertices into groups. In the clique cover problem, for two vertices 
to be in the same group, they must be adjacent to each other. In the 3-coloring problem, for two vertices 
to be in the same color group, they must not be adjacent. In some sense, the problems are almost the 
same but the adjacency requirements are exactly reversed. 

We claim that we can reduce the 3-coloring problem into the clique cover problem as follows: Given a 
graph G for which we want to determine its 3-colorability, output the pair (G, 3) where G denotes the 
complement of G. Feed the pair (G, 3) into a routine for clique cover. 

For example, the graph G in Figure 9.7 is 3-colorable and its complement (G, 3) is coverable by 3 
cliques. The graph G in Figure 9.8 is not 3-colorable; it is also not coverable by cliques. 



Figure 9.7: 3-colorable G and clique coverable (G, 3) 
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Figure 9.8: G is not 3 -colorable and (G,3) is not clique coverable 


9.6 Polynomial Time Reduction 

Definition: We say that a decision problem L] is polynomial-time reducible to decision problem L 2 
(written Li < p L 2 ) if there is polynomial time computable function f such that for all x, x G Li if and 
only if f (x) G L 2 . 

In the previous example we showed that 


3Col < P CCov 

In particular, we have f(G) = (G, 3). It is easy to complement a graph in 0(rt 2 ) (i.e., polynomial time). 
For example, flip the 0’s and l’s in the adjacency matrix. 

The way this is used in NP-completeness is that we have strong evidence that Li is not solvable in 
polynomial time. Hence, the reduction is effectively equivalent to saying that “since Li is not likely to be 
solvable in polynomial time, then L 2 is also not likely to be solvable in polynomial time. 


9.7 NP-Completeness 

The set of NP-complete problems is all problems in the complexity class NP for which it is known that if 
any one is solvable in polynomial time, then they all are. Conversely, if any one is not solvable in 
polynomial time, then none are. 

Definition: A decision problem L is NP-Hard if 


L' < P L for all L' G NP. 
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Definition: L is NP-complete if 

1. L G NP and 

2. L ; <p L for some known NP-complete problem L'. 

Given this formal definition, the complexity classes are: 

P: is the set of decision problems that are solvable in polynomial time. 

NP: is the set of decision problems that can be verified in polynomial time. 

NP-Hard: L is NP-hard if for all L' G NP, L' <p L. Thus if we could solve L in polynomial time, we 
could solve all NP problems in polynomial time. 

NP-Complete L is NP-complete if 

1. L G NP and 

2. L is NP-hard. 

The importance of NP-complete problems should now be clear. If any NP-complete problem is solvable 
in polynomial time, then every NP-complete problem is also solvable in polynomial time. Conversely, if 
we can prove that any NP-complete problem cannot be solved in polynomial time, the every 
NP-complete problem cannot be solvable in polynomial time. 


9.8 Boolean Satisfiability Problem: Cook’s Theorem 

We need to have at least one NP-complete problem to start the ball rolling. Stephen Cook showed that 
such a problem existed. He proved that the boolean satisfiability problem is NP-complete. A boolean 
formula is a logical formulation which consists of variables X|. These variables appear in a logical 
expression using logical operations 

1. negation of x: x 

2. boolean or: (x V p ) 

3. boolean and: (xAy) 

For a problem to be in NP, it must have an efficient verification procedure. Thus virtually all NP 
problems can be stated in the form, “does there exists X such that P(X)”, where X is some structure (e.g. 
a set, a path, a partition, an assignment, etc.) and P(X) is some property that X must satisfy (e.g. the set 
of objects must fill the knapsack, or the path must visit every vertex, or you may use at most k colors and 
no two adjacent vertices can have the same color). In showing that such a problem is in NP, the certificate 
consists of giving X, and the verification involves testing that P(X) holds. 
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In general, any set X can be described by choosing a set of objects, which in turn can be described as 
choosing the values of some boolean variables. Similarly, the property P(X) that you need to satisfy, can 
be described as a boolean formula. Stephen Cook was looking for the most general possible property he 
could, since this should represent the hardest problem in NP to solve. He reasoned that computers (which 
represent the most general type of computational devices known) could be described entirely in terms of 
boolean circuits, and hence in terms of boolean formulas. If any problem were hard to solve, it would be 
one in which X is an assignment of boolean values (true/false, 0/1) and P(X) could be any boolean 
formula. This suggests the following problem, called the boolean satisfiability problem. 

SAT: Given a boolean formula, is there some way to assign truth values (0/1, true/false) to the variables 
of the formula, so that the formula evaluates to true? 

A boolean formula is a logical formula which consists of variables Xt , and the logical operations x 
meaning the negation of x, boolean-or (x V y ) and boolean-and (x A p ). Given a boolean formula, we 
say that it is satisfiable if there is a way to assign truth values (0 or 1) to the variables such that the final 
result is 1. (As opposed to the case where no matter how you assign truth values the result is always 0.) 
For example 

(*i A (x 2 VxD) A ((x^AxI) Vxf) 

is satisfiable, by the assignment x ] = 1 , x 2 — 0 and x 3 = 0. On the other hand, 

(xf V (x 2 A x 3 ) ) A (xi V (x^Axj)) A (x 2 Vx 3 ) A (xi Vxl) 

is not satisfiable. Such a boolean formula can be represented by a logical circuit made up of OR, AND 
and NOT gates. For example, Figure 9.9 shows the circuit for the boolean formula 

( (xi A x 4 ) V x 2 ) A ( (x 3 A xfi] V xi) A xjj 



Figure 9.9: Logical circuit for a boolean formula 
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Cook’s Theorem: SAT is NP-complete 

We will not prove the theorem; it is quite complicated. In fact, it turns out that a even more restricted 
version of the satisfiability problem is NP-complete. 

A literal is a variable x or its negation x. A boolean formula is in 3-Conjunctive Normal Form (3-CNF) 
if it is the boolean-and of clauses where each clause is the boolean-or of exactly three literals. For 
example, 

(xi V X 2 V X 3 ) A (xf V x 3 V x 4 ) A (x 2 V X 3 V x^) 

is in 3-CNF form. 3SAT is the problem of determining whether a formula is 3-CNF is satisfiable. 3SAT 
is NP-complete. We can use this fact to prove that other problems are NP-complete. We will do this with 
the independent set problem. 

Independent Set Problem: Given an undirected graph G = (V, E) and an integer k, does G contain a 
subset V' of k vertices such that no two vertices in V' are adjacent to each other. 



Figure 9.10: 


The independent set problem arises when there is some sort of selection problem where there are mutual 
restrictions pairs that cannot both be selected. For example, a company dinner where an employee and 
his or her immediate supervisor cannot both be invited. 

Claim: IS is NP-complete 

The proof involves two parts. First, we need to show that IS € NP. The certificate consists of k vertices 
of V'. We simply verify that for each pair of vertices u, v £ V', there is no edge between them. Clearly, 
this can be done in polynomial time, by an inspection of the adjacency matrix. 

Second, we need to establish that IS is NP-hard This can be done by showing that some known 
NP-compete (3SAT) is polynomial-time reducible to IS. That is, 3SAT <p IS. 
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An important aspect to reductions is that we do not attempt to solve the satisfiability problem. Remember: 
It is NP-complete, and there is not likely to be any polynomial time solution. The idea is to translate the 
similar elements of the satisfiable problem to corresponding elements of the independent set problem. 

What is to be selected? 

3SAT: Which variables are to be assigned the value true, or equivalently, which literals will be true and 
which will be false. 

IS: Which vertices will be placed in V'. 

Requirements: 

3SAT: Each clause must contain at least one true valued literal. 

IS: V' must contain at least k vertices. 

Restrictions: 

3SAT: If Xi is assigned true, then Xi must be false and vice versa. 

IS: If u is selected to be in V' and v is a neighbor of u then v cannot be in V'. 

We want a function which given any 3-CNF boolean formula F , converts it into a pair (G, k) such that 
the above elements are translated properly. Our strategy will be to turn each literal into a vertex. The 
vertices will be in clause clusters of three, one for each clause. Selecting a true literal from some clause 
will correspond to selecting a vertex to add to V'. We will set k equal to the number of clauses, to force 
the independent set subroutine to select one true literal from each clause. To keep the IS subroutine from 
selecting two literals from one clause and none from some other, we will connect all the vertices in each 
clause cluster with edges. To keep the IS subroutine from selecting a literal and its complement to be 
true, we will put an edge between each literal and its complement. 

A formal description of the reduction is given below. The input is a boolean formula F in 3-CNF, and the 
output is a graph G and integer k. 

3SAT-TO-IS(F) 

1 k <— number of clauses in F 

2 for ( each clause C in F ) 

3 do create a clause cluster of 

4 3 vertices from literals of C 

5 for ( each clause cluster (xi , Xz, X 3 ) ) 

6 do create an edge (xt, Xj) between 

7 all pairs of vertices in the cluster 

8 for ( each vertex Xt ) 

9 do create an edge between X| and 
10 all its complement vertices x[ 
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1 1 return ( G , k) // output is graph G and integer k 


If F has k clauses, then G has exactly 3k vertices. Given any reasonable encoding of F , it is an easy 
programming exercise to create G (say as an adjacency matrix) in polynomial time. We claim that F is 
satisfiable if and only if G has an independent set of size k. 

Example: Suppose that we are given the 3-CNF formula: 


(xi V X 2 V x 3 ) A (xi V x 2 V x 3 ) A (xi V x 2 V x 3 ) A (xi V x 2 V x 3 ) 


The following series of figures show the reduction which produces the graph and sets k = 4. First, each of 
the four literals is converted into a three- vertices graph. This is shown in Figure 9.11 



Figure 9.1 1: Four graphs, one for each of the 3-terms literal 


Next, each term is connected to its complement. This is shown in Figure 9.12. 
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Figure 9.12: Each term is connected to its complement 


The boolean formula is satisfied by the assignment *i = 1 , x 2 = 1 x 3 = 0. This implies that the first 
literal of the first and last clauses are 1, the second literal of the second clause is 1, and the third literal of 
the third clause is 1. By selecting vertices corresponding to true literals in each clause, we get an 
independent set of size k = 4. This is shown in Figure 9.13. 



Figure 9.13: Independent set corresponding to Xi = 1 , x 2 = 1 x 3 = 0 


Correctness Proof: 

We claim that F is satisfiable if and only if G has an independent set of size k. If F is satisfiable, then each 
of the k clauses of F must have at least one true literal. Fet V' denote the corresponding vertices from 
each of the clause clusters (one from each cluster). Because we take vertices from each cluster, there are 
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no inter-cluster edges between them, and because we cannot set a variable and its complement to both be 
true, there can be no edge of the form (xt, xf) between the vertices of V' . Thus, V' is an independent set 
of size k. 

Conversely, G has an independent set V' of size k. First observe that we must select a vertex from each 
clause cluster, because there are k clusters, and we cannot take two vertices from the same cluster 
(because they are all interconnected). Consider the assignment in which we set all of these literals to 1. 
This assignment is logically consistent, because we cannot have two vertices labelled xt and x[ in the 
same cluster. 

Finally the transformation clearly runs in polynomial time. This completes the NP-completeness proof. 

Also observe that the the reduction had no knowledge of the solution to either problem. Computing the 
solution to either will require exponential time. Instead, the reduction simple translated the input from 
one problem into an equivalent input to the other problem. 


9.9 Coping with NP- Completeness 

With NP-completeness we have seen that there are many important optimization problems that are likely 
to be quite hard to solve exactly. Since these are important problems, we cannot simply give up at this 
point, since people do need solutions to these problems. Here are some strategies that are used to cope 
with NP-completeness: 

Use brute-force search: Even on the fastest parallel computers this approach is viable only for the 
smallest instance of these problems. 

Heuristics: A heuristic is a strategy for producing a valid solution but there are no guarantees how close 
it to optimal. This is worthwhile if all else fails. 


General search methods: Powerful techniques for solving general combinatorial optimization 
problems. Branch-and-bound, A*-search, simulated annealing, and genetic algorithms 


Approximation algorithm: This is an algorithm that runs in polynomial time (ideally) and produces a 
solution that is within a guaranteed factor of the optimal solution. 



